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1.1. 
Introducción 


La Mecánica de Fluidos estudia las leyes del movimiento de los ﬂuidos y sus procesos de interacción con los cuerpos 
sólidos. La Mecánica de Fluidos como hoy la conocemos es una mezcla de teoría y experimento que proviene por un 
lado de los trabajos iniciales de los ingenieros hidráulicos, de carácter fundamentalmente empírico, y por el otro del 
trabajo de básicamente matemáticos, que abordaban el problema desde un enfoque analítico. Al integrar en una única 
disciplina las experiencias de ambos colectivos, se evita la falta de generalidad derivada de un enfoque estrictamente 
empírico, válido únicamente para cada caso concreto, y al mismo tiempo se permite que los desarrollos analíticos 
matemáticos aprovechen adecuadamente la información experimental y eviten basarse en simpliﬁcaciones artiﬁciales 
alejadas de la realidad. 


La característica fundamental de los ﬂuidos es la denominada ﬂuidez. Un ﬂuido cambia de forma de manera continua 
cuando está sometido a un esfuerzo cortante, por muy pequeño que sea éste, es decir, un ﬂuido no es capaz de soportar 
un esfuerzo cortante sin moverse durante ningún intervalo de tiempo. Unos líquidos se moverán más lentamente que 
otros, pero ante un esfuerzo cortante se moverán siempre. La medida de la facilidad con que se mueve vendrá dada 
por la viscosidad que se trata más adelante, relacionada con la acción de fuerzas de rozamiento. Por el contrario en 
un sólido se produce un cambio ﬁjo γ para cada valor de la fuerza cortante aplicada. En realidad algunos sólidos 
pueden presentar en cierto modo ambos comportamientos, cuando la tensión aplicada está por debajo de un cierto 


⃗F 
⃗F 


γ 


Sólido 
Líquido 


Figura 1–1 Diferencias entre el comportamiento de un líquido y de un sólido frente a una fuerza cortante aplicada. 
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umbral presenta el comportamiento habitual, mientras que por encima de un cierto umbral el sólido puede plastiﬁcar, 
produciéndose una deformación más continua para una fuerza ﬁja, de forma parecida a como ocurre en un ﬂuido. Ésto 
es precisamente lo que ocurre en la zona de ﬂuencia. Si la fuerza persiste, se llega a la rotura del sólido. 


Así, mientras que un sólido experimenta un desplazamiento deﬁnido (o se rompe por completo) bajo la acción de 
una fuerza cortante, en los ﬂuidos pequeñas fuerzas producen grandes deformaciones no elásticas (en general no se 
recupera la forma) a volumen constante, que se realizan de forma continua. Mientras que para un sólido bajo una 
fuerza cortante constante se alcanza un ángulo de deformación determinado y constante, en un ﬂuido debemos hablar 
de una velocidad de deformación constante o no, ya que la deformación se produce de forma continua. 


Dentro de los ﬂuidos, la principal diferencia entre líquidos y gases estriba en las distintas compresibilidades de los 
mismos. 


Gases. Los gases presentan una gran compresibilidad, que inﬂuye sobre las características del ﬂujo, ya que tanto el 
volumen como la densidad varían con facilidad. En el caso de los gases el movimiento térmico vence a las 
fuerzas atractivas y, por tanto tienden a ocupar todo el volumen del recipiente que los contiene. 


Líquidos. En el caso de los líquidos, por el contrario, la compresibilidad es muy débil. Ésto es debido a que las 
fuerzas atractivas entre las moléculas del líquido vencen al movimiento térmico de las mismas, colapsando las 
moléculas y formando el líquido. Al contrario que en el caso de los gases, que tendían a ocupar todo el volumen 
que los contiene, los líquidos tienden a formar una superﬁcie libre. 


La noción de compresibilidad dada es la correspondiente a la estática de ﬂuidos. En dinámica de ﬂuidos, hay casos 
en los que la densidad no varía a lo largo del ﬂujo, incluso en un ﬂuido compresible, por lo que a ese ﬂujo le pode- 
mos aplicar las leyes de los ﬂuidos incompresibles. Tiene en este caso más sentido hablar de ﬂujo compresible o 
incompresible. 


Los ﬂuidos no conservan la forma. Al situarlos en un recipiente toman la forma del mismo (si lo llenan) o de parte 
del mismo. 


1.2. 
Descripción de un ﬂuido. Hipótesis del continuo. 


Para la descripción del movimiento de un ﬂuido recurriremos a las leyes generales de la Mecánica (leyes de Newton, 
leyes de conservación de la cantidad de movimiento y de la energía), junto con relaciones especíﬁcas condicionadas 
por la ﬂuidez. 


A escala microscópica la materia, y en particular un ﬂuido está compuesta de moléculas a cierta distancia promedio 
con espacio vacío entre ellas. Estas moléculas están continuamente moviéndose y colisionando entre sí. Un análisis 
exacto del problema debiera tener en cuenta la acción de cada molécula o grupo de moléculas en el ﬂuido. Este proce- 
dimiento, con algunas simpliﬁcaciones importantes es el que se adopta en Teoría Cinética y en Mecánica Estadística, 
pero es aún demasiado complejo para utilizarlo en el trabajo diario de hidráulica. 


En la mayor parte de los cálculos hidráulicos, el interés está realmente centrado en manifestaciones macroscópicas 
promedio que resultan de la acción conjunta de una gran cantidad de moléculas, manifestaciones como la densidad, 
la presión o la temperatura. En la práctica es posible hacer una simpliﬁcación importante, suponer que todas estas 
manifestaciones son el resultado de la acción de una hipotética distribución continua de materia, a la que denomi- 
naremos el continuo, o el medio continuo, en lugar de estudiar el conglomerado real de las moléculas discretas, de 
mucha mayor complejidad. De esta forma a la hora de estudiar nuestros problemas sustituiremos la materia real por 
este medio continuo ﬁcticio, cuyas propiedades varían de forma continua y reﬂejan las propiedades macroscópicas del 
medio real. 


Este concepto del medio continuo permite una gran simpliﬁcación en el análisis. Por supuesto, este enfoque debe 
utilizarse únicamente cuando arroje resultados razonablemente correctos. Por ejemplo, no puede utilizarse cuando el 
recorrido libre medio de las moléculas es del orden de las magnitudes características del problema. En estas condicio- 
nes, la acción de cada molécula individual es signiﬁcativa y debe estudiarse individualmente. 
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Por ejemplo, consideremos la acción sobre una superﬁcie de la pared en el caso de un depósito cerrado que contiene 
un gas a una cierta presión, en un estado estacionario. Incluso a baja presión, la gran cantidad de colisiones de molé- 
culas sobre la superﬁcie da lugar a una fuerza global que en la práctica puede considerarse independiente del tiempo, 
comportamiento que será correctamente simulado por nuestro hipotético medio continuo. Ahora bien, si la presión 
fuera tan baja que únicamente quedaran en el tanque unas pocas moléculas de forma que el recorrido libre medio de 
las mismas es del orden de magnitud del elemento considerado, se observará una actividad errática según las moléculas 
individuales o los grupos de moléculas bombardean la superﬁcie y no se podrá hablar de una fuerza constante, sino de 
una serie de choques aleatorios contra la superﬁcie. Este comportamiento no podría ser reﬂejado por nuestro medio 
continuo. Lo mismo ocurriría si considerando el gas discreto real, tomamos una superﬁcie muy pequeña, de forma que 
su tamaño es del orden del recorrido libre medio de las moléculas. 


Sin embargo, si ya estamos trabajando con el medio continuo y con magnitudes macroscópicas, un elemento de 
volumen inﬁnitesimal será un elemento de volumen del medio continuo y no de la materia real discreta, con lo que 
trabajaremos con la conﬁanza de que a partir del mismo obtendremos las magnitudes macroscópicas. 


Hipótesis del continuo: 


La materia y las propiedades físicas asociadas a la misma se consideran dispersas de forma continua 
en ella, y no concentradas en pequeñas fracciones (átomos y moléculas) de la misma. 


De este modo, sustituimos la materia real de carácter discreto por una materia ﬁcticia continua cuyas propiedades 
en cada punto vienen dadas por las propiedades promedio de la materia real en el entorno de ese punto. 


1.3. 
Propiedades de los ﬂuidos. 


1.3.1. 
Densidad ρ 


Se deﬁne como la masa por unidad de volumen. Sus unidades en el sistema internacional son [kg/m3]. 


Para un ﬂuido homogéneo, la densidad no varía de un punto a otro y puede deﬁnirse simplemente mediante 


ρ = V 


m 
(1–1) 


Por el contrario, para un ﬂuido inhomogéneo, la densidad ρ varía de un punto a otro. Por tanto tenemos que deﬁnir 
la densidad en un punto como la masa por unidad de volumen en un elemento diferencial de volumen† en torno a ese 
punto: 


ρ = ρ(x, y, z, t) = dm 


dV 
(1–2) 


Esto es posible gracias a la continuidad. En los líquidos, al tener baja compresibilidad, la densidad depende de la 
temperatura, pero apenas depende de la presión, ρ = ρ(T ). Para los ﬂuidos compresibles, la densidad depende en 
general tanto de la presión como de la temperatura, ρ = ρ(p, T ). Para el caso concreto de un gas ideal, con una 
ecuación de estado pV = nRT , la densidad tiene la forma concreta: 


ρ(p, T ) = Mp 


RT 
(1–3) 


†Nótese que, aunque sea un elemento inﬁnitesimal de volumen, contiene una gran cantidad de partículas y podemos deﬁnir propiedades 
promedio. 
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1.3.2. 
Peso especíﬁco γ 


El peso especíﬁco se deﬁne como el peso por unidad de volumen. En el sistema internacional sus unidades son 
[N/m3]. Para un ﬂuido homogéneo γ = mg/V = ρg, mientras que para un ﬂuido inhomogéneo, 


γ = γ(x, y, z, t) = g dm 


dV = ρg 
(1–4) 


donde g es la aceleración de la gravedad. 


1.3.3. 
Volumen especíﬁco v 


Se denomina volumen especíﬁco al volumen ocupado por la unidad de masa. Para un ﬂuido homogéneo se deﬁne 
como v = V/m = 1/ρ, mientras que en el caso general de un ﬂuido inhomogéneo tendremos que hablar de su valor 
en un punto, 


v = v(x, y, z, t) = dV 


dm = 1 


ρ. 
(1–5) 


En todos los casos, v = 1/ρ. Sus unidades en el sistema internacional son [m3/kg]. 


1.3.4. 
Viscosidad. 


Como se ha dicho en la introducción, la viscosidad reﬂeja la resistencia al movimiento del ﬂuido y tiene un papel 
análogo al del rozamiento en el movimiento de los sólidos. La viscosidad está siempre presente en mayor o menor 
medida tanto en ﬂuidos compresibles como incompresibles, pero no siempre es necesario tenerla en cuenta. En el caso 
de los ﬂuidos perfectos o no viscosos su efecto es muy pequeño y no se tiene en cuenta, mientras que en el caso de los 
ﬂuidos reales o viscosos su efecto es importante y no es posible despreciarlo. En el caso del agua a veces se habla del 
ﬂujo del agua seca para el ﬂujo no viscoso del agua y del ﬂujo del agua mojada para el ﬂujo viscoso. 


1.3.5. 
Presión. 


La presión en un punto se deﬁne como el valor absoluto de la fuerza por unidad de superﬁcie a través de una pequeña 
superﬁcie que pasa por ese punto y en el sistema internacional su unidad es el Pascal (1 Pa=1 N/m2). Mientras que 
en el caso de los sólidos en reposo, las fuerzas sobre una superﬁcie pueden tener cualquier dirección, en el caso de 
los ﬂuidos en reposo la fuerza ejercida sobre una superﬁcie debe ser siempre perpendicular a la superﬁcie, ya que 
si hubiera una componente tangencial, el ﬂuido ﬂuiría. En el caso de un ﬂuido en movimiento, si éste es no viscoso 
tampoco aparecen componentes tangenciales de la fuerza, pero si se trata de un ﬂuido viscoso sí que aparecen fuerzas 
tangenciales de rozamiento. 


De este modo, un ﬂuido en reposo a una presión p ejerce una fuerza −pd⃗S sobre cualquier superﬁcie plana arbitraria 
en contacto con el ﬂuido en el punto, deﬁnida por un vector unitario d⃗S, perpendicular a la superﬁcie. En general, la 
presión en un ﬂuido depende del punto, p = p(x, y, z). Así, para un ﬂuido en reposo la presión se deﬁne como la 
fuerza normal por unidad de superﬁcie. 


1.3.6. 
Compresibilidad. 


Se caracteriza por el coeﬁciente de compresibilidad, κ, deﬁnido como 


κ = − 1 


V 
dV 
dp 
(1–6) 


que representa la disminución relativa del volumen por unidad de aumento de presión. Sus unidades son de inversa de 
presión, en el sistema S.I. [m2/N]. Su inversa, 


K = 1 


κ 
(1–7) 


es el módulo de compresibilidad [N/m2]. Tanto κ como K dependen de la forma en que se realiza el proceso. 
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1.3.7. 
Dilatación térmica. 


Se caracteriza por el coeﬁciente de dilatación de volumen, que representa el aumento relativo del volumen producido 
por un aumento de la temperatura, y está deﬁnida como 


αV = 1 


V 
dV 
dT 
(1–8) 


donde V es el volumen inicial del líquido. Sus unidades son de inversa de grados [K−1] o [◦C−1] y depende de la 
forma en que realiza el proceso. 
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Comenzaremos estudiando el caso del ﬂuido estático, en el cual todas las partículas tienen (en movimiento prome- 
dio) una velocidad nula.† En este ﬂuido en reposo no aparecen fuerzas tangenciales que darían lugar a un movimiento 
tangencial. Así, las fuerzas de contacto ejercidas por un ﬂuido en reposo sobre la pared del recipiente que lo contiene 
son siempre perpendiculares a la misma. Asimismo, si introducimos un cuerpo cualquiera en el interior de un ﬂuido, 
se ejercerá una fuerza sobre toda la superﬁcie del cuerpo en contacto con el ﬂuido y esta fuerza será normal a la 
superﬁcie en cada punto. La fuerza por unidad de superﬁcie dependerá del punto de la superﬁcie, pero para un punto 
dado no dependerá de la orientación de la superﬁcie del cuerpo en ese punto. Denominaremos presión ejercida sobre 
la superﬁcie a la magnitud de esta fuerza normal por unidad de superﬁcie. La presión es una magnitud escalar. 


2.1. 
Ecuación general de la estática de ﬂuidos. Principio de Pascal. 


2.1.1. 
Expresión diferencial de la misma. 


El caso más general 


En el caso más general, la fuerza de presión por unidad de volumen ejercida por el ﬂuido exterior a un pequeño 
elemento de volumen sobre el contenido del mismo, vendrá dada por el gradiente de la presión, ⃗∇p (véase la sección 
3.3.2). Si no hay otras fuerzas como la de la gravedad, la presión debe compensarse para que haya equilibrio. 


†Rigurosamente, se mueven (en movimiento promedio) con una velocidad constante con respecto de un sistema de referencia inercial. 
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dx 


dy 


dz 


Figura 2–1 Obtención de la ecuación general de la estática de ﬂuidos. 


Consideremos una situación en la que estas fuerzas adicionales pueden describirse mediante una energía potencial, 
con Φ el potencial o energía potencial por unidad de masa. La fuerza por unidad de masa vendrá dada en esas circuns- 
tancias por −⃗∇Φ y si ρ es la densidad del ﬂuido, la fuerza por unidad de volumen por −ρ⃗∇Φ. La suma de las fuerzas 
que actúan sobre el elemento de volumen debe ser cero en el equilibrio estático, y por tanto se cumplirá la ecuación 


−⃗∇p − ρ⃗∇Φ = 0 
(2–1) 


que es la ecuación general de la hidrostática. Esta ecuación en el caso general no tiene solución. Si la densidad varía 
en el ﬂuido de forma arbitraria no es posible que las fuerzas estén equilibradas y el ﬂuido no puede estar en equilibrio 
estático, estableciéndose corrientes de convección. Esto es así ya que ⃗∇p es un gradiente puro mientras que el término 
ρ⃗∇Φ no lo es para una densidad ρ variable en general. El equilibrio hidrostático es posible únicamente en el caso de 
ρ constante o cuando ρ es función únicamente de p. Veremos a continuación algunos casos en los que esto es así. 


La ecuación de la hidrostática para ρ = ρ(z) 


Calculemos ahora la dependencia de la presión con la altura en un ﬂuido sometido a la atracción gravitatoria. Para 
ello consideraremos un pequeño elemento prismático inﬁnitesimal de volumen dV cuyas caras son perpendiculares a 
los ejes coordenados y tienen dimensiones dx, dy y dz, con su cara inferior a una altura z. La fuerza que actúa sobre 
la cara inferior del elemento prismático es obviamente 


d⃗F(z) = p(z)dxdy⃗uz 


mientras que la fuerza que actúa sobre la cara superior es 


d⃗F(z + dz) = p(z + dz)dxdy(−⃗uz) = − [p(z) + dp] dxdy⃗uz 


La fuerza de presión neta hacia arriba en la dirección vertical es: 


d⃗Fp = d⃗F (z) + d⃗F(z − dz) = −dpdxdy⃗uz 
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que a su vez debe ser igual al peso del ﬂuido contenido dentro del elemento prismático ρgdxdydz, de donde se obtiene 
la ecuación general de la estática de ﬂuidos en forma diferencial, 


dp = −ρgdz 


donde el signo negativo reﬂeja el hecho de que la presión disminuye al aumentar la altura. 


2.1.2. 
Caso particular: ﬂuido incompresible. 


En el caso de un ﬂuido incompresible la densidad no depende de la presión, y además se tiene una superﬁcie libre. 
Es posible describir el problema en función de la profundidad respecto de la superﬁcie libre h, que aumenta según z 
disminuye. La ecuación general de la estática quedaría entonces, en forma diferencial, como 


dp = −ρgdz = ρgdh 


Esta ecuación se puede integrar con facilidad entre la superﬁcie libre y un punto a una profundidad h, quedando 
h 


0 
dp = p − patm = 
h 


0 
ρgdh = ρgh 


que se escribe como 
p(h) = patm + ρgh = patm + γh 
(2–2) 


y que es la ecuación general de la estática en forma integral para un ﬂuido incompresible. En concreto, entre dos 
puntos del ﬂuido a distinta profundidad, la diferencia de presiones es 


p2 − p1 = ρg(h2 − h1) = −ρg(z2 − z1) 
(2–3) 


Estas expresiones subrayan el importante resultado de que la presión de un ﬂuido en reposo en un punto del mismo 
depende exclusivamente de la profundidad a la que se encuentra dicho punto, y de que la diferencia de presiones entre 
dos puntos depende exclusivamente de la diferencia de profundidades de los mismos (y, obviamente, de la densidad 
del ﬂuido). 


2.1.3. 
Principio de Pascal. 


El resultado anterior nos lleva de forma inmediata al principio de Pascal, 


En un ﬂuido incompresible, las variaciones locales de presión se transmiten íntegramente a todos los 
puntos del ﬂuido y en todos los sentidos, así como a las superﬁcies en contacto con el ﬂuido. 


o dicho con otras palabras, si la presión en un punto aumenta, por ejemplo por la aplicación de una fuerza externa, 
aumenta por igual en los demás puntos del ﬂuido, ya que según (2–3) la diferencia de presiones entre dos puntos 
determinados depende únicamente de ∆h para un ﬂuido estático. 


El hecho de que los líquidos sean compresibles (aunque muy poco) hace que este cambio de presión se transmita 
como una onda hasta que se restablece el equilibrio mecánico, y entonces ya se cumpliría el principio de Pascal. Este 
mismo principio se cumple también para ﬂuidos compresibles, una vez alcanzado el equilibrio. 


2.1.4. 
Paradoja de Pascal. 


A primera vista puede resultar algo sorprendente el hecho de que para un ﬂuido dado, la presión dependa exclusi- 
vamente de la profundidad y no de otras cosas como el tamaño y forma del recipiente, o pueden resultar algo extrañas 
situaciones como la siguiente, conocida como paradoja de Pascal. 


Consideremos los tres recipientes de idéntica base que se presentan en la ﬁgura 2–2. Los tres están llenos de agua 
hasta el mismo nivel, pero sus formas son muy distintas, uno tiene la parte superior muy cerrada, uno es cilíndrico y 
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S 
S 
S 


(a) 
(b) 
(c) 


Figura 2–2 Paradoja de Pascal. A pesar de que la cantidad de ﬂuido en los recipientes sea claramente distinta, la presión en la base (que es la 
misma para los tres) es la misma. 


el otro tiene la parte superior muy abierta, pero todos ellos tienen la misma base.† Puede resultar algo sorprendente el 
hecho de que en los tres casos la fuerza ejercida sobre la base sea la misma, como se desprende de la ecuación general 
de la estática de ﬂuidos. Sin embargo, no es tan sorprendente si se tiene en cuenta que las paredes del recipiente ejercen 
sobre el líquido una fuerza perpendicular a las mismas, que puede tener una componente vertical neta bien hacia abajo 
(primer caso), nula (segundo caso) o bien hacia arriba (tercer caso), componente vertical que es necesario tener en 
cuenta, y que en el primer caso añade un término adicional al peso del ﬂuido, en el segundo caso no afecta y en el 
tercer caso aminora el efecto del peso del ﬂuido. 


2.2. 
Tubos en U y manómetros. 


2.2.1. 
Fluidos miscibles. 


Consideremos el caso de un sistema de vasos comunicantes (o un tubo en U) que contiene un ﬂuido de densidad ρ. 
Si denominamos 0A y 0B a dos puntos en la base que se encuentran a la misma altura y sin ningún obstáculo entre 
medias, está claro que la presión en ambos puntos es la misma, p0A = p0B. Sin embargo queda la duda acerca de 
lo que ocurre con las presiones en los puntos M y N, que se encuentran a una altura d por encima de 0A y 0B y 
que no están comunicados de una forma tan clara y directa como lo están los puntos 0A y 0B. Sin embargo, se ve 
enseguida que si la diferencia de presiones entre M y 0A es pM − p0A = −ρgd y la diferencia de presiones entre M 
y 0B es pM − p0B = −ρgd, al tener en cuenta que p0A = p0B, se obtiene directamente que para un ﬂuido en reposo, 
pN = pM. 


2.2.2. 
Fluidos no miscibles. 


Las cosas no son tan claras sin embargo, si se tienen ﬂuidos no miscibles. Consideremos el caso de la ﬁgura en el 
que se tienen dos ﬂuidos no miscibles, de forma que el ﬂuido en la columna B es más denso que el ﬂuido en la columna 
A. Sea C la línea que pasa por la interfase entre los dos ﬂuidos y une dos puntos a la misma altura. Sean M y N dos 
puntos a la altura del nivel superior de la columna B. Si las dos columnas están abiertas a la atmósfera, las presiones 


†Lo mismo ocurre en el caso de la fuerza sobre una pequeña superﬁcie en la base de una presa llena de agua hasta un cierto nivel y la fuerza 
sobre una superﬁcie similar en un tubo vertical de sección la de la superﬁcie, lleno de agua hasta la misma altura de la presa. 
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replacements 


d 


A 
A 
B 
B 


M 


M 


N 


N 


p 
p 


dA 


pA 
pB 


CA 
CB 


0A 
0A 
0B 
0B 


h 
hA 


hB 
ρ 
ρA 
ρB 


Figura 2–3 Presión en un tubo en U. a) caso de ﬂuidos miscibles. b) Caso de ﬂuidos no miscibles. 


respectivas en M y N serán pN = patm y pM = patm + γAdA, por lo que claramente pM > pN. En líquidos no 
miscibles las presiones a la misma altura son en general distintas. 


El esquema de la ﬁgura nos permite obtener la densidad relativa de los dos ﬂuidos. Como las presiones en la línea 
C son las mismas para la columna A que para la B, al corresponder a un mismo ﬂuido, se tiene que de pCA = 
patm + γAhA y pCB = patm + γBhB se obtiene γAhA = γBhB, lo que permite obtener la densidad relativa en 
función de la altura de las distintas columnas como 


γA 
γB 
= ρA 


ρB 
= hA 


hB 
(2–4) 


2.2.3. 
Manómetros. 


Consideremos ahora el caso de un manómetro que queremos utilizar para medir la presión en un recipiente A en la 
forma que se muestra en la ﬁgura. Como a la altura MN hay un único ﬂuido, pN = pM. Como la presión en M viene 
dada por pM = pA + γAd y la presión en N por pN = patm + γh, se tiene que la presión en el depósito A vendrá 
dada por 
pA = patm + γh − γAd 
(2–5) 


que en el caso habitual de un ﬂuido manométrico mucho más denso que el ﬂuido en A se reduce a 


pA = patm + γAh 
(2–6) 


Esta aproximación es en general válida si γ ≫ γA. En caso contrario es necesario utilizar la primera de las expresio- 
nes. Normalmente con el mercurio como ﬂuido manométrico la aproximación va bastante bien, a no ser que se quiera 
medir la presión de un ﬂuido muy denso, ya que la densidad relativa del mercurio respecto del agua es alta, de 13,6. 


2.2.4. 
Manómetro diferencial. 


Un razonamiento análogo se puede hacer para el caso del manómetro diferencial de la ﬁgura, con una presión 
pM = pA + γAdA + γh en M, y una presión pN = pB + γB en N, lo que da una diferencia de presiones entre A y B 


pB − pA = γh + γAdA − γBdB 
(2–7) 


que para el caso habitual γ ≫ γA, γB queda simplemente como pB − pA = γh 
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M 
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h 
pA 


A 


M 
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h 


pA 


dA 


dB 


B 


pB 


Figura 2–4 Medida de la presión con un manómetro. a) Manómetro normal b) Manómetro diferencial. 


2.3. 
Variación de la presión con la altura en un gas perfecto en reposo. 


Hemos visto que la variación diferencial de la presión con la profundidad viene dada por dp = ρgdh = −ρgdz = 
−γdz y hemos resuelto el problema para un ﬂuido incompresible (γ y ρ no varían con la presión). 


Consideremos ahora el caso de un ﬂuido compresible, y nos restringiremos al caso en que este ﬂuido compresible es 
un gas que se comporta como un gas ideal (ésto ocurre con el aire atmosférico y con buena parte de sus componentes 
para un rango de presión y temperatura relativamente alto). 


Obtengamos en primer lugar la dependencia de la densidad del gas ideal con presión y temperatura. Si tenemos en 
cuenta que ρ = m/V y que n = m/ ¯ 
M, con ¯ 
M el peso molecular promedio del gas ideal, 


pV = nRT = pm 


ρ = m 


¯ 
M 


por lo que la densidad queda: 


ρ(p, T ) = 
¯ 
Mp 
RT . 
(2–8) 


La variación inﬁnitesimal de la presión con la altura se puede entonces escribir como 


dp = −ρgdz = −ρ(p, T )gdz = − 
¯ 
Mp 
RT gdz 
(2–9) 


Agrupando variables, 
dp 
p = − 
¯ 
Mg 
R 
dz 
T 
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que, al integrar queda 


ln p 


p0 
= − 
¯ 
Mg 
R 


z 


0 


dz 
T . 
(2–10) 


Ahora bien, en general T = T (z), es una función de z y por tanto, para poder obtener la variación de la presión con la 
elevación es necesario conocer la función T (z). 


2.3.1. 
Atmósfera isoterma 


Consideremos primero el caso más sencillo, el de una atmósfera isoterma con T = T0 = cte. En este caso, la 
expresión anterior queda 


ln p 


p0 
= − 
¯ 
Mg 
R 
z − z0 


T0 
(2–11) 


que da una variación exponencial de la presión con la altura 


p = p0 exp 
− 
¯ 
Mg 
RT0 
(z − z0) 
(2–12) 


que es una buena aproximación para pequeñas alturas. 


2.3.2. 
Dependencia lineal 


Para grandes alturas, la temperatura de la atmósfera disminuye de forma aproximadamentelineal hasta unos 11000m, 
expresándose esta dependencia en la forma 
T (z) = T0 − βterz 
(2–13) 


con T0 la temperatura absoluta al nivel del mar y βter el denominado gradiente térmico. Estas magnitudes varían de 
un día a otro, aunque como ejemplo, unos valores típicos serían T0 = 288, 16 K∼ 15◦C y βter ∼ 6, 50 Kkm−1. La 
integral queda entonces 


ln p 


p0 
= − 
¯ 
Mg 
R 


z 


0 


d 


T0 − βterz = − 
¯ 
Mg 
R 


− ln(T0 − βterz) 


βter 


z 


0 
= 
¯ 
Mg 
Rβter 
ln T0 − βterz 


T0 
(2–14) 


que se puede escribir como 


ln p 


p0 
= ln 
T0 − βterz 


T0 


¯ 
Mg 
Rβter , 


quedando la dependencia de la presión con la altura, para una dependencia lineal de la temperatura con la altura, en la 
forma 


p = p0 


T0 − βterz 


T0 


¯ 
Mg 
Rβter , 
(2–15) 


2.3.3. 
Atmósfera tipo 


De entre los distintos casos particulares considerados, la atmósfera de densidad homogénea tiene un interés esencial- 
mente teórico, mientras que la aproximación de atmósfera isoterma es aplicable únicamente a capas de la atmósfera 
con densidad constante. Esto último es aplicable a la baja estratosfera, especialmente a altas alturas, pero, salvo que se 
trate de capas delgadas, es difícilmente aplicable a la troposfera. La aproximación de atmósfera adiabática tiene una 
aplicación muy interesante, ya que da un límite superior para el gradiente de una atmósfera verticalmente estable y 
tiene la distribución de temperatura de una capa mezclada verticalmente. 


La atmósfera real tiene, en general, valores inferiores del gradiente térmico, y aunque en general éste no es constante 
con la altura, es posible deﬁnir una atmósfera de gradiente constante (lineal) que se aproxime al caso real. 
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h 
hcg 
hcp 


θ 


ξ 
ξcg 


ξcp 


x 


x 


y 


y 


O 


xcp 


ycp 


CP 
CG 
dA = dxdy 


h = ξ sen θ 
hCG = ξCG sen θ 
hCP = ξCP sen θ 


Figura 2–5 Cálculo de la fuerza sobre una superﬁcie plana. Consideramos una superﬁcie plana que forma un ángulo θ con la superﬁcie libre del 
líquido. 


Se han deﬁnido distintas atmósferas tipo que se utilizan como referencia para una aproximación a la atmósfera real 
y para calibrar y utilizar altímetros. Las condiciones en que se deﬁne la atmósfera tipo son las siguientes, 


• Atmósfera de aire seco puro con composición química constante en sentido vertical y peso molecular promedio 
¯ 
M = 28, 9644g. 


• El comportamiento del gas puede aproximarse por el de un gas ideal. Ésto es razonable para p ≈1 atm. 


• Se utiliza como valor estándar de la aceleración el de la gravedad a nivel del mar, g ≈ 9, 80 ms−2 


• El sistema se encuentra en equilibrio hidrostático. 


• Se toma como referencia la temperatura T0 = 15◦C a nivel del mar y la presión p0 = 1 atm. 


• Para alturas hasta los 11000m (tropopausa), se considera constante el gradiente térmico y de valor βter = 
−dT/dz = 6, 5 K·km−1 


• Para alturas mayores de 11000m (estratosfera) y hasta 20000m se considera la temperatura constante e igual 
a −56, 5◦C. A partir de esa altura se supone un gradiente térmico constante de valor −1, 0 K·km−1 hasta los 
32000 m. 


2.4. 
Fuerzas hidrostáticas sobre superﬁcies. 


2.4.1. 
Fuerza sobre superﬁcies planas. 


Consideremos una placa plana de forma arbitraria que se encuentra sumergida completamente en el líquido, como 
se muestra en la ﬁgura. 
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Valor de la fuerza resultante. 


Tomamos como sistema de referencia el sistema de coordenadas [x, y] centrado en el centro de gravedad de la 
superﬁcie. La fuerza total que actúa sobre la placa se expresa como 


FA = 


A 
pdA = 


A 
(patm + ρgh)dA = patmA + ρg 


A 
hdA 


Si tenemos en cuenta que h = ξ sen θ, y que θ se mantiene constante sobre la placa, la fuerza queda como 


FA = Apatm + ρg 


A 
ξ sen θdA = Apatm + ρg sen θ 


A 
ξdA 


Ahora bien, si recordamos la deﬁnición del centro de gravedad Aξcg = 


A ξdA, la expresión anterior se convierte en 


FA = Apatm + ρg sen θAξcg = Apatm + ρgAhcg = (patm + ρghcg)A 


que es simplemente 
FA = pcgA 
(2–16) 


La fuerza total que actúa sobre una superﬁcie plana cualquiera sumergida en un ﬂuido uniforme es 
igual a la presión que hay en el centro de gravedad de dicha superﬁcie multiplicada por su área, con 
independencia de la forma de la superﬁcie plana o de su inclinación θ. 


Punto de aplicación de la fuerza resultante. 


Imaginemos un recipiente que contiene un líquido y sobre éste una superﬁcie A libre que encaja perfectamente en 
la pared del recipiente. Al poder moverse la superﬁcie, ésta será empujada por el líquido y se saldrá. Para evitarlo 
habría que aplicar sobre la superﬁcie una fuerza normal a la misma de magnitud la fuerza que ejerce el líquido sobre 
la superﬁcie. Para que además la superﬁcie no gire, esta fuerza debe aplicarse en un punto determinado de forma que 
el momento total del sistema de fuerzas ejercido por el líquido sobre la superﬁcie se compense con el momento de la 
fuerza equivalente aplicada en ese punto. Este punto es el centro de presiones. 


Cuanto mayor es la profundidad, mayor es la presión. Por tanto, el punto de actuación de la fuerza total resultante 
(centro de presiones) no coincide con el centro de gravedad, sino que debe encontrarse más abajo. La línea de acción 
de esta fuerza debe pasar precisamente por este centro de presiones. Para calcular la posición (xcp, ycp) del centro de 
presiones se suman los momentos de las fuerzas elementales pdA respecto del centro de gravedad y se igualan con 
el momento (respecto del centro de gravedad) de la fuerza resultante aplicada en el centro de presiones. Obtengamos 
ésto para cada una de las componentes, 


Eje x: 
La componente del momento en la dirección x viene dada por 


F · ycp = 


A 
ypdA 
(2–17) 


Como p = patm + ρgξ sen θ se tiene a su vez que 


F · ycp = 


A 
ypdA = 


A 
y(patm + ρgξ sen θ)dA = patm 



A 
ydA + ρg sen θ 
yξdA 


El primer sumando es cero por la deﬁnición del centro de gravedad, y como y = ξcg − ξ se tiene, con ξ = ξCG − y 


F · ycp = ρg sen θ 
ξcg 



A 
ydA − 


A 
y2dA 


De nuevo el primer sumando es cero por la deﬁnición de centro de gravedad de una superﬁcie, por lo que la componente 
x del momento es 
F · ycp = −ρ sen θIxx 
(2–18) 
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con 
Ixx = 


A 
y2dA 
(2–19) 


el denominado momento de inercia de la sección plana (del área de la placa) respecto de su eje horizontal x, calculado 
en el plano de la placa. Como F = pcgA, la distancia en el plano de la placa a la que se encuentra el centro de presiones 
respecto del centro de gravedad, viene dada por 


ycp = −ρg sen θ 


F 
Ixx = −ρg sen θ 


pcgA Ixx 
(2–20) 


El signo − indica que el centro de presiones está por debajo del centro de gravedad, y como se puede ver, esta posición 
depende de la inclinación θ. 


Al aumentar la profundidad a la que se encuentra la placa, ycp se acerca al centro de gravedad, ya que todos los 
factores que intervienen en 2–20 permanecen constantes excepto pcg que aumenta. 


Eje y: 
La componente del momento en la dirección y viene dada por 


F · xcp = 


A 
xpdA 
(2–21) 


que, al igual que antes se convierte en 


F · xcp = 


A 
xpdA = 


A 
x(patm + ρgξ sen θ)dA = patm 



A 
xdA + ρg sen θ 
xξdA. 


El primer sumando también es cero por la deﬁnición de centro de gravedad y teniendo de nuevo en cuenta que y = 
ξcg − ξ se tiene, al sustituir ξ = ξcg − y 


F · xcp = ρg sen θ 
ξcg 



A 
xdA − 


A 
xydA 
, 


lo que da, para la componente x de la posición del centro de presiones 


F · xcp = −ρ sen θIxy 
(2–22) 


con 
Ixy = 


A 
xydA 
(2–23) 


conocido como el producto de inercia. Dividiendo por la fuerza F = pcg, la componente x de la posición del centro 
de presiones respecto del centro de gravedad en el sistema de coordenadas xy situado sobre la placa es: 


xcp = −ρg sen θ 


F 
Ixy = −ρg sen θ 


pcgA Ixy 
(2–24) 


El producto de inercia puede tener cualquier signo, al contrario que el momento de inercia, que puede ser sólo positivo. 
Dependiendo de este signo el centro de presiones se encontrará a un lado u otro del centro de gravedad. Si Ixy es < 0, 
xcp < 0, mientras que si Ixy es negativo, xcp > 0. En el caso particular en el que Ixy = 0, se tiene que xcp = 0 y, por 
tanto, el centro de presiones está directamente debajo del centro de gravedad sobre el eje y en el plano de la placa. 


Como normalmente la presión atmosférica pa no se tiene en cuenta ya que actúa sobre ambos lados de la placa (cara 
interna de un barco o cara seca de una compuerta o de una presa), se puede escribir pcg = ρghcg y la fuerza neta como 
F = ρghcgA. Las expresiones 2–20 y 2–24 se pueden escribir como 


ycp = −Ixx sen θ 


hcgA 
(2–25) 


xcp = −Ixy sen θ 


hcgA 
, 
(2–26) 


válidas únicamente cuando la superﬁcie superior del líquido está a la presión atmosférica. 
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mg 


ρl´ßqgVdes 


Figura 2–6 Flotación y principio de Arquímedes. 


2.4.2. 
Fuerza sobre superﬁcies curvas. 


Para calcular la fuerza ejercida por el agua sobre una superﬁcie curva hay que tener en cuenta la combinación de dos 
componentes, una horizontal y otra vertical. La componente horizontal se calcula obteniendo la fuerza que actuaría 
sobre la proyección de la superﬁcie curva en el plano vertical, con su valor y su línea de aplicación a través del 
correspondiente centro de presiones. 


La componente vertical se obtiene directamente a partir del peso del agua sobre la superﬁcie curva, o, en su caso del 
empuje vertical del agua sobre la misma, actuando a través del dentro de gravedad del líquido o del líquido desalojado, 
dependiendo del caso. 


2.5. 
Equilibrio de un cuerpo sumergido. 


2.5.1. 
Flotación. Principio de Arquímedes. 


Consideremos un cuerpo parcial o totalmente sumergido en un ﬂuido en reposo (líquido o gas). Este ﬂuido ejerce 
presiones sobre todas las partes de la superﬁcie del cuerpo, mayores cuanto mayor es la profundidad. Se denomina 
empuje sobre el cuerpo sumergido a la fuerza total hacia arriba ejercida por el agua. 


Calculemos el valor de esta fuerza. Si se considera la superﬁcie cerrada que delimita el cuerpo sumergido, la presión 
en cada punto de esa superﬁcie dará lugar a una distribución de fuerzas cuya resultante es precisamente el empuje 
que estamos calculando. Consideremos ahora que en el recipiente hubiera el mismo ﬂuido y hasta el mismo nivel 
que cuando estaba el cuerpo y dentro de este ﬂuido una superﬁcie cerrada imaginaria que coincida con la superﬁcie 
exterior del cuerpo que está sumergido en el caso real. Como la presión depende únicamente de la profundidad, la 
presión en todos los puntos de esta superﬁcie imaginaria es la misma que había en los puntos correspondientes de la 
superﬁcie real, por lo que la fuerza resultante que ejerce el ﬂuido de fuera de la superﬁcie sobre el ﬂuido de dentro 
de la misma será la misma que ejerce el ﬂuido sobre el cuerpo. Ahora bien, ahora es fácil hacer el cálculo, ya que el 
ﬂuido se encuentra en equilibrio mecánico y por tanto la fuerza que compensa el empuje es el peso del propio ﬂuido 
dentro de la superﬁcie imaginaria. 


¿Dónde estará aplicada esa fuerza? Volvamos de nuevo al caso del ﬂuido con la superﬁcie cerrada imaginaria. Si el 
sistema está en equilibrio mecánico no habrá rotación y para ello es necesario que los puntos de aplicación del empuje 
y del peso del ﬂuido desalojado sean el mismo. Por tanto, se tiene el principio de Arquímedes, 


Un cuerpo total o parcialmente sumergido en un ﬂuido es empujado hacia arriba con una fuerza igual 
al peso del ﬂuido desalojado, que actúa verticalmente a través del centro de gravedad del ﬂuido antes de 
ser desplazado. 
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Figura 2–7 Contraejemplo de porqué el centro de presiones está a lo largo de la línea que pasa por el centro de gravedad del ﬂuido desalojado. 


Frecuentemente se tiene la idea de que este empuje actúa a través del centro de gravedad del cuerpo sumergido. Esto 
es en general incorrecto y para verlo más gráﬁcamente consideremos un contraejemplo. Sea un cuerpo con forma de 
pesa de gimnasia con la particularidad de que cada bola está hecha de un material distinto, una de hierro, la otra de 
corcho blanco (poliestireno expandido), ambas del mismo volumen, con una masa y volumen de la barra prácticamente 
despreciables. El centro de gravedad de la pesa está a efectos prácticos en el centro de la bola pesada. ¿Qué ocurriría 
si el centro de aplicación del empuje estuviera también en el centro de gravedad del cuerpo? Ambas fuerzas estarían 
aplicadas sobre el mismo punto por lo que no se produciría rotación de la pesa hasta que la bola pesada quede abajo y 
la ligera arriba, contra lo que muestra la experiencia. Por el contrario, en el caso real en el que el empuje actúa sobre 
el centro de gravedad del ﬂuido desalojado, este par de fuerzas se produce en el sentido adecuado (como debe ser). 


2.5.2. 
Estabilidad de un cuerpo semisumergido 


Al estudiar un cuerpo semisumergido no solamente es importante conocer el valor del empuje, sino cómo actúa 
éste. Así, es de gran importancia conocer la capacidad del cuerpo ﬂotante de recuperar el estado de equilibrio una vez 
sacado del mismo. No nos sirve sólo que el barco ﬂote, sino la posición en que debe hacerlo. 


La estabilidad longitudinal (a lo largo) de los barcos es normalmente muy considerable, por lo que nos centraremos 
en la estabilidad transversal (frente a giros a izquierza o derecha). Consideremos el barco cuya sección transversal 
se esquematiza en la ﬁgura 2–8 y denominaremos eje de ﬂotación OO′ al eje de simetría del barco cuando éste se 
encuentra en la posición de equilibrio. En ésta, tanto el centro de gravedad como el de empuje están sobre este eje. 
Cuando como consecuencia de las fuerzas que actúan sobre el barco, éste se inclina un cierto ángulo α, una parte del 
barco ha salido del agua y otra parte que antes estaba fuera ha entrado. La posición del centro de gravedad del barco 
no ha variado en este giro, pero el centro de gravedad del líquido desplazado sí, quedando más hacia la izquierda. Así, 
el eje de empuje se desplaza, en este caso hacia la derecha. Si este eje de empuje queda a la derecha del centro de 
gravedad se produce un par recuperador que tiende a devolver el barco a su posición original, mientras que si queda 
a la izquierda el par fomentará el vuelco del barco. Denominamos metacentro M al punto de cruce del eje de empuje 
y el eje de ﬂotación del barco, y altura metacéntrica h a la distancia entre el metacentro y el centro de gravedad a lo 
largo del eje de ﬂotación. 


La altura metacéntrica es una magnitud característica de la sección transversal del cuerpo para un peso dado y su 
valor un indicador de la estabilidad del cuerpo. Valores típicos para barcos comerciales están entre 0, 3 y 0, 8 m. El 
cálculo de la altura metacéntrica es, en general, un trabajo pesado. 


Dependiendo de la posición relativa del metacentro y el centro de gravedad se tienen 3 casos: 


Equilibrio estable 
El metacentro se encuentra por encima del centro de gravedad. De esta forma el par de fuerzas es un par 
restaurador que lleva al barco a su posición inicial. 
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Figura 2–8 Estabilidad de un cuerpo semisumergido. 


Equilibrio indiferente 
El metacentro y el centro de gravedad coinciden. No hay par de fuerzas. 


Equilibrio inestable 
El metacentro se encuentra por debajo del centro de gravedad. El par de fuerzas es un par de vuelco. 


Así, la estabilidad del barco será tanto mayor cuanto más bajo se encuentre el centro de gravedad y cuanto mayor sea 
la altura metacéntrica. 


En condiciones de equilibrio estable el barco oscilará con una cierta frecuencia. El problema es similar a un péndulo 
con el metacentro como punto ﬁjo. Así, si MM es el momento de inercia respecto del eje perpendicular al plano de 
oscilación que pasa por el metacentro, el periodo será 


T = 2π 



MM 
mgh 
(2–27) 


donde h es la altura metacéntrica. Al disminuir el producto del peso mg por la altura metacéntrica h el periodo aumenta 
y por tanto el balanceo es algo menos desagradable. Sin embargo en la práctica ésto implica una disminución de la 
estabilidad. 
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3.1. 
Movimiento de un ﬂuido perfecto. 


3.1.1. 
Aproximaciones clásicas al estudio de los ﬂuidos. 


A la hora de abordar el estudio de los ﬂuidos, surgen dos aproximaciones clásicas, la de Euler y la de Lagrange. 


Método de Lagrange. 
El método de Lagrange constituye una generalización directa de la mecánica del punto mate- 
rial. Se estudia un volumen pequeño del ﬂuido y se sigue el movimiento de cada una de las partículas (de coordenadas 
x, y, z) en función del tiempo t, a través de la ecuación que describe la trayectoria de cada una de las partículas. El prin- 
cipal inconveniente de este sistema es que hacen falta una gran cantidad de ecuaciones para describir el movimiento 
del sistema, por lo que en la práctica no es útil. 


Método de Euler. 
En la aproximación de Euler se desiste de describir el movimiento del ﬂuido mediante la historia 
de cada una de las partículas individuales. En su lugar se especiﬁca el movimiento del ﬂuido por la densidad ρ(x, y, z) 
y la velocidad ⃗v(⃗r, t) de las partículas del mismo en ese punto, como una función del tiempo y del espacio. En otras 
palabras, en este método se estudia un punto del espacio y como es el movimiento del ﬂuido en ese punto en función 
del tiempo. Las herramientas de trabajo serán las típicas de la Teoría de Campos, con un campo de presiones, un 
campo de velocidades y un campo de densidades. 
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3.1.2. 
Tipos de ﬂujo. 


Atendiendo a las características del ﬂujo, éste puede clasiﬁcarse de acuerdo con distintos criterios. 


Flujo estacionario/no estacionario. 
Se dice que el ﬂujo es estacionario si la velocidad ⃗v(⃗r) y la densidad ρ(⃗r) del 
ﬂujo en un punto no dependen del tiempo y no estacionario en caso contrario. Esto no quiere decir que la velocidad 
y la densidad deban ser las mismas en dos puntos distintos, sino sólo que en un mismo punto no deben variar con el 
tiempo. 


Flujo irrotacional/rotacional. 
Se dice que el ﬂujo es irrotacional cuando el elemento del ﬂuido en un punto dado 
no tiene una velocidad angular neta alrededor de dicho punto y que es rotacional en caso contrario. Un ﬂuido que 
circula a través de una tubería recta de sección uniforme sería un ejemplo simple de ﬂujo irrotacional, mientras que un 
remolino de un río sería un ejemplo de ﬂujo rotacional. 


Flujo compresible/incompresible. 
Se dice que el ﬂujo es compresible si la densidad ρ en el mismo varía, como 
por ejemplo ocurre en los gases en el caso más general, mientras que se dice que el ﬂujo es incompresible cuando 
la densidad ρ apenas varía como es el caso de los líquidos. Nótese que es posible tener un ﬂujo aproximadamente 
incompresible aunque el ﬂuido en movimiento en sí sea un ﬂuido compresible siempre que a lo largo del ﬂujo en la 
región considerada la densidad ρ sea prácticamente la misma en todos los puntos. 


Flujo viscoso/no viscoso. 
Se dice que el ﬂujo es viscoso cuando aparecen en él importantes fuerzas de rozamiento 
que no se pueden despreciar. Como consecuencia de estas fuerzas de rozamiento aparecen unas fuerzas tangenciales 
entre las capas del ﬂuido en movimiento relativo y hay una disipación de energía mecánica. Por el contrario se dice 
que el ﬂujo es no viscoso cuando estas fuerzas de rozamiento son muy pequeñas o bien no se tienen en cuenta. 


3.1.3. 
Senda, línea de corriente y de traza. Tubo de ﬂujo. 


Vamos a deﬁnir ahora una serie de líneas que nos permitirán describir el movimiento de un ﬂuido. Estas son la 
senda, la traza o línea de traza y la línea de corriente o línea de ﬂujo. En general, cada una de éstas líneas es distinta 
de las otras, pero en el caso de un ﬂujo estacionario las tres coinciden. 


Senda. 


Se denomina al camino seguido realmente por una partícula de ﬂuido. La senda se deﬁne para una partícula a 
lo largo del tiempo. Su carácter es fundamentalmente experimental y se obtendría experimentalmente soltando una 
partícula marcador y haciendo una fotografía a obturador abierto durante el tiempo del estudio. 


Para obtener analíticamente la senda integraríamos el campo de velocidades para obtener las ecuaciones paramétri- 
cas (x(t), y(t), z(t)) de la senda para un elemento genérico del ﬂujo que en un cierto instante de referencia pasa por 
un punto de referencia dado (lo que nos permite obtener las constantes de integración). Si eliminamos el tiempo entre 
las distintas ecuaciones obtenemos la ecuación e la senda. 


Línea de traza. 


Se denomina línea de traza al lugar geométrico de las partículas que en instantes sucesivos pasaron por un punto 
dado. Su carácter es también fundamentalmente experimental y en la práctica se obtendría por ejemplo inyectando de 
forma continua en un punto ﬁjo del ﬂujo una serie de partículas marcadas y tomando una fotografía. 


Para obtener analíticamente la ecuación de la traza debemos previamente obtener la familia de constantes de inte- 
gración correspondientes a las sendas de las partículas que en instantes tt anteriores al instante genérico t considerado 
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1 
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1,5 
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2,5 
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3,5 


3,5 
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Traza 


Traza 


Senda 


Senda(s) 


Líneas de campo 


Líneas de campo 


Figura 3–1 Representación espacial (x, y) de un campo de velocidades y de algunas líneas de ﬂujo junto con varias sendas y una línea de traza. El 
campo de velocidades representado es el campo (x/(1 + t), y, 0) en el instante t = 1 s. La senda representada en trazo más grueso corresponde a 
la seguida por un elemento del ﬂuido que en t = 1 s pasa por el punto (2,7, 2,7) y la traza al lugar geométrico de los elementos del ﬂuido que en 
instantes sucesivos pasan por dicho punto. Junto a la senda anterior se muestran las sendas de las partículas que pasarán por el citado punto en los 
instantes t = 1,5 s, t = 2 s y t = 3 s, representadas con líneas ﬁnas de trazos. 


pasaron por el punto (x0, y0, z0) de referencia. Sustituyendo estos valores en las ecuaciones paramétricas de la senda 
tendríamos las sendas de cada una de estas partículas y eliminando tt tendríamos la ecuación de la traza, en el instante 
t, de las partículas que en instantes previos pasaron por (x0, y0, z0). 


Línea de corriente o línea de ﬂujo. 


Se deﬁne la línea de corriente como una línea que en un instante dado es tangente al vector velocidad en cada punto. 
Es importante recalcar que la línea de corriente está deﬁnida para un instante dado, mientras que la senda y la línea de 
traza contienen información de otros instantes. 


Al contrario de las dos anteriores, de carácter básicamente experimental, la línea de corriente tiene un profundo 
sustrato matemático. Para su obtención experimental se extenderían en el ﬂujo un conjunto de partículas marcadoras 
y se tomaría una fotografía del conjunto “casi” instantánea (obturador abierto durante un breve periodo de tiempo), 
de modo que en ésta se aprecian para cada partícula líneas pequeñas que marcan la dirección (ángulo de la línea) y el 
módulo de la velocidad (relacionado con la longitud de la línea). 


Como la velocidad debe ser tangente en cada punto a la línea de ﬂujo, debe cumplirse ⃗v × d⃗r = 0 con d⃗r = 
(dx, dy, dz). Así, éstas deben cumplir 
dx 
vx 
= dy 


vy 
= dz 


vz 
= dr 


v 
(3–1) 
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A1 
A2 


vdt 


Figura 3–2 Tubo de corriente utilizado para la obtención de la ecuación de continuidad en forma integral. 


que, si vx, vy vz son conocidas en función de la posición y el tiempo, pueden integrarse para dar la ecuación de la 
línea de corriente que en un determinado instante pasa por un punto dado. Nótese que esta integración puede ser muy 
compleja. 


Tubo de ﬂujo. 


Un tubo de ﬂujo o de corriente está formado por un haz de líneas de ﬂujo. Como las líneas de ﬂujo son tangentes 
al vector velocidad en cada punto, las líneas de ﬂujo no atraviesan las paredes del tubo de ﬂujo. Esto es así porque 
las paredes están a su vez formadas por líneas de ﬂujo y si las atravesaran dos líneas de ﬂujo se cortarían, estando la 
velocidad indeterminada en el punto de corte. Por tanto, a pesar de estar limitado por la superﬁcie ﬁcticia que envuelve 
el haz de líneas de ﬂujo, éste se comporta a todos los efectos como una superﬁcie impermeable. 


3.2. 
Ecuación de continuidad. 


3.2.1. 
Forma integral. 


Consideremos un tubo de corriente estrecho, de forma que se pueda considerar uniforme la velocidad en cualquier 
sección del tubo perpendicular al ﬂujo. En el interior del tubo la velocidad del ﬂujo es paralela a la línea de corriente 
en cada punto, pudiendo ser estas velocidades distintas en cada punto. 


Sea v1 la velocidad de la partícula en el punto 1, y v2 la velocidad de la partícula en el punto 2, con A1 y A2 
las secciones transversales de los tubos, perpendiculares a las líneas de corriente. Si el tubo es estrecho v1 y v2 son 
uniformes en A1 y A2 respectivamente. 


En un intervalo de tiempo dt, un elemento del ﬂuido recorrerá una distancia v dt, por lo que en el tiempo dt pasará 
por A1 la masa de ﬂuido 
dm1 = ρ1 · A1 · v1 dt 


donde ρ1 es la densidad del ﬂuido al pasar por la sección 1. El ﬂujo de masa o caudal másico se deﬁne como la masa 
que atraviesa una sección en la unidad de tiempo, y viene dado por 


Qm = dm1 


dt 
= ρ1A1v1, 
(3–2) 


donde se considera implícitamente que en ese intervalo inﬁnitesimal de tiempo ni A ni v varían apreciablemente en el 
recorrido del ﬂuido v dt. El caudal másico a través de la sección A1 es ρ1A1v1 y a través de la sección A2 es ρ2A2v2. 
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S2 


Figura 3–3 


Como las partículas del ﬂujo no pueden atravesar las paredes del tubo de ﬂujo debe cumplirse que, si el régimen es 
permanente (o estacionario) y no hay fuentes ni sumideros de partículas, ambos caudales másicos han de ser iguales 


Qm = ρ1A1v1 = ρ2A2v2 
⇒ 
Qm = ρAv = cte 


y análogamente para cualquier otra sección A perpendicular al tubo de ﬂujo, por lo que esta ley de conservación de la 
masa o ecuación de continuidad se puede escribir simplemente como 


ρAv = cte 
(3–3) 


a través de cualquier sección del tubo de ﬂujo perpendicular al mismo en régimen estacionario. 


Para el caso particular de ﬂujo incompresible ρ no depende del punto y esta ecuación de continuidad puede escribirse 
como 
A1v1 = A2v2 
⇒ 
Q = Av = cte 
(3–4) 


donde Q es el caudal o volumen que atraviesa la sección en la unidad de tiempo. 


Por ejemplo, en una canalización por la que circula un ﬂuido incompresible, se tiene la sencilla relación S1v1 = 
S2v2, que da, para la relación entre velocidades 


v2 = S1 


S2 
v1 
(3–5) 


En todo el cálculo anterior hemos considerado implícitamente que la velocidad v es uniforme en cada sección. Esto 
no es cierto en el caso general, pero la ecuación de continuidad sigue siendo válida en las mismas condiciones si la 
densidad es uniforme en la sección y en vez de la velocidad se utiliza la velocidad promedio en la sección, 


¯v = 1 


A 



A 
v dA 
(3–6) 


3.2.2. 
Forma diferencial. 


Estudiaremos la conservación de la masa en un prisma inﬁnitesimal de aristas dx, dy y dz ﬁjo en el espacio. El 
caudal másico que atraviesa la cara perpendicular al eje y de la izquierda es 


dQmy = ρvydAy = ρvydxdz 


mientras que el caudal másico que atraviesa la cara perpendicular al eje y de la derecha (a y + dy) es 


dQmy+dy = ρvydAy + ∂ 


∂y[ρvydAy]dy = [ρvy + ∂ 


∂y(ρvy)dy]dxdz 


De este modo, el caudal másico neto a través de las caras perpendiculares al eje y es 


Qneto 
my 
= Qmy − Qmy+dy = − ∂ 


∂y(ρvy)dxdydz 
(3–7a) 
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x 
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y 


y 


z 
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dx 
dy 


dz 


Figura 3–4 Obtención de la ecuación de continuidad en forma diferencial. 


Una expresión similar se obtiene para los caudales másicos netos a través de las caras perpendiculares al eje x, 


Qneto 
mx = Qmx − Qmx+dx = − ∂ 


∂x(ρvx)dxdydz, 
(3–7b) 


y otra para los caudales másicos netos a través de las caras perpendiculares al eje z, 


Qneto 
mz 
= Qmz − Qmz+dz = − ∂ 


∂z(ρvz)dxdydz. 
(3–7c) 


Si tenemos en cuenta que la masa total contenida en el prisma es dm = ρdV = ρdxdydz, el caudal másico total 
también debe ser 


dQneto 
total = ∂ 


∂t(ρdV ) 


que representa la velocidad con que varía la masa en el punto ﬁjo del espacio en el que estamos considerando el 
elemento de volumen.† Este caudal másico debe ser a su vez la suma de los caudales másicos a través de las caras 
perpendiculares a cada una de las direcciones x, y y z. 


dQneto 
total = ∂ 


∂t(ρdV ) = ∂ρ 


∂t dxdydz = 
− ∂ 


∂x(ρvx) − ∂ 


∂y (ρvy) − ∂ 


∂z (ρvz) 
dxdydz. 
(3–8) 


La parte entre corchetes no es más que menos la divergencia de ρ⃗v por lo que la expresión anterior se puede escribir 
como 
∂ρ 
∂t = −∇(ρ⃗v) 
(3–9) 


o 
∂ρ 
∂t + ⃗∇(ρ⃗v) = 0. 
(3–10) 


†Nótese que estamos estudiando un elemento del espacio ﬁjo en el mismo y por ello aparecen derivadas parciales. Véase la discusión sobre la 
velocidad en torno a la ecuación (3–23). 
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Figura 3–5 Variación de la cantidad de movimiento. 


que es la expresión de la ecuación de continuidad en forma diferencial. Para el caso particular de un ﬂuido incompre- 
sible, la densidad ρ no depende de la posición ni del tiempo y la ecuación de continuidad queda como 


⃗∇⃗v = 0. 
(3–11) 


En régimen estacionario, la densidad en un punto no cambia con el tiempo, por lo que ∂ρ/∂t = 0 y la ecuación de 
continuidad queda ⃗∇(ρ⃗v) = 0. 


3.3. 
Fuerza y aceleración en un elemento de ﬂuido. 


3.3.1. 
Fuerza neta sobre un elemento de un tubo de ﬂujo. Variación de la cantidad de movimiento. 


Calculemos la fuerza total que actúa sobre un elemento del ﬂuido que circula por un tubo de ﬂujo estrecho en 
régimen estacionario. Para ello, calcularemos la variación de la cantidad de movimiento de dicho elemento por unidad 
de tiempo, lo que nos dará la fuerza neta que actúa sobre el mismo. 


Consideremos ahora el elemento de ﬂuido que se muestra en la ﬁgura 3–5 y que inicialmente se encuentra entre las 
secciones 1 y 2. Un intervalo de tiempo dt después, el ﬂuido habrá avanzado hasta las secciones 1′ y 2′. La cantidad 
de movimiento inicial del elemento será la suma de las cantidades de movimiento de los subelementos 1 − 1′ y 1′ − 2, 


(m⃗v )inicial = dm1⃗v1 + (m⃗v )1′−2 


mientras que la cantidad de movimiento del elemento pasado un intervalo de tiempo dt será la suma de las cantidades 
de movimiento de los subelementos 1′ − 2 y 2 − 2′, 


(m⃗v )final = (m⃗v )1′−2 + dm2⃗v2, 


siendo la cantidad de movimiento de la parte 1′−2 la misma en ambos casos al ser el régimen permanente. La variación 
de la cantidad de movimiento en dicho intervalo es entonces 


d(m⃗v ) = (m⃗v )final − (m⃗v )inicial = dm2⃗v2 − dm1⃗v1. 


La ecuación de continuidad (3–3) nos dice que las masas dm1 y dm2 deben ser iguales, por lo que la variación de 
cantidad de movimiento por unidad de tiempo puede escribirse como 


d 
dt(m⃗v ) = ⃗Ftot = dm 


dt (⃗v2 − ⃗v1) = Qm(⃗v2 − ⃗v1) 
(3–12) 


donde Qm es el caudal másico (3–2), y esta variación de cantidad de movimiento por unidad de tiempo debe ser igual 
a la fuerza neta que actúa sobre el elemento. 
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Figura 3–6 Volumen de control utilizado para la obtención de la ecuación de Euler. 


3.3.2. 
Ecuaciones de movimiento microscópicas. Aceleración. Ecuación de Euler. 


Seguimos considerando el caso de un ﬂuido perfecto (no viscoso). Para resolver el problema del movimiento de 
un ﬂuido perfecto hay que escribir las ecuaciones de movimiento que relacionan las velocidades con las fuerzas que 
actúan en el ﬂuido. 


Recordemos que, para un ﬂuido perfecto, las fuerzas en el seno del ﬂuido en movimiento son siempre perpendicu- 
lares a la superﬁcie, son las fuerzas de presión. El ﬂuido que se encuentra del lado de la superﬁcie a donde apunta ⃗n 
ejerce una fuerza ⃗F = −p dS ⃗n sobre la superﬁcie dS. 


Consideremos el elemento de volumen diferencial de la ﬁgura,† de dimensiones dx dy dz y caras perpendiculares a 
los ejes coordenados. En la dirección y actúa, sobre la cara que se encuentra a una distancia y del origen, una fuerza 


dF(y) = p(y) dx dz 
(3–13) 


mientras que sobre la cara que se encuentra a una coordenada y + dy actúa una fuerza 


dF(y + dy) = −p(y + dy) dx dz = − 
p(y) + ∂p 


∂y dy 
dx dz. 
(3–14) 


La componente y de la fuerza de presión neta será la suma de las dos fuerzas anteriores, 


dFy = −p(y) dx dz − ∂p 


∂y dy dx dz + p(y) dx dz = −∂p 


∂y dy dx dz 
(3–15a) 


obteniéndose de la misma manera las expresiones para las componentes netas de la fuerza en las direcciones x y z, 


dFx = −∂p 


∂y dy dx dz 
(3–15b) 


dFz = −∂p 


∂y dy dx dz 
(3–15c) 


†Ahora, este elemento de volumen es un elemento de volumen del ﬂuido que se desplaza con el tiempo y por tanto debemos utilizar una 
derivada total. 
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La suma vectorial de estas componentes da la fuerza de presión total que actúa sobre el elemento de volumen, de masa 
ρ dV = ρ dx dy dz. 


d⃗Fpres. = − 
∂p 


∂x 
⃗i + ∂p 


∂y 
⃗j + ∂p 


∂z 
⃗k 
dx dy dz = −⃗∇p dx dy dz. 
(3–16) 


Si recordamos la segunda de las leyes del movimiento de Newton, 


⃗F = m⃗a = md⃗v 


dt 


la ecuación de movimiento para los ﬂuidos perfectos queda 


d⃗Ftot = d⃗Fpres. + ⃗fext dm = d⃗v 


dt dm = −⃗∇p dx dy dz + ⃗fextρ dx dy dz 


que se puede escribir como 
d⃗v 
dt = −1 


ρ 
⃗∇p + ⃗fext 
(3–17) 


donde ⃗fext representa la aceleración debida a una fuerza externa que actúa sobre el elemento de ﬂuido. Normalmente, 
⃗fext = ⃗g = −g⃗uz. 


La ecuación anterior no es más que la ecuación de movimiento ⃗F = m⃗a en la que, junto a la contribución de la 
fuerza de masa a través de ⃗fext aparece un término debido a las diferencias de presión entre dos puntos. Como ⃗∇p está 
dirigido hacia el sentido de presiones crecientes y el movimiento tiene lugar empujado por las altas presiones hacia 
las bajas presiones, aparece el signo negativo delante del gradiente. 


3.3.3. 
Aplicación a un ﬂuido en reposo 


Si además de que la única aceleración de masa a que está sometido el ﬂuido es la aceleración gravitatoria, éste está 


en reposo, ⃗v es cero en todos los puntos del ﬂuido así que d⃗v 


dt = 0 y la ecuación anterior queda 


−⃗∇p + ρ⃗g = 0 
(3–18) 


de componentes 


− ∂p 


∂x = 0 


−∂p 


∂y = 0 


−∂p 


∂z − ρg = 0. 


De las dos primeras expresiones de obtiene directamente que la presión en este caso no puede depender ni de la 
variable x ni de la variable y, es decir, debe ser uniforme en un plano horizontal. Como en este caso, p no depende 
ni de x ni de y, dp/dz = ∂p/∂z y de la última expresión se obtiene la conocida ecuación general de la estática de 
ﬂuidos, dp = −ρgdz. Esta ecuación queda, para un ﬂuido incompresible (cuya densidad ρ no varía), en la conocida 
forma p = p0 − ρgz. 


3.3.4. 
Aplicación a un ﬂuido en reposo en un sistema de referencia no inercial 


Superﬁcie libre en un depósito que se mueve con movimiento rectilíneo acelerado. 


Estudiemos ahora el caso que se presenta en la ﬁgura 3–7, donde un depósito que contiene un líquido de densidad ρ 
se mueve con movimiento acelerado. 
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x 


z 


⃗a 


⃗a 
⃗g 


(a) 
(b) 


Figura 3–7 Se tiene un depósito que se mueve con una aceleración ⃗a (a) y se quiere saber la distribución de presiones en el líquido y la ecuación 
de la superﬁcie libre. Visto el problema desde el sistema de referencia no inercial que se mueve con el depósito (b), sobre éste y su contenido actúa 
una fuerza ﬁcticia igual a −⃗a además de la fuerza gravitatoria. 


En vez de estudiar el problema desde el sistema de referencia inercial lo haremos desde un sistema de referencia que 
se mueve con el depósito. Visto desde este sistema de referencia, aparece una fuerza ﬁcticia con aceleración −⃗a en el 
sentido contrario a ⃗a, y el sistema se encuentra en reposo. Las ecuaciones de Euler aplicadas a este problema nos dan 


d⃗v 
dt = 0 = −1 


ρ 
⃗∇p + ⃗fext. 


En el sistema de referencia no inercial, ⃗fext = −a⃗ux − g⃗uz e igualando las distintas componentes se tiene 


−1 


ρ 
∂p 
∂x = a 
⇒ 
∂p 
∂x = −ρa 


−1 


ρ 
∂p 
∂y = 0 
⇒ 
∂p 
∂y = 0 


−1 


ρ 
∂p 
∂z = g 
⇒ 
∂p 
∂z = −ρg 


La variación inﬁnitesimal de presión entre dos puntos muy próximos vendrá dada por 


dp = ∂p 


∂xdx + ∂p 


∂ydy + ∂p 


∂z dz = −ρadx − ρgdz, 


que, tomando el origen de coordenadas como se indica en la ﬁgura da, al integrar, 


p(x, z) = −ρax − ρgz + p0 
(3–19) 


donde p0 es la presión en el origen de coordenadas, en este caso la presión atmosférica. Para un valor de x dado, según 
aumenta la profundidad (z cada vez más negativo), la presión aumenta linealmente con la profundidad, mientras que 
para un z dado, como la superﬁcie libre no es horizontal en este sistema de referencia, la presión disminuye según 
aumenta x. 


La superﬁcie libre será el lugar geométrico de los puntos que se encuentran a la presión atmosférica, es decir, 
aquéllos para los que al ser p = p0 se cumple 


p − p0 = 0 = −ρax − ρgz 


lo que da como ecuación de la superﬁcie libre 


z(x) = −a 


g x, 
(3–20) 


que es un plano inclinado un ángulo arctana/g. 


Agustin Martin Domingo 


3.3. Fuerza y aceleración en un elemento de ﬂuido. 
35 


r 


r 


z 


(a) 
(b) 


⃗g 


ω 


ωr2 


Figura 3–8 Superﬁcie libre en un depósito circular que gira sobre su eje. 


Superﬁcie libre en un depósito circular que gira sobre su eje. 


Estudiaremos ahora otro caso de sistema acelerado, pero en este caso con aceleración angular. Consideremos un 
depósito cilíndrico colocado en posición vertical que gira alrededor de su eje con una velocidad angular ω. 


En el sistema de referencia del cilindro que rota aparece una fuerza ﬁcticia de masa, la fuerza centrífuga, como 
consecuencia de que este sistema de referencia no es un sistema de referencia inercial. Así, la aceleración total debida 
a las fuerzas de masa en el sistema de referencia no inercial es ⃗fext = ω2r⃗ur − g⃗uz y las distintas componentes de las 
ecuaciones de Euler (en coordenadas cilíndricas) son: 


−1 


ρ 
∂p 
∂r + ω2r = 0 ⇒ ∂p 


∂r = ρω2r 


−1 


ρ 
1 
r 
∂p 
∂θ = 0 ⇒ ∂p 


∂θ = 0 


−1 


ρ 
∂p 
∂z − g = 0 ⇒ ∂p 


∂z = −ρg 


Una variación inﬁnitesimal de p en el entorno de un punto se representaría por 


dp = ∂p 


∂r dr + ∂p 


∂θdθ + ∂p 


∂z dz = ρω2rdr − ρgdz 


que, integrando en el sistema de ejes de la ﬁgura, da 


p(r, z) = p0 + 1 


2ρω2r2 − ρgz 
(3–21) 


donde p0 es la presión en el origen de coordenadas, en este caso la atmosférica. Para un z dado la presión aumenta 
cuadráticamente con r, ya que cada vez hay más ﬂuido encima del punto, mientras que para un r dado, la presión 
aumenta linealmente con la profundidad. 


La superﬁcie libre será el lugar geométrico de los puntos que se encuentren a la presión p0, en este caso la atmosfé- 
rica. Para estos puntos se tiene 


p − p0 = 0 = 1 


2ρω2r2 − ρgz 


lo que da, para la ecuación de la superﬁcie libre, 


z(r) = ω2r2 


2g 
(3–22) 


que es un paraboloide (y en el plano de la ﬁgura se ve como una parábola). 
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3.3.5. 
Vorticidad 


Volvamos sobre la ecuación de Euler (3–17) y discutamos sobre la aceleración d⃗v/dt. En principio podría pensarse 
que la expresión para la aceleración es de la sencilla forma ∂⃗v/∂t, pero esto no es así. ∂⃗v/∂t representa la rapidez con 
que la velocidad cambia en un punto ﬁjo del espacio, pero lo que necesitamos es la variación de la velocidad de un 
elemento particular del ﬂuido. Si la trayectoria tiene una curvatura o el ﬂuido se acelera de un punto a otro, aparece 
una componente de aceleración del elemento aunque la velocidad en un punto se mantenga constante. Así, la variación 
de velocidad de un punto a otro inﬁnitamente próximo viene dada por 


d⃗v = ∂⃗v 


∂t dt + ∂⃗v 


∂xdx + ∂⃗v 


∂y dy + ∂⃗v 


∂z dz 


de forma que la aceleración d⃗v/dt se escribe como 


d⃗v 
dt = ∂⃗v 


∂t + ∂⃗v 


∂xvx + ∂⃗v 


∂yvy + ∂⃗v 


∂z vz, 


que puede escribirse, representando ⃗∇ como un vector, en la forma 


d⃗v 
dt = (⃗v · ⃗∇)⃗v + ∂⃗v 


∂t . 
(3–23) 


Debe observarse que es posible una aceleración aunque ∂⃗v/∂t sea 0, es decir, aunque la velocidad del ﬂuido en un 
punto dado no esté variando. Por ejemplo, un elemento de agua que ﬂuye en círculo a una velocidad constante está 
acelerándose sometido a una aceleración centrípeta aunque la velocidad en un punto dado no cambie. 


Se denomina aceleración convectiva o aceleración de transporte al primero de los sumandos que aparecen en la 
ecuación (3–23) y corresponde al ritmo de cambio de la velocidad debido al cambio en la posición en que se mide. El 
segundo sumando es la denominada aceleración local y corresponde al ritmo de cambio de la velocidad en un punto 
concreto del ﬂujo. 


Utilizaremos la siguiente identidad del análisis vectorial, 


(⃗v · ⃗∇)⃗v = (⃗∇ × ⃗v) × ⃗v + 1 


2 
⃗∇(⃗v · ⃗v) 


y deﬁniremos un nuevo campo vectorial ⃗Ω como el rotacional de ⃗v, 


⃗Ω = ⃗∇ × ⃗v, 
(3–24) 


lo que nos permite escribir la identidad vectorial anterior como 


(⃗v · ⃗∇)⃗v = ⃗Ω × ⃗v + 1 


2 
⃗∇v2. 


Si la fuerza m⃗fext admite un potencial, la aceleración ⃗fext también admitirá uno. Así, la aceleración ⃗fext puede 
escribirse en función de su potencial Φ como ⃗fext = −⃗∇Φ. Si la única fuerza de masa que actúa es la gravitatoria, el 
potencial es simplemente Φ = gz. En función del potencial, la ecuación de movimiento (3–17) queda en la forma 


d⃗v 
dt = ∂⃗v 


∂t + ⃗Ω × ⃗v + 1 


2 
⃗∇v2 = −1 


ρ 
⃗∇p − ⃗∇Φ 
(3–25) 


Se denomina vorticidad al campo vectorial ⃗Ω. Si la vorticidad es cero en cualquier punto (es decir, si ⃗∇ × ⃗v = 0) se 
dice que el ﬂujo es irrotacional. 


Retomemos ahora un concepto habitual en electromagnetismo, la circulación de un campo vectorial. Ésta se deﬁne 
como la integral de línea, a lo largo de cualquier camino cerrado del campo vectorial correspondiente multiplicado 
por el elemento de longitud a lo largo de ese camino. En particular, la circulación del vector velocidad ⃗v a lo largo de 
cualquier camino cerrado en un ﬂuido se deﬁne como 
⃗v d⃗l. 
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Utilizando el teorema de Stokes esta circulación por unidad de área es igual a ⃗∇ × ⃗v. Así, la vorticidad ⃗Ω es la 
circulación del vector velocidad alrededor de un área unidad (perpendicular a la dirección de ⃗Ω). También puede 
deducirse que un elemento de ﬂuido en un punto del líquido rotará con una velocidad angular Ω/2. 


Particularicemos ahora al caso de un ﬂuido incompresible en el que nuestro interés se centra en el campo de velo- 
cidades ya que en estas condiciones es posible eliminar la presión de las ecuaciones anteriores. Para ello se toma el 
rotacional de los dos miembros de la ecuación (3–25). Como el rotacional de cualquier gradiente es cero y en un 
líquido la densidad es constante, se tiene, utilizando la ecuación de continuidad 


⃗∇⃗v = 0 
(3–11) 


para un ﬂuido incompresible, 
∂⃗Ω 
∂t + ⃗∇ × (⃗Ω × ⃗v) = 0, 
(3–26) 


ecuación que, junto con las ecuaciones 
⃗∇⃗v = 0 
(3–11) 


y 
⃗Ω = ⃗∇ × ⃗v 
(3–24) 


determina completamente el campo de velocidades ⃗v. Así, si conocemos la vorticidad ⃗Ω en un determinado instante 
se tiene el rotacional del vector velocidad y como además conocemos que la divergencia es cero, dada esa situación 
física es posible determinar la velocidad ⃗v en cualquier instante y lugar. Conociendo ⃗v la ecuación (3–26) permite 
obtener la vorticidad ⃗Ω en un instante de tiempo posterior. Utilizando de nuevo (3–24) se tiene la nueva velocidad y 
vuelta a empezar. Hay que notar que siguiendo este procedimiento se obtiene información únicamente sobre el campo 
de velocidades y se ha perdido por completo la información sobre la presión. 


Veamos un caso particular de lo anterior. Si ⃗Ω es cero en todo el ﬂuido en un determinado instante de tiempo, ∂⃗Ω/∂t 
también lo será, según la ecuación (3–26), y por tanto ⃗Ω será cero en cualquier instante de tiempo inmediatamente 
posterior. Así, el ﬂujo es permanentemente irrotacional. Si la rotación del ﬂujo es inicialmente cero, siempre será 
irrotacional y las ecuaciones a resolver serán 


⃗∇⃗v = 0 
⃗∇ × ⃗v = 0 


similares a las ecuaciones de los campos electrostático y magnetostático. 


3.4. 
Conservación de la energía. Ecuación de Bernoulli. 


Hemos obtenido la ecuación de continuidad para ﬂujo estacionario y hemos demostrado que, para un ﬂujo incompre- 
sible vA es una constante. Si además, el trabajo realizado por las fuerzas no conservativas es despreciable, la energía 
mecánica se conserva. 


Estudiaremos esto de dos formas, en primer lugar utilizando la ecuación de movimiento del ﬂuido perfecto escrita 
en función de la vorticidad, (3–25) y en segundo lugar mediante el análisis de las variaciones de energía en el ﬂuido 
contenido dentro de un tubo de ﬂujo estrecho y del trabajo mecánico realizado sobre el mismo, 


Obtención a partir de las ecuaciones de Euler 


Consideremos la ecuación (3–25) para el caso particular de ﬂujo estacionario en el que ∂v/∂t = 0 y recordemos 
que en ﬂujo estacionario, las trayectorias reales de las partículas (la senda) tienen lugar a lo largo de las líneas de ﬂujo. 
Multiplicando escalarmente la ecuación (3–25) por ⃗v el término ⃗v(⃗Ω × ⃗v) se anula y queda 


⃗v · ⃗∇ 
p 


ρ + 1 


2v2 + Φ 
= 0 
(3–27) 


Agustin Martin Domingo 


38 
Capítulo 3. Dinámica de ﬂuidos perfectos. 


Si consideramos un sistema de ejes coordenados que tiene uno de ellos (denominémosle xv) dirigido a lo largo de una 
línea de ﬂujo y los otros dos en direcciones perpendiculares, la única componente del gradiente que no se anula en el 
producto escalar es precisamente la componente en la dirección del ﬂujo, quedando la ecuación anterior en la forma 


v ∂ 


∂xv 


p 


ρ + 1 


2v2 + Φ 
= 0. 


o 
∂ 


∂xv 


p 


ρ + 1 


2v2 + Φ 
= 0 
(3–28) 


Esta ecuación nos dice que a lo largo de una línea de corriente el término entre corchetes debe ser constante, 


p 
ρ + 1 


2v2 + Φ = cte 
en una línea de ﬂujo 
(3–29) 


conociéndose a esta expresión como ecuación de Bernoulli. Esta ecuación es válida a lo largo de una línea de ﬂujo 
para un ﬂujo estacionario e incompresible que conserva la energía. La constante será en general distinta para distintas 
líneas de ﬂujo. 


Un caso particular interesante se tiene cuando el ﬂujo es además irrotacional. En este caso la vorticidad es nula en 
cualquier punto del ﬂujo y la ecuación (3–25) queda como 


⃗∇ 
p 


ρ + 1 


2v2 + Φ 
= 0. 


de donde se deduce que en este caso particular de ﬂujo irrotacional, 


p 
ρ + 1 


2v2 + Φ = cte 
en el ﬂujo 
(3–30) 


para puntos en cualquiera de las líneas de ﬂujo del ﬂuido. Si como es habitual la única fuerza de masa que actúa sobre 
el ﬂuido es la atracción gravitatoria, Φ = gz y la ecuación de Bernoulli queda 


p 
ρ + 1 


2v2 + gz = cte 
(3–31) 


o, escrito en la forma habitual, 


p + ρgz + 1 


2ρv2 = cte, 
(3–32) 


o también como 
p 
ρg + z + v2 


2g = cte, 
(3–33) 


En esta última forma cada uno de los elementos de la ecuación de Bernouilli tiene dimensiones de longitud. Física- 
mente eso representa energías por unidad de peso del ﬂuido. Además, cada uno de los sumandos de (3–33) recibe un 
nombre especial, 


p 
ρg es la altura de presión. Representa la energía por unidad de peso debida a las fuerzas de presión. 


z es la altura geométrica y representa la energía potencial gravitatoria del ﬂuido por unidad de peso. 


v2 
2g es la altura de velocidad o altura cinética y representa la energía cinética por unidad de peso del ﬂuido en un punto. 


La suma de las alturas de presión y geométrica p 


ρg + z recibe el nombre de altura piezométrica y corresponde a la 


altura que alcanzaría el líquido en un tubo vertical abierto a la atmósfera (piezómetro) colocado en ese lugar del ﬂujo. 
Asimismo se denomina altura total a la suma de los tres elementos de la ecuación de Bernoulli, correspondiendo a la 
energía total por unidad de peso del ﬂuido. 
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v1dt 
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p1A1 
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z2 − z1 


Figura 3–9 Obtención de la ecuación de Bernoulli para una línea de ﬂujo. En la ﬁgura se representa un tubo de ﬂujo estrecho de forma que la 
velocidad pueda considerarse uniforme en cada sección. 


Obtención a través de la conservación de la energía 


Veamos cómo la ley de conservación de la energía relaciona la altura del ﬂuido, el módulo de la velocidad del ﬂujo 
y la presión para puntos que se encuentran a lo largo de una línea de ﬂujo. Para ello consideraremos un tubo de ﬂujo 
estrecho alrededor de la línea de ﬂujo considerada, como se muestra en la ﬁgura 3–9. De este modo la velocidad en 
las secciones transversales del tubo de ﬂujo puede ser considerada uniforme. 


En primer lugar hallaremos el trabajo realizado en un intervalo de tiempo dt sobre el ﬂuido que está en la región 
limitada por las secciones A1 y A2 y el tubo de corriente. El ﬂuido que se encuentra a la izquierda de la superﬁcie A1 
ejerce sobre ésta una fuerza p1A1, perpendicular a la superﬁcie. En el intervalo de tiempo dt esta fuerza realizará un 
trabajo 
F1dx1 = p1A1v1dt 


sobre el volumen de ﬂuido considerado, aumentando en esa cantidad la energía del mismo. Análogamente, el ﬂuido 
que se encuentra a la derecha de la superﬁcie A2 ejerce sobre ésta una fuerza p2A2, por lo que el trabajo realizado por 
esta fuerza sobre el elemento de volumen considerado es 


F2dx2 = −p2A2v2dt, 


por lo que el trabajo total realizado sobre el ﬂuido es 


δW = p1A1v1dt − p2A2v2dt, 
(3–34) 


expresión válida para cualquier tipo de ﬂuido. Para el caso concreto de un ﬂuido incompresible, se cumple v1A1 = 
v2A2, por lo que la expresión anterior puede escribirse como 


δW = (p1 − p2)v1A1dt = (p1 − p2)v2A2dt 


o simplemente como 
δW = (p1 − p2) dV 
(3–35) 
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donde dV es el volumen del elemento inﬁnitesimal. 


Este trabajo realizado sobre el ﬂuido debe traducirse en un aumento de su energía total. Así, hay un cambio de 
la energía cinética en la región A1-A′ 
1 que se acaba de abandonar y en la región A2-A′ 
2 que se acaba de ocupar, y 
también un cambio en la energía potencial. En el volumen dV , que al ser el ﬂuido incompresible, es el mismo en las 
dos regiones la masa es ρdV , la energía potencial es gzdm (con z la altura del punto) y la energía cinética 1 


2v2dm. El 
cambio en la energía total queda entonces como 


dU = 
gz2dm + 1 


2v2 
2dm 
− 
gz1dm + 1 


2v2 
1dm 
(3–36) 


y debe ser igual al trabajo realizado sobre el ﬂuido, dU = δW, ya que no se produce ni disipación ni otro tipo de 
intercambio de calor. 


dU = 
gz2dm + 1 


2v2 
2dm 
− 
gz1dm + 1 


2v2 
1dm 
= δW = (p1 − p2)dV = (p1 − p2)dm 


ρ 
(3–37) 


que, al simpliﬁcar queda 


p1 − p2 = ρgz2 + 1 


2ρv2 
2 − ρgz1 − 1 


2ρv2 
1 


o, escrito en la forma habitual, 


p1 + ρgz1 + 1 


2ρv2 
1 = p2 + ρgz2 + 1 


2ρv2 
2 = cte, 
(3–38) 


o también como 
p1 
ρg + z1 + v2 
1 


2g = p2 


ρg + z2 + v2 
2 


2g = cte, 
(3–39) 


expresión conocida como ecuación de Bernoulli, válida a lo largo de una línea de ﬂujo para un ﬂujo que conserva la 
energía, estacionario en un ﬂujo incompresible. 


Ya se ha comentado en la sección anterior el signiﬁcado físico de cada uno de los términos que aparecen en la 
Ecuación de Bernoulli (3–33) o (3–39). 


Cuando el ﬂujo es irrotacional, esta expresión es válida entre dos puntos de un mismo tubo de ﬂujo, y no solamente 
entre puntos de la misma línea de ﬂujo. 


3.5. 
Aplicaciones de las ecuaciones anteriores. 


3.5.1. 
Velocidad de salida por un oriﬁcio (teorema de Torricelli). 


1 


2 


h 


h0 


Aplicando la ecuación de Bernoulli entre los puntos 1 
y 2 


1 
2ρv2 
1 + pa + ρgh = 1 


2ρv2 
2 + pa + ρgh 


Como S1 ≫ S2, de v1S1 = v2S2 se tiene que v1 ≃ 0, 
por lo que la velocidad de salida por un oriﬁcio en un 
depósito abierto a la atmósfera es: 


v2 = 


2gh 
(3–40) 


Obsérvese que la expresión obtenida es válida única- 
mente cuando el depósito está abierto a la atmósfera 
y no cuando el depósito está a presión o hay una capa 
de otro líquido encima de la del líquido que está desaguando. 
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3.5.2. 
Tiempo de vaciado de un depósito. 


A 


a 


y 
y0 


dy 


En un intervalo de tiempo dt, el nivel del líquido en el 
depósito desciende en dy. La variación de volumen será 
entonces: 


−Ady = Q dt = 


2gy a dt 


Agrupando e integrando, 


y 


y0 


dy 
y1/2 = −a√2g 


A 


t 


0 
dt = 
2y1/2y 


y0 = −at 


A 



2g 


de donde queda 


t = A 


a 


2 


g 


y1/2 
0 
− y1/2 
(3–41) 


Para que el depósito se vacíe completamente se necesi- 
taría un tiempo 


tvac = A 


a 


2 


g y0 
(3–42) 


3.5.3. 
Presión dinámica y presión estática. Tubo de Pitot 


El tubo de Pitot se utiliza para medir la velocidad de un ﬂuido. Para ello se utilizan dos aberturas, A y B, perpendi- 
cular y paralela a la dirección del ﬂujo. Si pC es la presión en el gas sobre el ﬂuido (que puede regularse mediante la 
llave), claramente se tiene 
pA = ρg(H + h) + pC 
pB = ρgH + pC 
⇒ 
pA − pB = ρgh 
(3–43) 


En B, la presión será la presión estática de la corriente ﬂuida, pB, mientras que en A el ﬂuido está estancado con 
velocidad 0 y la presión es la presión dinámica. Como ambos puntos están a la misma altura, se tiene 


pA + 0 + 0 = pB + 1 


2ρv2 + 0 


(dicho en otras palabras, a la hora de calcular la altura total en A hay que tener en cuenta el aumento de altura debido 
al ﬂuido que intenta entrar y no puede). Así, se tiene 


v2 = 2(pA − pB) 


ρ 
⇒ 
v = 



2(pA − pB) 


ρ 
. 
(3–44) 


Esto se puede expresar en función de la diferencia de alturas h en la forma 


v = 


2gh 
(3–45) 


Los tubos de Pitot para gases tienen la forma que se muestra en la ﬁgura, y en ellos se utiliza un ﬂuido manométrico 
de densidad ρ′. En este caso, la velocidad del ﬂuido (de densidad ρ) toma la forma 


v = 



2ρ′gh 


ρ 
(3–46) 
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A 
B 


C 


h 


H 
⃗v 


Figura 3–10 Tubo de Pitot. 


A 


B 


h 


⃗v 


Figura 3–11 Tubo de Pitot para gases. 


3.5.4. 
Efecto Venturi. Venturímetro. 


El efecto Venturi es una consecuencia inmediata de la ecuación de continuidad y de la ecuación de Bernoulli. 
Consideremos dos puntos A y B a la misma altura, pero con distintas secciones SA y SB. La ecuación de Bernoulli 
nos da 
pA + 1 


2ρv2 
A = pB + 1 


2ρv2 
B 


que implica, al ser SA > SB 


pB − pC = 1 


2ρ(v2 
B − v2 
C) = 1 


2ρv2 
B 


1 − S2 
B 


S2 
A 


> 0. 
(3–47) 


Como se ve, la presión es más débil donde la velocidad es mayor, es decir, donde la sección es menor. 


Una aplicación inmediata de este efecto es el denominado venturímetro. Consiste en un calibrador colocado en una 
tubería para medir la velocidad de ﬂujo de un líquido. Despejando la velocidad de la ecuación (3–47) se tiene 


v2 
B = 2(pB − pC) 


ρ 
S2 
A 


S2 
A − S2a 


que da 


v2 
A = S2 
B 


S2 
A 
v2 
B = S2 
B 


S2 
A 


2(pB − pC) 


ρ 
S2 
A 


S2 
A − S2a 
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pA 
pB 
pC 


⃗v 
⃗v 


⃗v′ 


Figura 3–12 Efecto Venturi en un estrechamiento. 


1 
2 
A 
a 


⃗h 


ρ 


ρ′ 


Figura 3–13 Relación entre la velocidad y el estrechamiento. 


resultando ﬁnalmente, para la velocidad del ﬂujo 


vA = 



2S2 
B(pA − pB) 


ρ(S2 
A − S2 
B) 
= SB 



ρ′gh2 


ρ(S2 
A − S2 
B) = 
(3–48) 


Otra aplicación de este efecto es la denominada trompa de agua, que se utiliza para hacer bajos vacíos de forma limpia. 
Este mismo efecto se utiliza también en pulverizadores. 
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4.1. 
Introducción. 


A menudo, las fuerzas no conservativas en un ﬂuido no pueden despreciarse. Estas fuerzas dan lugar a una disi- 
pación de la energía mecánica del ﬂuido en energía interna del mismo, igual que las fuerzas de rozamiento disipan 
energía de un bloque deslizante en energía interna del bloque y de la superﬁcie sobre la que se desliza (aumentan sus 
temperaturas). Como consecuencia de ésto, cuando un ﬂuido viscoso circula por una tubería horizontal uniforme, la 
presión a lo largo de una línea de ﬂujo disminuye, al contrario de lo que ocurre en un ﬂuido perfecto (ﬁg. 4–1). Debe 
notarse que al no variar la sección de la tubería, la velocidad del ﬂujo en la misma no varía, al contrario de la presión. 


4.1.1. 
Concepto de viscosidad. 


Todo ﬂuido real posee la propiedad de la viscosidad. La viscosidad es la causa de la resistencia que aparece en el ﬂujo 
de un ﬂuido en tuberías y canales y está relacionada con la oposición de un ﬂuido a las deformaciones tangenciales, 
en las que se produce rozamiento entre distintas partes del ﬂuido. Por ejemplo, a causa de las fuerzas viscosas de 
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Fluido perfecto 


Fluido real 


∆h 


Figura 4–1 La presión disminuye a lo largo del ﬂujo en un ﬂuido real, al contrario de lo que ocure en un ﬂuido perfecto 


laminar 
turbulento 


Figura 4–2 Régimen laminar y régimen turbulento. En régimen laminar el movimiento de las capas del ﬂuido es ordenado, en régimen turbulento 
no. 


rozamiento en una tubería, las capas de ﬂuido más próximas a la pared (más lentas) frenan a las capas más interiores 
(más rápidas), apareciendo una distribución de velocidades. 


El efecto de estas fuerzas de rozamiento se traduce en una disipación de energía que debe tenerse en cuenta en la 
ecuación de Bernoulli (3–33) en la forma siguiente, si el ﬂujo va de 1 a 2, 


p1 
ρg + z1 + v2 
1 


2g − ∆h12 = p2 


ρg + z2 + v2 
2 


2g , 
(4–1) 


donde ∆h12 es la energía disipada entre 1 y 2 por unidad de peso del ﬂuido, a la que se denomina pérdida de carga. 


4.1.2. 
Régimen laminar y régimen turbulento. 


Dependiendo de la importancia relativa entre las fuerzas viscosas y las fuerzas de inercia, el régimen de movimiento 
de un ﬂuido puede ser de dos tipos, 


Régimen laminar o de Poiseuille: El régimen del ﬂujo es laminar cuando el ﬂujo tiene un movimiento ordenado, en 
el que las partículas del ﬂuido se mueven en líneas paralelas, sin que se produzca mezcla de materia entre las 
distintas capas. 


Régimen turbulento o de Venturi: El régimen de movimiento de un ﬂuido es turbulento cuando el ﬂuido presenta 
un movimiento desordenado con mezcla intensiva entre las distintas capas. 


Para distinguir conceptualmente ambos régimenes es conveniente recordar la forma en que experimentalmente se 
observa el paso de un régimen de movimiento a otro. Consideremos una tubería cilíndrica de sección constante por la 
que circula un ﬂuido a una cierta velocidad. Si introducimos en el ﬂujo una aguja ﬁna por la que lentamente inyectamos 
tinta como se muestra en la ﬁgura 4–2 podemos tener dos comportamientos: 


• Cuando el régimen del ﬂujo es laminar, la tinta avanza con el agua como un ﬁno hilo, mezclándose lentamente 
con ésta por procesos de difusión. La mezcla entre el agua y la tinta es muy lenta. Esta situación se tiene para 
diámetros pequeños y velocidades bajas del ﬂuido. 


• A medida que aumenta la velocidad del ﬂujo, llega un momento en que el régimen de movimiento pasa a ser 
turbulento. Cuando el régimen de movimiento del ﬂujo es turbulento, la tinta se mezcla de forma inmediata con 
el agua, diluyéndose en ésta. 
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n 


n 


dn 


τ 


τ + dτ 


x, v 


pared 


β 


⃗v(n) 


v 


v + dv 


Figura 4–3 En un ﬂuido real en movimiento aparecen fuerzas viscosas de rozamiento ⃗Froz entre las distintas capas del ﬂuido, con una fuerza por 
unidad de superﬁcie τ. Esto hace que la velocidad en una sección del ﬂuido no sea uniforme. 


Dependiendo del tipo de ﬂujo, laminar o turbulento, las fuerzas de rozamiento que aparecen serán de distintos tipos. 
En el caso de régimen laminar, las fuerzas que se ejercen entre las distintas capas del ﬂuido son tangentes a la dirección 
del movimiento. Por el contrario, cuando se está en régimen turbulento aparece una segunda contribución a la fuerza 
de rozamiento debida a la mezcla entre las distintas capas. 


4.2. 
Efecto de la viscosidad en los ﬂuidos reales. 


4.2.1. 
Fluidos newtonianos. 


Comencemos considerando el caso en el que el movimiento de un ﬂuido tiene lugar de forma bien ordenada, 
moviéndose éste en líneas paralelas, en lo que acabamos de denominar régimen laminar, teniéndose que las fuerzas de 
rozamiento son tangentes a la dirección del movimiento. Estas fuerzas de rozamiento presentan en general distintos 
comportamientos dependiendo de la rapidez con que varía la velocidad en la dirección perpendicular a la superﬁ- 


cie, rapidez representada por ∂v 


∂n. Una situación muy frecuente se tiene cuando esta relación es de proporcionalidad, 
denominándose ﬂuidos newtonianos a los ﬂuidos que satisfacen esta condición. 


Así, en los ﬂuidos newtonianos, las fuerzas de rozamiento que aparecen entre las capas de ﬂuido en movimiento una 
con respecto a la otra, para una superﬁcie de contacto dS, vienen dadas por la fórmula de Newton, 


dFroz = −µ ∂v 


∂ndS, 
(4–2) 


que en función de la tensión cortante (fuerza por unidad de superﬁcie) queda en la forma 


τ = −∂Froz 


∂S 
= µ ∂v 


∂n, 
(4–3) 


conocida como ley de Newton. En esta expresión ∂v 


∂n es la derivada que caracteriza la variación de la velocidad en 
la dirección transversal al movimiento y la constante de proporcionalidad µ es la viscosidad dinámica o coeﬁciente 
de viscosidad dinámica que depende de las propiedades del ﬂuido, indicando los signos negativos que se trata de una 
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fuerza de rozamiento que se opone al movimiento. Las unidades de la viscosidad dinámica en el sistema internacional 
son [kg m−1s−1] o [N s m−2≡Pa s]. La viscosidad dinámica µ depende mucho de la temperatura y poco de la presión, 
como se trata en más detalle en el apartado 4.2.2. Solamente los ﬂuidos newtonianos (por ejemplo agua, aire, Hg, 
productos petrolíferos típicos, etc...) cumplen esta ley de Newton. Algunos ejemplos de ﬂuidos no newtonianos son 


• Fluidos biológicos. 


• Pinturas al aceite. 


• Soluciones coloidales. 


• Soluciones acuosas de combinados de gran masa molecular. 


La ley de Newton, que en principio fue propuesta como hipótesis para una primera aproximación, está ampliamente 
comprobada experimentalmente en buen número de ﬂuidos y de hecho es la más usual en Mecánica de Fluidos cuando 
se quieren calcular tensiones tangenciales en el movimiento laminar de un ﬂuido. 


Se deﬁne la viscosidad cinemática ν en la forma ν = µ/ρ, y sus dimensiones en el sistema internacional son [ν] = 
m2s−1. La viscosidad cinemática no depende de las características dinámicas del ﬂuido. 


En los sistemas CGS e internacional las unidades de las viscosidades dinámica y cinemática son, respectivamente, 


CGS 
MKS 


µ 
1 poise = 1 dy·s·cm−2 = 1 g/cm·s = 10−1 kg/m·s 
1 kg/m·s = 1 Pa·s = 10 poises 
ν 
1 stokes = 1 cm2·s−1 = 10−4 m2s−1 
1 m2s−1 = 104 stokes 


4.2.2. 
De qué depende la viscosidad 


En general la viscosidad no depende apenas de la presión. Sin embargo, sí depende de la temperatura, y de una forma 
distinta para líquidos y gases. La viscosidad de un líquido disminuye con un aumento de la temperatura, mientras que 
la viscosidad de un gas tiene el comportamiento contrario, aumenta con la temperatura. Esto es debido a los distintos 
orígenes de la viscosidad en ambos casos. 


La movilidad de las moléculas de un líquido es más limitada, con una presencia importante de fuerzas de cohe- 
sión, relacionada con las fuerzas de rozamiento y por tanto con la viscosidad, ya que son estas fuerzas de cohesión 
las que hacen que durante el movimiento del líquido unas capas frenen a otras. Cuando la temperatura aumenta la 
cohesión entre las moléculas disminuye, al aumentar la energía cinética de las mismas. Esto da como consecuencia 
una disminución de la viscosidad. 


La viscosidad de los líquidos se aproxima por la siguiente ley empírica, 


µ = Ae 
−BT 
(4–4) 


donde tanto A como B son magnitudes que se obtienen al ajustar la expresión anterior con los resultados experimen- 
tales de un líquido en particular. 


En el caso de los gases las moléculas tienen una movilidad mucho mayor y el origen de la viscosidad está relacionado 
con el paso de moléculas de diferente energía entre distintas capas. Imaginemos dos capas adyacentes, cada una de 
ellas con una distinta velocidad de movimiento. En promedio las moléculas de la capa más rápida se mueven a mayor 
velocidad que las moléculas de la capa más lenta. Como en el sentido del ﬂujo las moléculas de la capa más rápida se 
mueven en promedio a una velocidad mayor que las de la capa más lenta, el paso de moléculas de la capa más rápida 
a la más lenta se traducirá en una aceleración de la misma, mientras que el paso de moléculas de la capa más lenta a 
la más rápida se traducirá en un frenado de estas últimas, dando este movimiento aleatorio lugar al comportamiento 
viscoso del gas. Como superpuesto a este movimiento se tiene un movimiento aleatorio de las moléculas relacionado 
con la temperatura del gas, pasarán moléculas de la capa más rápida a la más lenta y viceversa. Este paso será tanto 
más rápido cuanto mayor sea la temperatura, y esto dará lugar a una mayor intensidad de la aceleración y el frenado a 
que hacíamos mención en el párrafo anterior, aumentando por tanto la viscosidad con la temperatura. 
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Tabla 4–1 Viscosidades dinámica y cinemática del agua a distintas temperaturas con sus correspondientes densidades. La viscosidad dinámica se 
da en unidades de 10−6 Pa s=10−5 poise, la cinemática en unidades de 10−6 m2/s y la densidad en kg/m3. 


T (◦C) 
Viscosidad 
dinámica 
Densidad 
Viscosidad 
cinemática 
T (◦C) 
Viscosidad 
dinámica 
Densidad 
Viscosidad 
cinemática 


0 
1793 
999,841 
1,793 
60 
469 
983,21 
4,770 
10 
1309 
999,7 
1,309 
70 
406 
977,79 
4,152 
20 
1006 
998,2 
1,008 
80 
357 
971,8 
3,673 
30 
800 
995,65 
0,803 
90 
315 
965,31 
3,263 
40 
657 
992,21 
0,662 
100 
284 
958,35 
2,963 
50 
550 


Tabla 4–2 
Viscosidades dinámica y cinemática de varios ﬂuidos junto con su densidad. La viscosidad dinámica se da en unidades de 
10−6 Pa s=10−5 poise, la cinemática en unidades de 10−6 m2/s y la densidad en kg/m3. 


Fluido 
T (◦C) 
Viscosidad 
dinámica 
Densidad 
Viscosidad 
cinemática 


Mercurio 
0 
1690 
13595 
0,124 
50 
1410 
13472 
0,105 
100 
1220 
13351 
0,091 
Etanol 
20 
1192 
791 
1,506 
Metanol 
20 
591 
810 
0,730 
Benceno 
20 
649 
899 
0,722 
Gasolina 
20 
290 
680-780 
0,427 
Glicerina 
20 
83×104 
1260 
658,7 
Aceite de oliva 
20 
84×104 
920 
913,0 
Aceite lubricante 
20 
26×104 
900-920 
279 
Aire 
20 
18 
1,293 
13,92 
Dióxido de carbono 
20 
16 
1,977 
8,09 
Nitrógeno 
20 
18,4 
1,251 
14,71 
Oxígeno 
20 
20,9 
1,429 
14,62 


La dependencia de la viscosidad de los gases con la temperatura se aproxima por una de las siguientes leyes empí- 
ricas, conocidas como ley de Sutherland 


µ = µ0(T/T0)3/2(T0 + S) 


T + S 
(4–5) 


o por la ley de la potencia 


µ = µ0 
T 
T0 


n 
(4–6) 


En ambos casos µ0 es una viscosidad dinámica conocida a una cierta temperatura, mientras que los valores S y n se 
ajustan con los datos experimentales de un líquido concreto. 


4.2.3. 
La condición de no deslizamiento 


Resumiendo, el esfuerzo cortante en un ﬂuido viscoso newtoniano en régimen laminar es proporcional a la pendiente 
de la velocidad y es máximo en la pared. Además, justo en la pared, la velocidad del ﬂuido es cero con respecto a la 
pared. A este hecho se le denomina condición de no deslizamiento y es una característica de los ﬂuidos viscosos. Un 
ejemplo sencillo de ésto se tiene en las aspas de un ventilador, que tienen polvo ﬁno aunque el ventilador se encuentre 
en movimiento, ya que justo en la superﬁcie del ventilador el aire se encuentra en reposo. Al girar las aspas, éstas 
pueden mover las partículas grandes de polvo pero no las pequeñas. Cuanto más potente sea el ventilador podrá mover 
más partículas, pero siempre quedarán algunas pequeñas que no se mueven. 
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α 


β 


β 


τ 
τ 


dv/dn 
dv/dn 


viscoplástico 


pseudoplástico 


dilatantes 
newtonianos 


Figura 4–4 Reogramas para distintos tipos de ﬂuidos. 


4.2.4. 
Fluidos no newtonianos. 


Acabamos de ver que denominamos ﬂuidos newtonianos a aquéllos en los que la dependencia entre las fuerzas de 
rozamiento y el perﬁl de velocidades es lineal. Sin embargo, no todos los ﬂuidos tienen el mismo comportamiento, 
aunque la mayor parte de los gases y líquidos simples se comportan como ﬂuidos newtonianos. Por ejemplo, pastas, 
lodos y polímeros de alta densidad son ejemplos de ﬂuidos que no tienen un comportamiento newtoniano. Deno- 
minaremos ﬂuidos no newtonianos o anómalos a aquellos ﬂuidos que no siguen la ley de Newton del rozamiento 
interno. 


El comportamiento de los ﬂuidos no newtonianos se representa mediante las curvas de corriente, curvas reológicas† 


o reogramas, como se muestra en la ﬁgura 4–4, donde se representa la tensión cortante producida por las fuerzas de 
rozamiento frente a la tasa de variación de la velocidad en una sección dv/dn. Estas curvas se construyen a partir de 
datos experimentales obtenidos con viscosímetros de distintos tipos en distintas condiciones de funcionamiento. 


Las curvas de corriente correspondientes a los ﬂuidos newtonianos son lineales y se representan en el diagrama por 
rectas que pasan por el origen de coordenadas, con la viscosidad igual a la pendiente de la recta, 


µ = tan α = 
τ 


dv/dn 


que caracteriza completamente las propiedades reológicas del ﬂuido a temperatura y presión dadas, con independencia 
del gradiente de velocidad. 


Las curvas de corriente de los líquidos no newtonianos son muy diversas y en el caso general no lineales. Depen- 
diendo del comportamiento de estos líquidos en el diagrama reológico se tienen distintos comportamientos no newto- 
nianos: 


Fluidos pseudoplásticos: En éstos, al aumentar la velocidad de variación de la velocidad la fuerza cortante viscosa 
aumenta menos que con los newtonianos (véase la ﬁgura). Un ejemplo de materiales pseudoplásticos son 
muchos lodos no newtonianos o fundidos de polímeros. 


Fluidos dilatantes: Los ﬂuidos dilatantes tienen el comportamiento contrario al de los ﬂuidos pseudoplásticos. Muchas 
pinturas son ﬂuidos dilatantes. 


Las curvas de corriente de los ﬂuidos pseudoplásticos y dilatantes se describen muy bien en régimen laminar mediante 
la expresión 


τ = κ 
dv 


dn 


n 
(4–7) 


†Precisamente se denomina Reología a la rama de la Física de Materiales que estudia la respuesta de los ﬂuidos a los esfuerzos. 
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donde κ y n son magnitudes constantes para un líquido dado. κ es una medida del grado de consistencia del líquido 
(teniéndose mayor κ a mayor viscosidad) mientras que n representa el grado de comportamiento no newtoniano. Los 
ﬂuidos newtonianos corresponden al caso particular κ = µ y n = 1, los ﬂuidos pseudoplásticos al caso n < 1 y los 
dilatantes a n > 1. 


Materiales no lineales de Bingham o viscoso-plásticos Otro tipo de ﬂuido no newtoniano corresponde no descrito 
por la ecuación anterior corresponde a los materiales lineales de Bingham o ﬂuidos viscoplásticos. En estos 
materiales el comportamiento es similar al de un sólido para tensiones cortantes por debajo de un determinado 
umbral, es decir presenta una deformación ﬁja para un cierto valor de la tensión cortante. Sin embargo, cuando 
el esfuerzo cortante es mayor que ese umbral, el material presenta el comportamiento de un ﬂuido viscoso 
newtoniano, en la forma 


τ = τ0 + µ dv 


dn 
(4–8) 


donde se conoce a τ como tensión de desplazamiento, tensión inicial de desplazamiento o límite de ﬂuidez. 


En el trabajo común con ﬂuidos no newtonianos se introduce el concepto de viscosidad efectiva, que permite con 
ciertas limitaciones realizar los cálculos mediante las fórmulas de ﬂuidos newtonianos. Para un líquido dado ésta 
no es una magnitud constante, sino que depende de la tasa de cambio de la velocidad dv/dn y de la tensión de 
desplazamiento τ0. Su valor se calcula en el diagrama reológico como la pendiente de las rectas que unen el origen de 
coordenadas con los puntos de la curva de corriente para las condiciones del ﬂujo, 


µe = tan β = 
τ 


dv/dn 
(4–9) 


Hay que insistir en que los valores de la viscosidad efectiva se determinan solamente por el estado instantáneo del 
ﬂujo. 


En los ﬂuidos pseudoplásticos la viscosidad efectiva µe disminuye con el aumento de τ o dv/dn, dando la apariencia 
de que se diluyen al ﬂuir (al ﬂuir aparece una τ y un dv/dn). Por el contrario en los ﬂuidos dilatantes µe aumenta con 
τ o dv/dn y por tanto parece que se espesan al ﬂuir. 


Otro tipo de materiales presentan una combinación de características viscosas y elásticas especial y se les denomina 
ﬂuidos viscoelásticos. Estos materiales se comportan como sólidos ante variaciones rápidas del esfuerzo aplicado, pero 
se comportan como ﬂuidos ante pequeños esfuerzos continuados. Los plásticos a temperatura ambiente sometidos a 
una carga presentan un comportamiento viscoelástico, así como de forma muy espectacular una papilla espesa de 
harina de maíz. 


4.3. 
Ecuaciones de movimiento de los ﬂuidos reales. 


¿Qué forma tienen las ecuaciones del movimiento para un ﬂuido viscoso? En el caso de un ﬂuido viscoso tendremos 
que tener en cuenta que, además de las fuerzas de presión y de las fuerzas de masa como la gravitatoria, sobre un 
ﬂuido viscoso en movimiento actúan las fuerzas viscosas de rozamiento. 


El estado de tensiones en un punto del ﬂuido puede expresarse mediante el tensor de tensiones 


[σij] = 


 


 
 


−p + τxx 
τxy 
τxz 


τyx 
−p + τyy 
τyz 


τzx 
τzy 
−p + τzz 


 


 
 
(4–10) 


donde cada elemento ij del tensor representa la fuerza por unidad de superﬁcie que actúa en la dirección i sobre la 
cara perpendicular al eje j. Esta tensión es en general inhomogénea, es decir, varía de un punto a otro del ﬂujo. Si 
consideramos el elemento de volumen inﬁnitesimal de la ﬁgura, las fuerzas que actúan sobre las caras x y x + dx en 
la dirección x serán 
dFxx(x) = [p(x) − τxx(x)] dydz 
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x 


y 


z 


dx 


dy 


dz 


−p(x) + τxx(x) 


−p(x) + τxx(x) − ∂p 


∂xdx + ∂τxx 


∂x dx 


τxy(y) 


τxy(y) + ∂τxy 


∂y dy 


τxz(z) 


τxz(z) + ∂τxz 


∂z dz 


Figura 4–5 Componente x de las fuerzas que actúan sobre un ﬂuido real en movimiento en las distintas caras de un elemento de volumen. 


sobre la cara que se encuentra a una distancia x del plano yz y 


dFxx(x + dx) = [−p(x + dx) + τxx(x + dx)] dydz 


sobre la cara que se encuentra a una distancia x + dx del plano yz. Si desarrollamos en serie esta última expresión, se 
tiene 


dFxx(x + dx) = 
−p(x) − ∂p 


∂xdx + τxx(x) + ∂τxx 


∂x dx 
dydz 


quedando la fuerza neta que actúa en la dirección x sobre las caras perpendiculares a esta misma dirección 


dF neta 
xx 
= dFxx(x) + dFxx(x + dx) = 
− ∂p 


∂x + ∂τxx 


∂x 


dxdydz 
(4–11a) 


Un cálculo análogo para las fuerzas en la dirección x que actúan sobre las caras perpendiculares a las direcciones y y 
z da 


dF neta 
xy 
= dFxy(y) + dFxy(y + dy) = ∂τxy 


∂y dxdydz 
(4–11b) 


dF neta 
xz 
= dFxz(z) + dFxz(z + dz) = ∂τxz 


∂z dxdydz 
(4–11c) 


quedando la fuerza neta en la dirección x como 


dF neta 
x 
= 
− ∂p 


∂x + ∂τxx 


∂x + ∂τxy 


∂y + ∂τxz 


∂z 


. 
(4–12a) 
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Un cálculo análogo da, para las fuerzas netas en las direcciones y y z 


dF neta 
y 
= 
−∂p 


∂y + ∂τyx 


∂x + ∂τyy 


∂y + ∂τyz 


∂z 


(4–12b) 


dF neta 
z 
= 
−∂p 


∂z + ∂τzx 


∂x + ∂τzy 


∂y + ∂τzz 


∂z 


. 
(4–12c) 


Multiplicando estas ecuaciones respectivamente por los vectores unitarios⃗i, ⃗j y ⃗k y sumando vectorialmente se tiene, 
para la fuerza neta que actúa sobre las distintas superﬁcies 


dF neta 
sup 
= − ∂p 


∂x 
⃗i + −∂p 


∂y 
⃗j + −∂p 


∂z 
⃗k+ 


∂τxx 


∂x + ∂τxy 


∂y + ∂τxz 


∂z 


⃗i+ 


∂τyx 


∂x + ∂τyy 


∂y + ∂τyz 


∂z 


⃗j+ 


∂τzx 


∂x + ∂τzy 


∂y + ∂τzz 


∂z 


⃗k. 


(4–13) 


De una forma más compacta esta expresión se puede escribir en la siguiente forma 


dF neta 
sup 
= −⃗∇p + 
∂ 


∂x, ∂ 


∂y , ∂ 


∂z 


· 


 


 
 


τxx τyx τzx 


τxy τyy τzy 


τxz τyz τzz 


 


 
 
(4–14) 


o en notación tensorial 
dF neta 
sup 
= −⃗∇p + ⃗∇ˆτij. 
(4–15) 


De esta forma, la ecuación de conservación de la cantidad de movimiento para un ﬂuido viscoso se escribe, en su 
forma diferencial, como 


ρd⃗v 


dt = ρg(−⃗k) − ⃗∇p + ⃗∇ˆτij 


= ρ 
∂⃗v 


∂t + (⃗v · ⃗∇)⃗v 
(4–16) 


donde en la parte derecha se ha descompuesto la aceleración en sus componentes local y convectiva (3–23). 


4.3.1. 
Fluido perfecto o ﬂujo no viscoso. 


Obviamente para el caso de un ﬂuido perfecto, todas las componentes del tensor ˆτij son nulas y se recupera la 
ecuación de Euler para ﬂujo no viscoso 


ρd⃗v 


dt = ρg(−⃗k) − ⃗∇p 
(3–17) 


4.3.2. 
Fluido newtoniano. Las ecuaciones de Navier-Stokes. 


Para un ﬂuido newtoniano los esfuerzos viscosos son proporcionales a la derivada de la velocidad y al coeﬁciente 
de viscosidad. Para el caso de ﬂujo incompresible, las componentes del esfuerzo viscoso para un ﬂuido newtoniano 
quedan 


τxx = 2µ∂vx 


∂x 
(4–17a) 


τyy = 2µ∂vy 


∂y 
(4–17b) 


τzz = 2µ∂vz 


∂z 
(4–17c) 
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para las componentes normales, y 


τxy = τyx = µ 
∂vx 


∂y + ∂vy 


∂x 


(4–17d) 


τxz = τzx = µ 
∂vz 


∂x + ∂vx 


∂z 


(4–17e) 


τyz = τzy = µ 
∂vy 


∂z + ∂vz 


∂y 


(4–17f) 


para las cortantes. Sustituyendo estas expresiones en la ecuación anterior se obtienen las ecuaciones de Navier-Stokes 


ρdvx 


dt = ρgx − ∂p 


∂x + µ 
∂2vx 


∂x2 + ∂2vx 


∂y2 + ∂2vx 


∂z2 


(4–18a) 


ρdvy 


dt = ρgy − ∂p 


∂y + µ 
∂2vy 


∂x2 + ∂2vy 


∂y2 + ∂2vy 


∂z2 


(4–18b) 


ρdvz 


dt = ρgx − ∂p 


∂z + µ 
∂2vz 


∂x2 + ∂2vz 


∂y2 + ∂2vz 


∂z2 


(4–18c) 


para un ﬂuido newtoniano. Nótese que éste es un sistema de ecuaciones diferenciales en derivadas parciales no lineales 
de segundo orden que tiene cuatro incógnitas, vx, vy, vz y p. Para resolverlo, necesitaremos una cuarta ecuación que 
no es sino la ecuación de continuidad (3–11) para un ﬂujo incompresible. 


4.4. 
Distribución de velocidades y tensiones cortantes. Fórmula de Poiseuille. 


Determinaremos ahora el caudal o gasto Q en una tubería cilíndrica de radio R como la representada en la ﬁgura 4–6 
en condiciones de ﬂujo laminar. A través de la tubería de sección constante ﬂuye un líquido (homogéneo e incom- 
presible) de viscosidad ν y las presiones en dos puntos 1 y 2 separados una distancia L son respectivamente p1 y 
p2. 


Si consideramos régimen estacionario, al ser v ̸= v(t) no hay aceleración en la dirección axial y la suma de fuerzas 
en esta dirección es cero. Observando el equilibrio de fuerzas en un cilindro imaginario de radio r (distancia al eje del 
cilindro) y longitud L dirigido según el eje mostrado en la ﬁgura, que forma un ángulo θ con la horizontal se tiene, en 
la dirección axial, 
F = 0 ⇒ p1(πr2) − p2(πr2) − τSlat. − mg sen θ = 0. 


Como la masa m es ρSL y L sen θ = z, se tiene 


τSlat. = (p1 − p2)S − ρgzS. 


Podemos expresar la diferencia de presiones p1 − p2 en función de la pérdida de carga ∆h12 a partir de la ecuación 
(4–1), que teniendo en cuenta que v1 = v2 al ser S1 = S2 y que z = z2 − z1 da 


p1 − p2 = ρg∆h12 + ρgz. 


Como la sección de la tubería es constante, la pérdida de carga ∆h12 es igual a la variación de la altura piezométrica 
entre los dos puntos, ∆hpiezo 
12 
= ∆h12, por lo que podemos escribir también 


p1 − p2 = ρg∆hpiezo 
12 
+ ρgz. 


Sustituyendo en la ecuación anterior se van los términos en z y se obtiene la distribución de tensiones cortantes 


τ = ρg∆hpiezo 
12 
r 


2L 
(4–19) 
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r 


⃗v 
z 


R 


L 


mg 


θ 
θ 


p1 


p2 


p1S 


p2S 


τ 


τ 


Figura 4–6 Tubería cilíndrica por la que circula un ﬂuido newtoniano. 


Si el régimen del ﬂujo no fuera laminar habría dos componentes de la tensión cortante, la mostrada arriba y además 
una originada en la mezcla de ﬂuido. En régimen laminar sólo aparece ésta y aplicando a este ﬂujo la ley de Newton,† 


τ = −µ∂v 


∂r = ρg∆hpiezo 
12 
r 


2L 
⇒ − 
v 


vm´ax 
dv = 
r 


0 


ρg∆hpiezo 
12 
rdr 


2µL 


donde vm´ax es la velocidad máxima sobre el eje, quedando 


vm´ax − v = ρgr2 


4µL ∆hpiezo 
12 


Si tenemos en cuenta que para r = R (es decir sobre la pared) debe cumplirse v = 0 se tiene 


vm´ax = ρgR2 


4µL ∆hpiezo 
12 


quedando la ley de distribución de velocidades para un ﬂuido newtoniano en régimen estacionario y laminar en la 
forma parabólica 


v(r) = ρg∆hpiezo 
12 


4µL 
(R2 − r2) 
(4–20) 


coincidiendo el eje de la parábola con el eje de la tubería. 


El gasto o caudal a través de una sección elemental dS en la cual la velocidad es constante viene dado por dQ = vdS, 
siendo el caudal total la integral de los caudales elementales 


¯vS = Q = 


S 
vdS = 
R 


0 
v2πrdr = πR4 


8µLρg∆hpiezo 
12 
(4–21) 


Si comparamos con la expresión anterior para la velocidad máxima se observa que Q = 1 


2πR2vm´ax y si tenemos 
además en cuenta que la velocidad media está relacionada con el caudal a través de Q = ¯vS y que la sección es 
precisamente S = πr2 se tiene de forma inmediata que 


¯v = Q 


S = ρgR2 


8µL ∆hpiezo 
12 
= 1 


2vm´ax 
(4–22) 


para esta distribución parabólica de velocidades. 


†El signo − aparece al tomarse los ejes respecto de r = 0 
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v0 


v0 


v = v(y) 


capa límite 


zona de ﬂuido libre 


Figura 4–7 Consideremos una pequeña parte de la sección de un ﬂujo viscoso alejada de las paredes, en la que podamos considerar la velocidad 
uniforme. Si introducimos una cuchilla en el ﬂujo y observamos la distribución de velocidades en las cercanías de la cuchilla se observan dos 
comportamientos bien diferenciados. Lejos de la cuchilla el efecto de las fuerzas de inercia predomina sobre el de las fuerzas viscosas y se tiene la 
zona de ﬂuido libre, mientras que cerca de la cuchilla el efecto predominante es el de las fuerzas de inercia y se tiene la denominada capa límite. 


4.5. 
Número de Reynolds. Concepto de capa límite. 


Resulta de gran interés conocer la inﬂuencia relativa de las fuerzas de inercia y viscosas en el movimiento del ﬂuido 
y caracterizar esta inﬂuencia. 


El número de Reynolds es un parámetro adimensional que surge a partir de la relación entre las fuerzas de inercia 
y las fuerzas viscosas en las ecuaciones de Navier-Stokes (4–18). El término de inercia ρ(⃗v · ⃗∇)⃗v es el asociado con 
la componente convectiva de la aceleración y dimensionalmente va como ρv2/L, mientras que el término viscoso 
µ∇2v es el asociado con la parte viscosa, y va dimensionalmente como µv/L2, siendo L y v una longitud y una 
velocidad característica respectivamente. Como ambos términos son dimensionalmente equivalentes, su cociente será 
una magnitud adimensional. 


Así, se deﬁne el número de Reynolds en en la forma 


ρv2/L 
µv/L2 
⇒ 
Re = ρ¯vD 


µ 
= ¯vD 


ν 
(4–23) 


donde ¯v es la velocidad media característica del ﬂujo, D es una longitud característica (típicamente el diámetro), ρ es 
la densidad del ﬂuido, µ es la viscosidad dinámica y ν = µ/ρ es la viscosidad cinemática. 


Se denomina velocidad característica a la velocidad por debajo de la cual toda turbulencia es amortiguada por la 
acción de la viscosidad del ﬂuido. Experimentalmente se observa que de forma aproximada para números de Reynolds 
menores que 2300† el régimen de movimiento del ﬂuido es laminar (predominan las fuerzas de rozamiento) mientras 
que para números de Reynolds mayores que 2300 el régimen de movimiento es turbulento (predominan las fuerzas de 
inercia). Experimentalmente se observa el paso de uno a otro régimen introduciendo un colorante en el ﬂujo como se 
mencionó anteriormente. 


Para ﬂuidos de pequeña viscosidad o en general para valores grandes del número de Reynolds, el movimiento del 
ﬂuido suele estudiarse por separado en dos zonas, como se muestra en la ﬁgura 4–7. 


1. Zona de espesor pequeño (capa límite) en la que la inﬂuencia de las fuerzas viscosas es tan importante como la 


†En realidad, el cambio de régimen tiene lugar entre 2000 y 2400, dependiendo de las condiciones concretas en que se realiza el proceso. 
Algunos autores dan otros valores en ese rango como referncia para el cambio de régimen 
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∆h = v2 


2g 


Válvula cerrada 


Válvula abierta (ﬂuido perfecto) 


Válvula abierta (ﬂuido real) 


Figura 4–8 Se comparan los niveles del agua en los piezómetros cuando la válvula está cerrada, y cuando, estando abierta, circula por ella un ﬂuido 
perfecto y un ﬂuido real. 


de las otras fuerzas. 


2. Zona donde la inﬂuencia de la viscosidad es despreciable (zona de ﬂuido libre). 


4.6. 
Pérdidas de carga. Generalidades. 


4.6.1. 
Pérdidas de carga lineales. 


Las pérdidas de carga lineales son aquellas que tienen lugar de forma continua a lo largo del ﬂujo, por el hecho 
de que el ﬂuido avance, aunque no cambien ni la dirección del ﬂujo ni la sección de la tubería. Como un ejemplo, 
consideremos un depósito que desagua a través de un tubo de sección constante en el cual se han colocado una serie 
de piezómetros, como se muestra en la ﬁgura. 


Si la válvula de salida está cerrada, el nivel en todos los piezómetros es el mismo e igual al nivel de agua en el 
depósito. Al abrir la válvula se observan distintos comportamientos dependiendo de si tenemos un ﬂuido perfecto o un 
ﬂuido real. Si lo que tenemos es un ﬂuido perfecto, el nivel de agua en todos los piezómetros baja la misma cantidad, 
y ésta es igual a la altura de velocidad 


∆h = v2 


2g 


Por el contrario, si lo que tenemos es un ﬂuido real, la pérdida de carga se reparte de forma continua a lo largo del 
tubo. Esto hace que el nivel del agua en los piezómetros vaya disminuyendo según se avanza en el tubo. 


Estas pérdidas de carga son, en general, una función compleja de la geometría del sistema, de las propiedades 
del ﬂuido, y de las características del ﬂujo. Experimentalmente se observa que las pérdidas de carga son en general 
proporcionales al cuadrado de la velocidad promedio en la mayor parte de los casos de interés técnico. Esta observación 
se expresa en función de la altura de velocidad en la conocida ley de Darcy-Weiesbach para las pérdidas de carga 
lineales debidas al rozamiento, 


∆hl = f L 


D 
¯v2 


2g . 
(4–24) 


El factor de proporcionalidad f es el coeﬁciente de fricción que en general varía poco para pequeños cambios del ﬂujo 
y es constante para el caso particular de un ﬂujo fuertemente turbulento. 


En general el coeﬁciente de fricción f depende del número de Reynolds, de la rugosidad relativa de la tubería ǫ/D 
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y de la sección de la tubería, 
f = f(Re, ǫ 


D, sección de la tubería) 


donde la rugosidad relativa representa la relación entre la rugosidad de un elemento típico y las dimensiones de la 
tubería representadas por el diámetro de la misma D 


No sólo es importante la magnitud, sino también la forma de la sección de la canalización, ya que diferencias de la 
forma circular dan lugar a ﬂujos secundarios que aumentan la caída de presión. Cuando se trabaja con canalizaciones 
no circulares se utiliza generalmente el radio hidráulico o el diámetro hidráulico 


Rh = 2A 


per. = 2 × área de la sección 


perímetro mojado 
Dh = 2Rh = 4A 


per. 
(4–25) 


en vez del radio o el diámetro, y se trata la canalización como si fuera de sección circular, con el número de Reynolds 
calculado en función del diámetro hidráulico. 


Régimen laminar. 


Para tuberías de sección circular uniforme por las que circula un ﬂuido newtoniano en régimen laminar, esta pérdida 
de carga que aparece tanto en tuberías lisas como rugosas puede deducirse de la fórmula de Poiseuille. A partir de la 
expresión para el caudal (4–21) y teniendo en cuenta que el diámetro es dos veces el radio, se tiene que 


∆hl = 8µLQ 


ρgπR4 = 8µL¯v 


ρgR2 = 32µL¯v 


ρgD2 = 32νL¯v 


gD2 
(4–26) 


Si tenemos en cuenta la deﬁnición del número de Reynolds, Re = ¯vD/ν, esta expresión queda como 


∆hl = 32νL¯v 


gD2 
= 64 


Re 
L 
D 
¯v2 


2g = f L 


D 
¯v2 


2g 
(4–27) 


donde f = 64/Re es el coeﬁciente de fricción (adimensional), L es la longitud, D el diámetro y ¯v2 


2g la altura de 


velocidad. Esta expresión coincide con la ley empírica (4–24) de Darcy-Weiesbach 


∆h = f L 


D 
¯v2 


2g 
(4–24) 


que se utiliza para un ﬂujo general, tomando para el caso particular del ﬂujo laminar, la relación f = 64/Re, conocida 
también como ecuación de Hagen-Poiseiuille. 


Estas relaciones permiten obtener la tensión cortante en la pared de la tubería circular para un ﬂujo laminar. De las 
expresiones anteriores (4–19), se tiene 


τ0 = ∆hl 
γR 
2L = ∆hl 
γR 
2L 
Al despejar ∆hl se obtiene 


∆hl = 2τ0L 


γR = 4τ0L 


γD 
que, al comparar con la ecuación de Darcy-Weiesbach da la tensión en las paredes de la tubería 


τ0 = f γ¯v2 


8g = f ρ¯v2 


8 
(4–28) 


Cuando las canalizaciones no son cilíndricas la relación entre el número de Reynolds y el coeﬁciente de fricción es, 
para régimen laminar y ﬂuido newtoniano, de la forma 


f = k 


Re 
con 48 ≤ k ≤ 96 
(4–29) 


en función de un factor geométrico k. Algunos valores típicos del factor geométrico k son los siguientes 
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Cuadrado 
56,91 
Rectángulo 2:1 
62,19 
Rectángulo 5:1 
76,28 
Láminas paralelas 
76,28 


Régimen turbulento 


La ley experimental de Darcy-Weiesbach es válida para todo tipo de ﬂujo con el coeﬁciente de fricción adecuado. 
Como acabamos de ver, para un ﬂuido newtoniano en régimen laminar y tuberías circulares, el coeﬁciente de fricción 
viene dado por f = 64/Re. 


Tuberías lisas. 
En régimen turbulento (Re >2300) y tuberías lisas se utilizan diversas expresiones f(Re) depen- 
diendo del rango del número de Reynolds. Para números de Reynolds en el rango (3000<Re<100000) se utiliza la 
fórmula de Blausius 
f = 0,316/Re0,25 
(4–30) 


mientras que para números de Reynolds superiores a éstos, pero menores que 3 × 106 se utiliza la fórmula de von 
Karman 
1 
√f = 0,869 ln 
Re 
√f 


− 0,8 = 2 log 
Re 
√f 


− 0,8 
(4–31) 


Tuberías rugosas 
En régimen turbulento y tuberías rugosas se utilizan distintas funciones del tipo f = f(ǫ/D, Re) 
donde ǫ es el coeﬁciente de rugosidad y d el diámetro de la conducción. 


Para régimen turbulento en tuberías circulares la más utilizada de todas es la fórmula de Colebrook (1939) 


1 
√f = −0,869 ln 
ǫ/D 


3,7 + 2,523 


Re√f 


= −2 log 
ǫ/D 


3,7 + 2,523 


Re√f 


. 
(4–32) 


Ésta es una ecuación implícita que converge bien en pocas iteraciones cuando se aborda mediante métodos numéricos. 
El diagrama de Moody que se presenta en la ﬁgura 4–9 es una representación logarítmica de la ecuación de Colebrook 
en ejes coeﬁciente de fricción/número de Reynolds para distintos valores de la rugosidad relativa, combinado con la 
relación f = 64/Re para ﬂujo laminar. 


En el intervalo de rugosidades relativas 10−6 ≤ ǫ/D ≤ 10−2 y de números de Reynolds 5000 ≤ Re ≤ 108 se 
puede hacer la aproximación explícita 


f = 
1, 325 


ln 
ǫ/D 


3,7 + 5,74 


Re0,9 


2 = 
0,25 


log 
ǫ/D 


3,7 + 5,74 


Re0,9 


2 
(4–33) 


que da valores aproximados de la ecuación de Colebrook en buena parte del rango de interés. 


Para rugosidades relativas altas y números de Reynolds altos, f depende fundamentalmente de la rugosidad relativa 
y poco del número de Reynolds, como se observa en el diagrama de Moody y puede ajustarse por 


1/ 


f = 1,14 − 0,869 ln(ǫ/D) = 1,14 − 2 log(ǫ/D). 
(4–34) 
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Figura 4–9 Diagrama de Moody para tuberías de sección circular obtenido a partir de la relación f = 64/Re para régimen laminar y la ecuación 
de Colebrook (4–32). 


f0 
f 


∆h 


v 


Re 


D. Moody 


Cuando se conocen las características del ﬂujo y de la tubería es sencillo 
obtener el coeﬁciente de fricción f y la pérdida de carga directamente desde 
el diagrama de Moody. Cuando la situación es la contraria, es decir se cono- 
cen las características de la tubería y de alguna forma la pérdida de carga (p.ej, 
midiendo las presiones entre dos puntos si se conoce el diámetro de la tubería) 
y se quiere obtener el caudal, no es posible hacer el proceso inverso salvo que 
se conozca f para esas características del ﬂujo (f depende tanto de las carac- 
terísticas de la tubería como de las características del ﬂujo). Para resolver el 
problema, se utilizan métodos iterativos. Se parte de un valor de f estimado (p.ej, si no se tiene más información, algo 
en torno a 0,02), se calcula la velocidad a partir de la ecuación de Darcy-Weiesbach. De aquí y de las características 
de la tubería se obtiene el número de Reynolds y con la rugosidad relativa se extrae del diagrama de Moody un nuevo 
valor del coeﬁciente de fricción. Repitiendo el proceso hasta que converja suﬁcientemente se tiene un buen valor de f 
y el caudal como Q = vS. 


4.6.2. 
Pérdidas de carga singulares. 


Representan disipaciones adicionales de energía en el ﬂujo, normalmente causadas por los ﬂujos secundarios (turbu- 
lencias) inducidos en cambios de sección o de dirección de la tubería. Aunque en algunos sistemas constituyen una 
parte importante de la pérdida de carga total, se les suele denominar denominar pérdidas secundarias o menores. 
Corresponden a: 
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• 
Secciones de rejilla. 
• 
Ensanchamientos. 
• 
Estrechamientos. 
• 
Codos. 
• 
Válvulas. 
• 
Bifurcaciones. 


Para obtener las pérdidas de carga locales o singulares, en general se utilizan expresiones experimentales de la forma 


∆h = K v2 


2g 
(4–35) 


donde K es un coeﬁciente que se obtiene a veces de forma teórica y a veces de forma experimental y que se encuentra 
en los manuales de hidráulica. En el caso concreto de los ensanchamientos esta expresión aparece en función de (∆v)2 


(que en el fondo es similar a 4–35). 


No hay que olvidar que, aunque el coeﬁciente K parece constante, en general tiene alguna dependencia adicional 
de las condiciones del ﬂujo. En particular, los siguientes factores afectan a K 


• La geometría exacta del componente en cuestión. 


• El número de Reynolds. 


• Proximidad a otros elementos, etc... (Los valores tabulados son para componentes “aislados”, con una cierta 
longitud de tubería entre el elemento analizado y el siguiente). 


Por ejemplo, para los ensanchamientos graduales la dependencia de K con las secciones inicial y ﬁnal del tubo es tal 
que la pérdida de carga viene dada por† 


∆h = K (∆v)2 


2g 
(4–36) 


donde en K aún queda incluida una dependencia en el ángulo del cono. 


A continuación se da un resumen de algunos de los coeﬁcientes de pérdidas más habituales. En todas las expresiones 
que siguen se ha utilizado el subíndice 1 para representar los valores aguas arriba y 2 para representar la situación aguas 
abajo. Para enfatizar a que parte corresponden las velocidades se ha intentado utilizar las mayúsculas para los valores 
correspondientes a la sección mayor y las minúsculas para los valores correspondientes a la sección menor. 


1. 
De depósito a tubería (pérdida a la entrada) 
Conexión a ras de la pared 
tubería entrante 
conexión abocinada 


0,50 v2 
2 


2g 
1,00 v2 
2 


2g 
0,05 v2 
2 


2g 


2. De tubería a depósito (pérdida a la salida) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1,00 v2 
1 


2g 


3. Ensanchamiento brusco . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (v1 − V2)2 


2g 


4. Ensanchamiento gradual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K (v1 − V2)2 


2g 


†Esto no es más que un caso particular de 4–35, ya que por la ecuación de continuidad, v1 S1 = v2 S2 y por tanto, 


∆h = K (∆v)2 


2g 
= K (v2 − v1)2 


2g 
= K 
1 − S2 


S1 


2 v2 
2 


2g 
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ángulo | d1/d2 
1,2 
1,4 
1,6 
1,8 
2,0 
2,5 
3,0 
4,0 
5,0 


4◦ 
0,02 
0,03 
0,03 
0,04 
0,04 
0,04 
0,04 
0,04 
0,04 
10◦ 
0,04 
0,06 
0,07 
0,07 
0,07 
0,08 
0,08 
0,08 
0,08 
15◦ 
0,09 
0,12 
0,14 
0,15 
0,16 
0,16 
0,16 
0,16 
0,16 
20◦ 
0,16 
0,23 
0,26 
0,28 
0,29 
0,30 
0,31 
0,31 
0,31 
30◦ 
0,25 
0,36 
0,42 
0,44 
0,46 
0,48 
0,48 
0,49 
0,50 
50◦ 
0,35 
0,50 
0,57 
0,61 
0,63 
0,65 
0,66 
0,67 
0,67 
60◦ 
0,37 
0,53 
0,61 
0,65 
0,68 
0,70 
0,71 
0,72 
0,72 


5. Venturímetros, boquillas y oriﬁcios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 
1 


cr 
− 1 
v2 
2 


2g 


6. Contracción brusca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K v2 
2 


2g 


D1/d2 
1,2 
1,4 
1,6 
1,8 
2,0 
2,5 
3,0 
4,0 
5,0 


K 
0,08 
0,17 
0,26 
0,34 
0,37 
0,41 
0,43 
0,45 
0,46 


7. Codos, accesorios, válvulas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K v2 


2g 


codo a 45◦ . . . . . . . . . . . . . . . . . . 
0,35 − 0,45 
codo a 90◦ . . . . . . . . . . . . . . . . . . 
0,50 − 0,75 
bifurcaciones en T . . . . . . . . . . . 
1,50 − 2,00 
válvulas de compuerta (abierta) 
≃ 0,25 
válvulas de control (abierta). . . 
≃ 3,0 


8. Salida de tubería a la atmósfera (desagüe libre a la atmósfera). 


La pérdida en el desagüe libre de una tubería a la atmósfera es generalmente muy pequeña al no producirse 
apenas turbulencias en la tubería en ese proceso. 


4.6.3. 
Variación de altura debida a una máquina. 


Las máquinas hidráulicas pueden 


• Aportar energía al ﬂujo (bombas). 


• Extraer energía del ﬂujo (turbinas). 


El aporte de energía que realiza una bomba al ﬂujo se representa habitualmente a través de la energía aportada por 
unidad de peso del ﬂuido, a la que se denomina altura de la bomba o altura aportada por la bomba HB. Veamos ahora 
la relación entre la potencia de la bomba y la altura de la misma. Para ello consideraremos una bomba colocada en 
una tubería horizontal de sección constante y estudiaremos el aumento de energía en un elemento dm del ﬂuido que 
atraviesa la bomba en un intervalo de tiempo dt. Si denominamos h1 a la energía por unidad de peso a la entrada de 
la bomba y h2 a la energía por unidad de peso a la salida de la bomba, la energía de un elemento de masa dm antes 
de entrar en la bomba será gh1dm mientras que la energía del elemento de masa dm al salir de la bomba será gh2dm, 
por lo que la variación de energía del ﬂuido en el intervalo de tiempo dt es g(h2 − h1)dm = gHBdm donde HB es 
la altura de la bomba o energía por unidad de peso aportada por la bomba. La potencia de la bomba será entonces, en 
función de la altura de la misma 


P = gHBdm 


dt 
= gHB 
dm 
dt = gHBQm = ρgQHB 
(4–37) 
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1 


2 


h1 


h2 


Figura 4–10 


A 
A 
B 
B 


Bomba 
Turbina 


pA + ρgHbomba = pB 
pA − ρgHturbina = pB 


Figura 4–11 Expresiones para incluir correctamente el efecto de las bombas y de las turbinas. 


La misión de la turbina es extraer energía del ﬂuido, disminuyendo su energía. Análogamente al caso de las bombas, 
se denomina altura HT de la turbina a la energía que ésta extrae por unidad de peso del ﬂuido. De forma semejante al 
caso de la bomba, la potencia extraída por la turbina está relacionada con la altura HT de ésta mediante 


Ptur = ρgHT Q = gHT Qm 
(4–38) 


Para considerar correctamente los signos de las alturas de bomba y turbina, es conveniente seguir el sentido del ﬂujo, 
añadiendo a la altura total inicial todas las alturas aportadas por bombas y restándole todas las alturas de pérdidas y 
todas las alturas extraídas por turbinas, llegándose a la ecuación de Bernoulli generalizada. 


p1 
ρg + v2 
1 


2g + h1 + 
HB − 
HT − ∆htot = p2 


ρg + v2 
2 


2g + h2 
(4–39) 


4.6.4. 
Representación gráﬁca de la pérdida de carga. 


Consideremos el caso de una tubería alimentada por un gran depósito cuyo nivel se mantiene constante, que descarga 
en otro gran depósito cuyo nivel se mantiene también constante. Las líneas con que se trabaja habitualmente son la 
línea de alturas totales y la línea de alturas piezométricas. 


1. La línea de altura total en cada punto es el resultado de restar a la del depósito (nivel de energía del depósito) 
todas las pérdidas de carga encontradas en el camino (en este caso todas las pérdidas de carga a la izquierda del 
punto) y sumar todas las ganancias añadidas por las bombas. 


2. La línea de altura piezométrica se obtiene restando la altura de velocidad v2 


2g de la línea de altura total 
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nivel 


Htotal 


Hpiezo 


∆Ht = 
i ∆hi 


Figura 4–12 Línea de alturas totales y línea de alturas piezométricas. 


4.7. 
Redes de tuberías 


Analizaremos ahora algunos ejemplos sencillos de redes de tuberías y los elementos fundamentales a tener en cuenta 
para su cálculo. 


El caso más sencillo es el de tuberías en serie. En este caso, en régimen permanente el caudal másico Qm que circula 
a través de cada tubería es el mismo, Qm1 = Qm2... = Qmn, mientras que la pérdida de carga total a lo largo del 
sistema de tuberías en serie se obtiene como la suma de las pérdidas de carga en cada uno de los tramos (más las 
pérdidas en los elementos de unión de las tuberías), 


∆htot = 


i 
∆hi 


En las tuberías en paralelo, la corriente que originalmente circulaba por una tubería se distribuye en un conjunto de 
tuberías para posteriormente reunirse de nuevo en una sola. En este caso el caudal en la línea principal es la suma de 
los caudales en las tuberías en que se divide, Qprincipal = Qm1 + ... + Qmn. Como los valores de la altura de carga 
en los puntos de bifurcación deben ser los mismos para todas las tuberías colocadas en paralelo, las pérdidas de carga 
en todas ellas deben ser iguales. 


Cuando tenemos redes de tuberías más complejas hay dos reglas fundamentales que deben tenerse en cuenta: 


A 


B 


C 


Q1 


Q1 


Q2 


Q2 
Q3 


Q3 


Figura 4–13 Reglas fundamentales para trabajar con redes de tuberías 


1. En cualquier nodo (unión de dos o más tuberías) debe conservarse la masa. Es decir, el caudal másico neto debe 
ser cero en régimen estacionario, todo lo que entra en un intervalo de tiempo debe salir en el mismo intervalo, 
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no puede producirse acumulación de materia. Para el caso de un ﬂuido incompresible, ésto es cierto también 
para el caudal. 


i 
Qmi = 0 


2. No debe de haber pérdida de carga neta a lo largo de un bucle. Si comenzamos en una posición del bucle y 
sumamos las variaciones de carga hidráulica, cada una con su correspondiente signo, al llegar a la posición 
original la carga hidráulica debe ser la misma que al principio, 


i 
∆hrozi + 
∆Hbombas + 
∆Hturbinas = 0 


en el bucle. 


4.8. 
Cavitación 


Se denomina cavitación a la formación de cavidades en el líquido en movimiento llenas de vapor. Se produce cuando 
la presión en algún lugar del ﬂuido se hace menor que la presión de vapor saturante del líquido, es decir la presión que 
corresponde a la ebullición del líquido a una temperatura dada. 


Por ejemplo, acelerando agua desde el reposo hasta velocidades de 15 m/s la presión desciende en aproximadamente 
una atmósfera y produciéndose cavitación. De hecho esto ocurre a menudo en los elementos en los que se producen 
pérdidas locales cuando se produce un gran aumento de velocidad al estrecharse el ﬂujo después de dilatarse. 


Así, en las partes de baja presión puede comenzar una intensa evaporación de líquido, formándose las cavidades de 
cavitación formadas por burbujas de vapor. Al aumentar posteriormente la presión el vapor se condensa, normalmente 
con un brusco chasquido y las cavidades de cavitación desaparecen. 


La cavitación va acompañada de un notable aumento de las pérdidas locales y se caracteriza a través un parámetro 
adimensional denominado número de cavitación 


Ca = 2(p − pv) 


ρv2 
(4–40) 


donde p y v son la presión y la velocidad antes de la resistencia local pv la presión de vapor saturante del agua. 


Para que en el elemento local no se produzca cavitación el número de cavitación debe ser menor que un valor crítico 
que se determina de forma experimental, Ca < Cacr.. 


También se produce cavitación en hélices, bombas y turbinas. La cavitación produce grandes oscilaciones, con 
aumento de las pérdidas de carga y además ayuda a la corrosión de los metales. Ésto último es consecuencia de los 
procesos casi instantáneos de contracción de las burbujas que producen grandes presiones que inicialmente desmigajan 
lentamente el material y posteriormente extienden su efecto al interior del mismo, combinados a menudo con efectos 
químicos. Como consecuencia de todo ésto el metal se corroe y llega a adquirir una estructura esponjosa. 


4.9. 
Golpe de ariete o choque hidráulico 


Se denomina golpe de ariete o choque hidráulico al brusco aumento de la presión que se produce en las tuberías 
cuando se detiene bruscamente el líquido que circula por ellas. Ésto ocurre, por ejemplo, cuando se cierran rápidamente 
alguno de los distintos dispositivos (válvulas, grifos) presentes en la tubería o cuando se produce una parada imprevista 
de las bombas que impelen el líquido. 


El golpe de ariete es especialmente peligroso en las tuberías largas por las que circulan grandes cantidades de líquido 
a altas velocidades. En estos casos, el golpe de ariete puede producir graves daños en las juntas de unión, en la propia 
tubería o en las máquinas. 


Para estudiar este fenómeno consideremos una tubería horizontal de sección circular constante S y radio R por la que 
circula un líquido de velocidad v constante. Si se cierra la válvula de forma idealmente instantánea la capa de líquido 
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⃗v 


p0 
p0 + ∆p 


dl 


Figura 4–14 Golpe de ariete. El líquido que circula por la tubería se detiene súbitamente. 


en contacto con esta se detendrá de forma inmediata produciéndose un aumento de la presión como consecuencia de 
la transformación de energía cinética en energía potencial de presión. 


Consideremos primero un caso extremo, el caso de un ﬂuido completamente incompresible que circula por una 
tubería totalmente rígida. En este caso, si la válvula se cierra en un tiempo ∆t todo el ﬂuido se detiene simultáneamente 
en ese tiempo la bajo la acción de la fuerza ∆pπR2 en la superﬁcie en contacto con la válvula y la cantidad de 
movimiento de éste cambia de ρLπR2v a 0 en el intervalo ∆t. El impulso de la fuerza mencionada debe ser igual a la 
variación de la cantidad de movimiento del líquido, 


∆pπR2∆t = ρLπR2 


Esta ecuación da un buena estimación del aumento de presión asociado al cierre lento de una válvula. Sin embargo 
da un resultado pequeño y no demasiado interesante. También hay que recalcar que da como resultado una presión 
inﬁnita si el cierre de la válvula se produce de forma instantánea, pero para presiones muy altas las aproximaciones de 
líquido totalmente incompresible y tubería totalmente rígida son malas. 


Cuando la válvula se cierra de forma instantánea el ﬂuido se frenará rápidamente y la presión aumentará, pero este 
aumento de presión producirá una pequeña compresión en el ﬂuido y una pequeña expansión de la tubería (que espera- 
mos elástica). Así, el líquido en contacto con la válvula se detendrá instantáneamente, pero al ser éste compresible, el 
resto del líquido en la tubería no se detendrá de forma instantánea. El límite del volumen que incluye el líquido parado 
se desplaza alejándose de la válvula con una cierta velocidad c a la que denominamos velocidad de propagación de la 
onda de presión. 


Examinemos el elemento de volumen del líquido Scdt = Sdl junto a la válvula. En un tiempo dt este volumen se 
ha detenido y su cantidad de movimiento ha pasado de ρSdl v a 0. El impulso de la fuerza que actúa en el tiempo dt 
deteniendo el ﬂuido es ∆pSdt. Al igualarlo a la variación de la cantidad de movimiento del elemento en ese tiempo 
dt tenemos 
∆pSdt = ρSdl v 


Como c = dl/dt se obtiene el valor del aumento de presión en el golpe de ariete: 


∆p = ρcv 
(4–41) 


Después se detiene la siguiente capa y así sucesivamente, y la presión elevada que inicialmente surge en las proxi- 
midades de la válvula se propaga poco a poco por toda la tubería con una velocidad c, en contra del ﬂujo original del 
líquido. 


Si la presión al principio de la tubería no cambia (por ejemplo, cuando el líquido viene de un depósito abierto a la 
atmósfera) una vez que la onda de presión alcanza la sección inicial de la tubería, ésta comienza a propagarse en el 
sentido contrario con la misma velocidad c, ahora como onda de disminución de presión. Cuando esta onda alcanza la 
válvula la presión cae y se hace menor que antes del cierre y de nuevo una onda de presión sale en sentido contrario. 
Así, se alternan ciclos de aumento y disminución de la presión en intervalos de tiempo iguales al tiempo empleado por 
la onda en recorrer completamente en los dos sentidos el tramo de tubería entre la válvula y el principio. 
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Así, el líquido realiza un movimiento oscilatorio que se va amortiguando como consecuencia de la disipación de 
energía que tiene lugar en el proceso. 


La velocidad de propagación de la onda de presión depende de las características del líquido y de la tubería y viene 
dada por la expresión:† 


c = 



K 
ρ 



1 + KD 


E∆R 


−1 


(4–42) 


donde K es el módulo de compresibilidad del líquido (inversa del coeﬁciente de compresibilidad), ρ la densidad del 
líquido, E el módulo de Young del material de la tubería, D el diámetro de la misma y ∆R el espesor de la pared de 
la tubería. 


Si se considera totalmente la tubería totalmente rígida (E = ∞) la velocidad de propagación de la onda de presión 
queda en la forma 


c = 



K 
ρ , 
(4–43) 


igual a la velocidad de propagación del sonido en un líquido. 


Para evitar los efectos del golpe de ariete se utilizan diversos sistemas: 


1. Dispositivos que aumentan el tiempo de cierre de todo tipo de válvulas aminorando así el efecto. 


2. Uso de un diámetro grande en la tubería para que las velocidades sean bajas. 


3. Dispositivos de seguridad en las bombas para que en caso de corte de energía reduzca lentamente la velocidad 
del motor. 


4. Colocación de válvulas de seguridad automáticas, campanas de aire y chimeneas de gran sección. 


†que se obtiene de la condición de igualdad entre la energía cinética del líquido que se mueve por la tubería y la suma de los trabajos de 
compresión del líquido y alargamiento del tubo. 
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5.1. 
Fenómenos de superﬁcie: 
tensión superﬁcial y capilaridad. 


5.1.1. 
Tensión superﬁcial. 


Las superﬁcies libres de los líquidos presentan una serie de fenómenos cuyo origen se reduce a la tendencia de la 
superﬁcie a permanecer lo más pequeña posible. Este hecho, conectado con la presencia de una tensión en la película 
superﬁcial de forma análoga a la tensión en una membrana elástica ﬁna que se encuentra tensa puede observarse de 
forma sorprendente mediante una serie de experimentos simples. 


Algunos ejemplos de tensión superﬁcial. 


Consideremos un bastidor como el de la ﬁgura 5–1, que se sumerge en una solución jabonosa, quedando al sacarlo 
un ﬁna película de la misma. Si sobre esta película depositamos un hilo cerrado en un lazo, se observa que éste no 
adopta una posición especial, sino la posición que haya tomado al caer. Sin embargo, si se pincha la parte de la película 
dentro del hilo, éste adopta la forma de una circunferencia por la acción de la tensión superﬁcial. De esta forma, la 
película de ﬂuido ha alcanzado su mínima superﬁcie. 


Otro efecto de la tensión superﬁcial es el siguiente: ¿Qué ocurre cuando un objeto de mayor densidad que el agua se 
deposita sobre la misma? La respuesta más inmediata sería que éste se hundiría. Sin embargo, ésto no siempre es así. 
Consideremos el caso de una aguja sobre el agua. Si ésta se deposita con el suﬁciente cuidado es posible hacerla ﬂotar 
sobre la superﬁcie del líquido. Esto es debido a que las fuerzas que soportan la aguja son debidas no sólo al empuje 
del agua sobre la misma, que es obviamente menor que su peso, sino también a la tensión superﬁcial. 
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Figura 5–1 Efecto de la tensión superﬁcial en una película delgada. 


Otro efecto de la tensión superﬁcial es el hecho de que un líquido que no se encuentre sometido a fuerzas externas 
tiende a tener una forma esférica, minimizando de esta forma su superﬁcie para un volumen dado. 


La tensión superﬁcial es utilizada por algunos insectos como el zapatero (gerris lacustris) que se muestra en la 
ﬁgura 5–2. Éste se mantiene muy por encima del agua gracias a la tensión superﬁcial en sus patas, que utiliza para 


Figura 5–2 El zapatero se desliza sobre la superﬁcie del agua gracias a la tensión superﬁcial. Fotografía de ISAKA Yoji bajo licencia CC-BY-SA. 
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Tabla 5–1 
Valores de la tensión superﬁcial medidos en dy/cm. En todos los casos se asume que el líquido está en contacto con aire. 
(1dy/cm=10−3Nm−1). Salvo indicación en sentido contrario, los datos están tomados a 20◦C 


Sustancia 
tensión 
superﬁcial 
Sustancia 
tensión 
superﬁcial 


Acetona 
23,7 
Agua a 5◦C 
74,92 
Metanol 
22,6 
Agua a 10◦C 
74,22 
Etanol 
22,3 
Agua a 15◦C 
73,49 
Aceite de oliva 
33 
Agua a 20◦C 
72,75 
Aceite de paraﬁna 
26 
Agua a 25◦C 
71,97 
Benceno 
28,9 
Agua a 30◦C 
71,18 
Glicerina 
64 
Mercurio 
475 


remar (las centrales) o a modo de timón (las posteriores), mientras que caza a sus presas con las delanteras, pero todas 
ellas contribuyen a que se mantenga por encima del agua. Obsérvese la curvatura de la superﬁcie del agua bajo las 
patas. 


Principios físicos de la tensión superﬁcial. 


Cualitativamente ésto tiene su origen en la fuerza de atracción molecular no compensada en las moléculas de la 
superﬁcie. Mientras que en el interior de un líquido, cada molécula está rodeada de otras moléculas por todos lados, 
y por tanto las fuerzas intermoleculares se compensan unas a otras, en la superﬁcie no hay moléculas por encima de 
las moléculas superﬁciales. Así, las moléculas de la superﬁcie experimentan una fuerza neta dirigida hacia el interior 
de modo que éstas no escapan. De este modo la superﬁcie intenta hacerse lo menor posible. Asimismo, si se eleva 
ligeramente una molécula superﬁcial, los enlaces moleculares con las moléculas adyacentes se alargan, produciéndose 
una fuerza restauradora que tira de la molécula desplazada de nuevo hacia la superﬁcie. 


Esto es lo que ocurre cuando se coloca cuidadosamente una aguja sobre la superﬁcie. Las moléculas superﬁciales 
descienden ligeramente, separándose de las demás moléculas de la superﬁcie, que ejercen una fuerza atractiva sobre 
las mismas, soportando la aguja. Así, la superﬁcie de un líquido presenta una cierta semejanza con un membrana 
elástica tensa. 


Este concepto de tensión en la superﬁcie de un ﬂuido, conocido como tensión superﬁcial fue introducido por J. A. 
von Segner como la fuerza que actúa en la superﬁcie de un ﬂuido a cada lado de una línea de longitud unidad dibujada 
en la superﬁcie. 


La fuerza necesaria para romper la superﬁcie se puede medir levantando un alambre ﬁno. La fuerza necesaria para 
romper la superﬁcie resulta ser proporcional a la longitud de la superﬁcie que se rompe, que es el doble de la longitud 
del alambre, si hay una película de superﬁcie en ambos lados del mismo, o a la longitud del mismo, si sólo hay 
película superﬁcial en un lado. Si el alambre posee una masa m y una longitud l, y es necesario ejercer una fuerza F 
para separarlo de la superﬁcie (de dos caras), el coeﬁciente de tensión superﬁcial σ es: 


σ = F − mg 


2l 
(5–1) 


Algunos valores de éste coeﬁciente de tensión superﬁcial se dan en la tabla 5–1. 


Relación entre la tensión superﬁcial y la diferencia de presiones a través de una superﬁcie. 


Estas no son las únicas consecuencias de la tensión superﬁcial. Hasta el momento se ha dicho que, una vez alcanzado 
el estado de equilibrio, las presiones de dos cuerpos en contacto han de ser iguales. Sin embargo, esta aﬁrmación es 
cierta únicamente cuando el efecto de los fenómenos de superﬁcie es despreciable. Cuando se tiene en cuenta la tensión 
superﬁcial, las presiones de dos medios en contacto pueden ser diferentes. 
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α/2 
α/2 


α/2 


r 


r1 


r2 


dα1 


dα2 


dl1 
dl2 
σdl1 


σdl1 


σdl2 


σdl2 


Figura 5–3 Fuerza ejercida por la tensión superﬁcial sobre un elemento inﬁnitesimal rectangular de una superﬁcie curva. 


Como ejemplo de lo que se acaba de aﬁrmar, consideremos el caso de una gota de líquido en el aire. La tendencia 
de ésta a disminuir su superﬁcie conduce a una compresión de la gota y a un aumento de la presión interna, siendo la 
presión del líquido en la gota mayor que la del aire circundante. A la diferencia entre ambas presiones se le denomina 
diferencia, variación de presión o presión diferencial y la designaremos por pdif. 


Para calcular esta magnitud, debemos tener en cuenta que el trabajo realizado por las fuerzas superﬁciales al dismi- 
nuir la superﬁcie de la gota en la magnitud dS viene dado por la disminución de la energía superﬁcial σdS. Este trabajo 
es realizado por las moléculas de la superﬁcie sobre el líquido del interior de la gota, y puede expresarse también como 
el trabajo volumétrico realizado sobre el líquido, en función de la variación de volumen dV de la gota, en la forma 
(pgota − pext)dV = pdifdV . Ambas expresiones se reﬁeren al mismo trabajo, en un caso realizado sobre el líquido 
de la gota y en el otro caso realizado por el líquido de la gota, y por tanto son iguales: 


σdS = pdifdV 
(5–2) 


Veamos ahora que ocurre con una gota esférica de radio r, con una superﬁcie S = 4πr2 y un volumen V = 4πr3/3. 
Sustituyendo ésto en la expresión anterior, se tiene: 


σ8πrdr = pdif4πr2dr 


con lo que tenemos que la presión diferencial en el caso de una gota esférica viene dada por: 


pdif = 2σ 


r 
(5–3) 


Un cálculo análogo puede realizarse para una superﬁcie cilíndrica, como la que aparece cuando entre dos planos 
verticales próximos hay agua. En este caso, el volumen es V = πr2L y la superﬁcie S = 2πrL, con L la altura del 
cilindro. Así, para el caso de una superﬁcie cilíndrica, la presión diferencial es: 


pdif = σ 


r 
(5–4) 


Para obtener la relación entre la tensión superﬁcial σ y la fuerza resultante perpendicular a la superﬁcie consideremos 
el elemento diferencial de superﬁcie rectangular que se observa en la ﬁgura 5–3. Si σ es la tensión superﬁcial, r1 
y r2 los radios de curvatura de la superﬁcie en las direcciones indicadas y dl1 y dl2 las longitudes de los lados del 
rectángulo, actuarán unas fuerzas σdl1 y σdl2 en los lados del elemento de superﬁcie. Si dα1 y dα2 son los ángulos 
subtendidos por los lados dl1 y dl2 respectivamente, dα1/2 y dα2/2 serán los ángulos que forman las fuerzas σdl1 y 
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σdl2 con el plano tangente a la superﬁcie, y se cumplirá que dl1 = r1dα1 y dl2 = r2dα2. Por tanto, la componente 
perpendicular de la fuerza de tensión superﬁcial que actúa sobre el lado 1 viene dada por 


2σdl1 sen dα2/2 ≈ 2σdl1dα2/2 = σdl1 
dl2 
r2 
(5–5) 


y la componente perpendicular de la fuerza que actúa sobre el lado 2 por 


2σdl2 sen dα1/2 ≈ 2σdl2dα1/2 = σdl2 
dl1 
r1 
(5–6) 


Para que se esté en condiciones de equilibrio mecánico, la suma de estas fuerzas debe ser compensada por una fuerza 
surgida de la diferencia de presiones entre los dos medios de la forma 


∆F = ∆pdl1dl2 
(5–7) 


Así, en el caso general de una superﬁcie de radios principales de curvatura (mínimo y máximo) r1 y r2, la diferencia 
de presiones viene dada por la ecuación de Laplace (o Laplace-Young). 


pdif = σ 
1 


r1 
+ 1 


r2 


(5–8) 


La diferencia de presiones entre las dos caras de una superﬁcie líquida es proporcional a la suma de las inversas 
de los radios de curvatura en dos direcciones mutuamente ortogonales, siendo la presión mayor siempre hacia la 
parte cóncava de la superﬁcie. Hay que recalcar que la suma 1/r1 + 1/r2 es un invariante respecto de la rotación de 
coordenadas, de forma que la ecuación (5–8) es independiente de la dirección en la que se toman los ejes, obviamente 
de acuerdo con el punto de vista físico. 


La ecuación (5–8) se puede particularizar fácilmente al caso cilíndrico con r1 = ∞ y al caso esférico con r1 = r2 


5.1.2. 
Ángulo de contacto 


Se denominan fuerzas de cohesión a las fuerzas de atracción entre una molécula del líquido y otras moléculas del 
mismo. La fuerza que tiene lugar entre una molécula del líquido y otra sustancia, como la pared de un tubo delgado, 
es una fuerza adhesiva. 


Cuando la fuerza adhesiva es menor que la fuerza de cohesión el ángulo se hace grande, mientras que cuando la 
fuerza de cohesión es menor que la fuerza adhesiva, el ángulo se hace pequeño. 


Estudiemos ahora ésto en función de las tensiones superﬁciales implicadas en el problema. Por el borde de la 
superﬁcie de un líquido que se halla en un recipiente concurren tres medios: (1) sólido, que es la pared del recipiente, 
(2) líquido, y (3) gaseoso. 


En la línea de contacto están aplicadas tres fuerzas de tensión superﬁcial, estando cada una de ellas dirigida según la 
tangente hacia el interior de la superﬁcie de contacto de los dos medios correspondientes. El ángulo entre la superﬁcie 
del líquido y la pared sólida se denomina ángulo de contacto θ. 


La superﬁcie del líquido adquiere una forma tal que la resultante de las tres fuerzas σSL, σSG y σLG no tenga 
componente a lo largo de la pared del recipiente (la reacción de la pared elimina la componente perpendicular a la 
misma), quedando la condición de equilibrio: 


σSG = σSL + σLG cos θ 


donde no se ha incluido la longitud de la línea de contacto, que es la misma en todos los casos. Así, el coseno del 
ángulo de contacto queda: 


cos θ = σSG − σSL 


σLG 
(5–9) 


Como se ve, el ángulo de contacto depende exclusivamente de la naturaleza de los tres medios en contacto (de las 
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Figura 5–4 Equilibrio de fuerzas en la línea de contacto sólido-líquido-gas 
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Figura 5–5 Distintos comportamientos de una gota sobre una superﬁcie. 


tensiones superﬁciales en los límites) y no depende ni de la forma del recipiente ni del valor de g. Estas tensiones 
superﬁciales y las fuerzas asociadas son muy sensibles al estado de las superﬁcies divisorias y a sus purezas. 


Si σSG > σSL, es decir si la tensión superﬁcial sólido-gas es mayor que la tensión superﬁcial sólido-líquido, se 
tiene que cos θ > 0 y que θ < 90◦. El ángulo es agudo y el menisco tiene la forma cóncava de la parte izquierda de 
las ﬁguras 5–4 y 5–5. En este caso se dice que el líquido “moja” a la superﬁcie sólida, una gota del líquido depositada 
sobre la superﬁcie se extiende sobre la misma. 


Cuando σSG < σSL, la tensión superﬁcial sólido-gas es menor que la tensión superﬁcial sólido-líquido, teniéndose 
que en este caso cos θ < 0◦ y que θ > 90◦. Este caso corresponde al menisco convexo, y una gota depositada sobre la 
superﬁcie se recoge sobre sí misma, como se muestra en la parte derecha de las ﬁguras 5–4 y 5–5. En este caso se dice 
que el líquido “no moja” la superﬁcie del sólido. Ejemplos de ángulos de este tipo se tienen en los casos de mercurio 
sobre vidrio (150◦) y de agua sobre paraﬁna (105◦). Las gotas de estos líquidos se elevan intentando disminuir el área 
de contacto con la superﬁcie. 


5.1.3. 
Capilaridad 


La elevación de un líquido por un tubo delgado, llamado tubo capilar, es otro fenómeno de superﬁcie íntimamente 
ligado con la tensión superﬁcial. En este proceso aparecen fuerzas tanto de cohesión como de adhesión. Si sumergimos 
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un tubo delgado de vidrio en agua, hemos visto que, ya que el agua moja bastante bien al vidrio, tendremos un ángulo 
de contacto muy agudo, de forma que la superﬁcie estará muy curvada si el tubo es lo suﬁcientemente delgado. El 
resultado de todo ésto es que se producirá una fuerza ascendente, subiendo el agua por el tubo hasta que el peso de la 
columna de agua iguala esta fuerza. 


r 


h 


⃗F 


θ 


Figura 5–6 Capilaridad en un tubo ﬁno 


Consideremos el caso de un líquido que ha alcanzado una altura h en un tubo capilar de radio r. El tubo está abierto 
en su parte superior y se encuentra, por tanto, a la presión atmosférica. La componente vertical de la fuerza que sostiene 
al líquido es F · cos θ que, al ser la longitud de la línea de contacto l = 2πr, queda como σ2πr cos θ. 


Si despreciamos la ligera curvatura de la superﬁcie del líquido, el volumen del mismo en la columna es πr2h 
quedando, al hacer el equilibrio de las fuerzas: 


σ2πr cos θ = ρghS = ρghπr2 


dando para la altura de la columna capilar: 


h = 2σ cos θ 


ρgr 
(5–10) 
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5.2. 
Fluidos en medios permeables. 


En los capítulos anteriores hemos estudiado el movimiento de ﬂuidos “libres”. Sin embargo, los ﬂuidos también 
pueden moverse a través de ciertos materiales, como ocurre con las aguas subterráneas que se ﬁltran a través de los 
terrenos, o el petróleo que se mueve en las capas petrolíferas a los pozos petrolíferos. Este movimiento es consecuencia 
de la acción conjunta de fuerzas gravitatorias, de presión, y otras de diverso origen, que dan lugar a un desplazamiento 
que se produce siempre hacia las zonas de menor energía potencial. 


Se denominan materiales o medios permeables a aquéllos que permiten el paso de ﬂuidos a través de los mismos, ya 
se trate de gases o de líquidos. Los materiales permeables pueden ser de muy distintos tipos, esponjosos como algunas 
espumas, ﬁbrosos, granulares como la arena de los terrenos, etc... En todos ellos la característica común es que el 
ﬂujo tiene lugar a través de una intrincada red de microcanales, con unas elevadas pérdidas de carga dada la pequeña 
sección de éstos. En el caso de líquidos aparecerán también fuerzas de ascensión capilar junto a las fuerzas que dan 
lugar al ﬂujo. 


Naturaleza de los terrenos como medios permeables. 


Los distintos terrenos consisten en una gran cantidad de partículas pequeñas en contacto, con poros o espacios 
huecos entre las mismas que pueden estar ocupados por agua, gas o vapor. El ﬂuido puede circular de una parte a otra 
del medio a través de multitud de microcanales entre las partículas que comunican los distintos poros, comportándose 
como medios permeables. En el caso de un ﬂuido en reposo en un medio permeable, éste asciende por el medio 


Tabla 5–2 Valores típicos de la subida capilar para una serie de materiales 


Material 
Subida capilar 


Arena gruesa 
2 − 5 cm 
Arena normal 
12 − 35 cm 
Arena ﬁna 
35 − 70 cm 
Arcillas 
2 − 4 m y superior 


permeable por encima de la altura de la superﬁcie libre del líquido de fuera del medio permeable. Esto es debido a 
lo siguiente: El medio permeable está constituido por un complejo conjunto de microcanales, a través de los cuales 
asciende el ﬂuido por capilaridad. En la tabla 5–2 de dan algunos valores típicos de esta subida capilar para una serie 
de materiales. A la hora de modelizar el comportamiento de este medio permeable frente a este efecto, se suele sustituir 
en los cálculos la red real de capilares por un sistema de capilares verticales de un cierto diámetro con una determinada 
densidad de capilares por unidad de superﬁcie. De esta forma se puede tratar el problema, así como el problema que 
surge cuando hay dos medios permeables de distintas características superpuestos. 


La naturaleza del sistema de poros es variada dependiendo del tipo de terreno, aunque hay dos casos límites intere- 
santes, el de las arenas y el de las arcillas. 


Las arenas se componen de partículas macroscópicas de forma redondeada o angulosa. Estas partículas drenan 
fácilmente, no absorben apenas agua, tienen un pequeño potencial capilar y cuando se secan no encogen de tamaño. 
Por el contrario, las arcillas están compuestas de partículas microscópicas de forma laminar. Estas absorben agua con 
facilidad, son bastante impermeables una vez mojadas, tienen un elevado potencial capilar y disminuyen de volumen 
de forma considerable al secarse. 


Las arenas se aproximan más al medio poroso ideal y son representativas del tipo de medio poroso que trataremos. 
Aunque ésto parece imponer serias restricciones al estudio que vamos a hacer, en la mayoría de los problemas reales la 
baja permeabilidad de las arcillas las hace ser relativamente impermeables en comparación con los terrenos de grano 
más grueso. 


Denominaremos nivel freático al lugar geométrico de los puntos que se encuentran a la presión atmosférica. Corres- 
ponde al nivel que se alcanzaría en cada punto de la superﬁcie del terreno si en el mismo se introdujera un piezómetro 
la profundidad suﬁciente. Si tomamos como referencia de medida de presiones la presión atmosférica, la presión del 
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nivel capilar 


nivel freático 


terreno seco 


Figura 5–7 Nivel freático y capilar en un terreno. Por encima del nivel capilar no se encuentran granos de arena húmedos. 


agua es positiva por debajo del nivel freático. Por encima del nivel freático pueden encontrarse granos de arena moja- 
dos debido al efecto de la capilaridad hasta el nivel capilar. La presión del agua capilar es negativa y dará lugar al 
denominado potencial de matriz. 


Porosidad. 


Asumamos que las partículas del medio poroso son todas de forma esférica uniforme. Si denominamos V al volumen 
total y Vh al volumen de los huecos, se deﬁne la porosidad como 


n = Vh 


V 
(5–11) 


que corresponde a la fracción del volumen total ocupada por los poros, y el índice de poros como 


e = 
Vh 


V − Vh 
(5–12) 


que corresponde al cociente entre el espacio ocupado por los poros y el espacio ocupado por las partículas sólidas. 
Ambos están relacionados por las expresiones 


n = 
e 


1 + e 
e = 
n 


1 − n 
(5–13) 


Para un empaquetamiento cúbico de esferas de diámetro d el volumen total (por esfera) es V = d3, mientras que el 
volumen hueco es Vh = d3 − πd3/6 quedando una porosidad de 


n = 1 − π 


6 = 0,476. 


Para un empaquetamiento romboédrico, que representa la forma más compacta de empaquetar una serie de esferas 
uniformes, la porosidad es: 


n = 1 − 
√ 


2 
6 π = 0,26. 


En la tabla 5–3 se dan valores para las porosidades de algunos terrenos naturales. Como hemos visto en los ejemplos 
de empaquetamiento cúbico y romboédrico de esferas uniformes, incluso en un medio poroso ideal el espacio hueco 
no es regular, sino que está compuesto por toda una serie de celdas cavernosas interconectadas por canales estrechos. 
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Tabla 5–3 Porosidades de algunos terrenos naturales. 


Medio 
Porosidad 


Arena uniforme, suelta 
46 
Arena uniforme, densa 
34 
Arcilla glacial blanda 
55 
Arcilla glacial dura 
37 
Arcilla orgánica blanda 
75 
Bentonita blanda 
84 


El medio poroso como un medio continuo. 


Los terrenos naturales contienen partículas que se desvían de forma considerable del comportamiento esférico (como 
hemos visto por ejemplo en las arcillas) y que además están lejos de ser de tamaño uniforme. Por tanto la red de 
microcanales es extremadamente complicada y de hecho es imposible en la práctica tratar de una forma razonable un 
terreno atendiendo de forma detallada a la naturaleza microscópica de su estructura de dichos microcanales. Por esta 
razón no fue posible tratar de forma sistemática el ﬂujo permeable hasta el tratamiento de Darcy. 


Afortunadamente, en los problemas de ﬂujo a través de medios permeables no es necesario tener en cuenta el 
ﬂujo a través de los canales individuales. Nuestro interés principal es el ﬂujo macroscópico, entendido como el ﬂujo 
global a través de una sección que contiene multitud de microcanales, y que puede ser considerada uniforme a efectos 
macroscópicos, en contraste con la distribución aproximadamente parabólica que se tiene para un ﬂuido newtoniano en 
régimen laminar a través de un microcanal individual de sección circular. De este modo, en vez de trabajar con el medio 
permeable real con toda su estructura microscópica, trabajaremos con un medio permeable ﬁcticio de propiedades 
continuas en cada punto que coinciden con el valor medio de las propiedades microscópicas del mismo en el entorno 
de dicho punto. Esto es de nuevo la aproximación del continuo o hipótesis del continuo que se vio para el análisis del 
movimiento de los ﬂuidos, pero ahora extendida al medio permeable a través del cuan ﬂuyen éstos. 


Matemáticamente el ﬂujo permeable es una situación análoga a la que se tiene en los procesos de conducción del 
calor en un cuerpo entre puntos que se encuentran a distintas temperaturas, actuando en ese caso la temperatura como 
magnitud potencial. De esta forma, se produce un ﬂujo de calor a través del medio entre puntos que se encuentran 
a distintas temperaturas, dependiendo este ﬂujo de las características del medio. En el estudio del problema de las 
ﬁltraciones a través de un medio permeable, la magnitud potencial que regula el proceso será la energía por unidad de 
peso del ﬂuido, produciéndose entre puntos en los que el ﬂuido tiene distinta energía por unidad de peso una ﬁltración 
que depende de las características del medio. 


5.2.1. 
Velocidad de descarga y velocidad microscópica promedio. 


La velocidad de descarga o velocidad de ﬂujo se deﬁne como la cantidad de ﬂuido que se ﬁltra a través del medio 
poroso por unidad de área y de tiempo. La velocidad de descarga no corresponde a la velocidad real del ﬂujo del agua 
subterránea ya que en realidad el agua en realidad circula por el medio permeable en una compleja trayectoria no 
rectilínea, recorriendo grandes distancias, sino que tiene un signiﬁcado únicamente macroscópico. 


Obviamente, como el ﬂujo puede tener lugar únicamente entre poros de terreno saturado interconectados entre sí, 
la velocidad a través de cualquier sección puede únicamente considerarse en un sentido estadístico, lo que hacemos 
a través del tratamiento del medio real como un medio continuo. Si denominamos xp a la relación efectiva entre el 
área Ap correspondiente a los poros y el área total A de la sección, está claro que xp = Ap/A. El caudal ﬁltrado es 
entonces 
Q = Sv = xpA¯vm, 
(5–14) 


donde xp¯vm es la velocidad de descarga y a ¯vm se le denomina velocidad microscópica promedio† de ﬁltración. Así, 
esta velocidad es la velocidad promedio a través de los canales individuales, mientras que la velocidad de descarga 


†En algunos textos en castellano se utiliza la denominación velocidad de ﬁltración para la velocidad de descarga. Como ésto no es generalizado 
y en textos sajones es más frecuente lo contrario (denominar velocidad de ﬁltración a la velocidad microscópica promedio), evitaremos en lo posible 
el término velocidad de ﬁltración. 
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es la velocidad promedio a través de todo el medio, incluyendo los granos de terreno. Realmente ni siquiera esta 
velocidad es la velocidad real, sino sólo el promedio de las componentes de la velocidad microscópica en la dirección 
del ﬂujo. Dependiendo de la importancia de las componentes en las direcciones perpendiculares al ﬂujo la velocidad 
microscópica real puede ser aún más alta. 


Indaguemos algo más en xp. Si consideramos una sección A del medio que se encuentra a una altura z y denomina- 
mos Ap a la superﬁcie de esa sección correspondiente a los poros, se tiene que a la altura z: 


xp(z) = Ap(z) 


A 
. 
(5–15) 


Promediando a un cilindro de altura h, el valor promedio de xp será 


¯xp = 1 


h 


h 


0 
xp(z)dz 
(5–16) 


que, dividiendo por A y operando, da 


¯xp = 
1 
hA 


h 


0 
Axp(z)dz = 1 


V 


h 


0 
Ap(z)dz 
(5–17) 


donde V es el volumen total del cilindro, mientras que la integral es el volumen hueco total en el cilindro. Así, el valor 
medio de xp es precisamente la porosidad n. 


En general se trabaja con velocidades de descarga, y todo lo que sigue se referirá a velocidades de descarga. Sólo en 
algunos casos muy concretos en los que se busca el comportamiento microscópico se trabaja con la velocidad real. 


5.2.2. 
La ley de Darcy en una dimensión. 


Como es bien sabido de la mecánica de ﬂuidos libres, para ﬂujo estacionario no viscoso e incompresible, la ecuación 
de Bernoulli 
p 
ρg + ¯v2 


2g + z = Htot = cte 
(3–33) 


establece que la suma de la alturas de presión p 


ρg , de velocidad o cinética ¯v2 


2g y geométrica z es constante en cualquier 


punto a lo largo de una línea de ﬂujo.† En el ﬂujo a través de medios permeables como en el caso del permeámetro que 
se muestra en la ﬁgura 5–8 es necesario tener en cuenta las pérdidas de energía producidas por la resistencia viscosa 
en los poros. Por tanto es necesario considerar la ecuación de Bernoulli generalizada en la forma ya conocida 


pA 
ρg + ¯v2 
A 
2g + zA − ∆h = pB 


ρg + ¯v2 
B 
2g + zB 
(4–1) 


donde ∆h representa la pérdida total de carga (energía por unidad de peso del ﬂuido) en la distancia ∆l. Se denomina 
gradiente hidráulico promedio al cociente 
¯i = −∆h 


∆l . 
(5–18) 


Tomando distancias ∆l cada vez más pequeñas respecto del punto inicial, se tiene, en el límite inﬁnitesimal, el 
gradiente hidráulico en cada punto 


i = − l´ım 
∆l→0 
∆h 
∆l = −dh 


dl 
(5–19) 


En la inmensa mayoría de los problemas de ﬁltraciones, la altura de velocidad es tan pequeña que puede despreciarse. 
Por ejemplo, una velocidad de ﬁltración de 30 cm/seg, que es grande comparada con las velocidades de ﬁltración 
típicas, daría una altura de velocidad de 0,5 cm. Por tanto, la ecuación de Bernoulli para ﬂujo en medios permeables 
puede escribirse a efectos prácticos en la forma 


pA 
ρg + zA − ∆h = pB 


ρg + zB 
(5–20) 


†Si el ﬂujo es no rotacional, ésto es cierto en general entre dos puntos cualesquiera del tubo de ﬂujo. 
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Figura 5–8 Flujo a través de un permeámetro. 


a través de las alturas piezométricas, que es lo que realmente podemos medir. 


Por tanto, en los problemas normales de permeación, es indiferente hablar de altura piezométrica o total y se deno- 
mina a ésta potencial hidráulico Φ que tiene, en un punto genérico del ﬂujo permeable, la forma 


Φ = Htot = p 


ρg + z 
(5–21) 


Debido a la compleja naturaleza de los terrenos, antes de 1856 el estudio del ﬂujo permeable del agua en los mismos 
desaﬁaba a todos los análisis. En ese año, Henry Darcy publicó una relación simple basada en sus experimentos en 
“Les fontaines publiques de la ville de Dijon”[1]. Para ello utilizó el permeámetro que se muestra en la ﬁgura 5–9 y, 
para distintos valores del caudal, analizó la pérdida de carga entre el principio y el ﬁnal del permeámetro. De esta forma 
observó que para caudales no demasiado grandes se satisface una relación de proporcionalidad entre la velocidad de 
descarga v = Q/S y el gradiente hidráulico o pérdida de carga por unidad de longitud del medio permeable. Esta 
ecuación, 
Q 
S = v = k ∆h 


∆l 
⇒ 
v = ki = −k dΦ 


dl 
(5–22) 


conocida como ley de Darcy, muestra una dependencia lineal entre el gradiente hidráulico y la velocidad de descarga v. 
En Mecánica del Suelo, se denomina al coeﬁciente de proporcionalidadk coeﬁciente de permeabilidad o conductividad 
hidráulica que como se ve claramente en la ecuación (5–22) tiene dimensiones de velocidad. Este coeﬁciente representa 
la facilidad que ofrece el medio para dejarse atravesar por el ﬂuido. Su valor depende del tamaño promedio de los poros 
y está relacionado también con la distribución de tamaños y formas de las partículas y con la estructura del medio. Así, 
éste puede cambiar en un factor de hasta 106 entre el de una grava y el de una arcilla. Incluso una pequeña proporción 
de material ﬁno en uno más grueso puede cambiar de forma importante su magnitud. 


Aunque en la ecuación (5–22) la ley de Darcy se presenta en forma diferencial es necesario recalcar que en modo 
alguno describe la situación microscópica en la red de poros del medio. De forma estricta, la ley de Darcy representa 
el equivalente estadístico macroscópico de las ecuaciones del movimiento de Navier-Stokes para el ﬂujo viscoso a 
través de un medio permeable. Es precisamente esta equivalencia la que permite el tratamiento del ﬂujo permeable en 
el marco de una teoría de campos potenciales. A nuestros efectos, todas las contribuciones viscosas serán consideradas 
a través de la ley de Darcy, y el ﬂujo a todos los demás efectos será considerado como no viscoso o perfecto. 
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Figura 5–9 Permeámetro original utilizado por Darcy. Extraído de su informe Les fontaines publiques de la ville de Dijon, 1856[1]. 


Rango de validez de la ley de Darcy. 


Como se ha visto en la sección 4.1.2 la observación visual de colorantes inyectados en los líquidos llevó a Reynolds 
(1883) a concluir que el orden dentro del ﬂujo depende de la velocidad del mismo. A velocidades pequeñas el ﬂujo era 
ordenado, en capas, es decir, laminar. Por el contrario, al aumentar la velocidad se produce una mezcla entre el agua y 
el colorante, el ﬂujo se hace irregular o turbulento. 


Las distintas experiencias realizadas muestran de forma concluyente que para pequeñas velocidades (en el rango 
de ﬂujo laminar), la ley de Darcy (5–22) representa de forma adecuada el ﬂujo a través de un medio permeable, con 
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una relación lineal entre el gradiente hidráulico y la velocidad del ﬂujo. Según el ﬂujo se va haciendo turbulento, el 
comportamiento del gradiente hidráulico se aproxima mejor mediante una expresión de la forma 


i = −dΦ 


dx = av + bvn 
(5–23) 


donde a y b son constantes y n varía entre 1 y 2, dependiendo de las características del ﬂujo. Por ejemplo, para el ﬂujo 
turbulento de agua a través de columnas de perdigones de tamaño uniforme, se obtiene n = 2 (Lindquist). 


Al igual que en el caso de los ﬂuidos libres (véase la sección 4.5), el criterio para caracterizar el tipo de régimen, 
laminar o turbulento, es utilizar el número de Reynolds (número adimensional que relaciona las fuerzas de inercia con 
las viscosas), pero adaptado al ﬂujo a través de medios permeables, 


Re = vDρ 


µ 
= vD 


ν , 
(4–23) 


donde ahora v es la velocidad de descarga, D el diámetro promedio de las partículas, ρ la densidad del ﬂuido, µ su 
viscosidad dinámica y ν su viscosidad cinemática. 


El valor crítico del número de Reynolds en el cual el ﬂujo permeable cambia de laminar a turbulento está entre 1 y 
12. Sin embargo para nuestros efectos, será suﬁciente dar por válida la ley de Darcy cuando el número de Reynolds es 
menor o igual que la unidad. 
vDρ 


µ 
= vD 


ν 
≤ 1 
(5–24) 


Sustituyendo los valores de ρ y de µ para el agua en la ecuación anterior y asumiendo un valor de la velocidad de 
descarga de 0,25 cm/seg, se obtiene D = 0,4 mm. Este es un valor típico del tamaño medio de partículas para la arena 
gruesa. Por tanto parece bastante razonable que para medios permeables normales, en los cuales el tamaño medio de 
grano es considerablemente menor, la validez de la ley de Darcy sea clara. De hecho, hay que subrayar que el carácter 
laminar del ﬂujo del agua a través de medios permeables en las condiciones dadas es uno de los pocos ejemplos válidos 
de este tipo de ﬂujo en situaciones reales de ingeniería hidráulica. 


La conductividad hidráulica. 


Hemos visto que para ﬂujo laminar a través de un medio permeable se satisface la Ley de Darcy 


v = −k dΦ 


dl 


donde k es la conductividad hidráulica. La densidad y la viscosidad del ﬂuido ejercen una inﬂuencia en la velocidad 
con la que éste circula a través del medio permeable, por lo que es interesante aislar la parte de k que depende de 
esas propiedades. Para ello, introducimos la permeabilidad física o permeabilidad intrínseca k0 (con unidades de 
superﬁcie), que es un parámetro que caracteriza las propiedades del medio y es independiente de las propiedades del 
ﬂuido. La relación entre la permeabilidad física y la conductividad hidráulica viene dada por 


k = k0 
ρg 
µ = g 


ν k0 
(5–25) 


Sustituyendo esta expresión en la ecuación de Darcy (5–22), se tiene 


v = −k0 
ρg 
µ 
dΦ 
dl = −k0 
g 
ν 
dΦ 
dl 
(5–26) 


que indica que la velocidad de descarga es inversamente proporcional a la viscosidad del ﬂuido. Es interesante utilizar 
esta última expresión cuando se trabaja con más de un ﬂuido o cuando las propiedades de éste varían por ejemplo 
a causa de cambios de temperatura. Aunque la permeabilidad física es idealmente la misma para distintos ﬂuidos y 
depende únicamente del medio permeable, en la práctica éste puede incluso cambiar de dimensiones con distintos 
ﬂuidos obteniéndose en ese caso distintos valores de la permeabilidad física 
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Tabla 5–4 Valores aproximados de la conductividad hidráulica y el tamaño medio de grano para distintos medios porosos. 


Tipo de terreno 
Conductividad hidráulica 
k (cm/seg) 
Tamaño medio 
de grano (mm) 


Gravas 
≥ 1, 0 
≥ 2 
Arena gruesa 
1,0 − 10−1 
0,6 − 2 
Arena media 
10−1 − 10−2 
0,2 − 0,6 
Arena ﬁna 
10−2 − 10−3 
0,06 − 0,2 
Arena limosa 
10−3 − 10−4 
0,02 − 0,06 
Limo 
10−4 − 10−5 
0,006 − 0,02 
Turba o limo ﬁno 
10−5 − 10−6 
0,0002 − 0,006 
Arcilla meteorizada 
10−6 − 10−7 


Arcilla no meteorizada 
≤ 10−7 


En los problemas típicos de ﬂujo a través de medios permeables, en los cuales estaremos interesados en el ﬂujo de 
un único ﬂuido relativamente incompresible es más sencillo utilizar la ley de Darcy con k. Algunos valores típicos de 
la conductividad hidráulica para diversos medios se dan en la tabla 5–4. 


Puede relacionarse la permeabilidad física k0 con la porosidad n del material para ﬂujo en régimen laminar como 
el que se tiene en la región de validez de la ley de Darcy a partir de la fórmula de Poiseuille (4–21) y de la relación 
(5–14). En la drástica aproximación de tuberías cilíndricas de radio R el correspondiente al tamaño promedio de los 
poros, la velocidad de descarga queda en la forma 


¯vm = nQ 


A = nQ 


8 
ρg 
µ 
∆hpiezo 


L 


que al comparar con la ley de Darcy en función de la permeabilidad física (5–26) da un valor para la permeabilidad 
física en la aproximación que utilizamos de 


k0 = nR2 


8 
(5–27) 


5.2.3. 
Potencial hidráulico y sus componentes. El campo hidráulico. 


Como hemos dicho anteriormente, la resultante de las fuerzas que actúan en cada punto del ﬂujo sobre el ﬂuido 
puede obtenerse como el gradiente de la energía potencial en dicho punto del medio poroso. Cuando se realiza un 
trabajo sobre el ﬂuido que circula por el medio poroso, esta energía potencial cambia en una cantidad igual al trabajo 
realizado. Se denomina potencial hidráulico Φ a esta variación de la energía potencial del ﬂuido respecto de la que 
tenía en un estado de referencia, medida por unidad de peso del ﬂuido, es decir, la energía total por unidad de peso del 
ﬂuido. Como hemos visto, hay distintas contribuciones a este potencial hidráulico, el potencial de presión o altura de 
presión Φp = p/(ρg), el potencial gravitacional o altura geométrica Φz = z, y el potencial de velocidad o altura de 
velocidad o cinética Φv = v2/(2g), que como hemos visto es despreciable en los problemas habituales de movimiento 
de ﬂuidos a través de un medio permeable. 


Hay otra contribución al potencial hidráulico de la que aún no hemos hablado de forma explícita, y es el potencial 
de matriz o de succión, Φm. Esta contribución está ligada a la presencia de fuerzas de ascensión (o descenso) capilar. 
Éstas aparecen como consecuencia de la interacción del agua con el esqueleto del medio poroso o matriz de partículas 
sólidas en que está embebida (véase la sección 5.1.3). Estas fuerzas capilares son en general atractivas para los suelos 
y cerámicas porosas, pero pueden ser negativas para otros materiales porosos. Cuando las fuerzas son atractivas el 
potencial de matriz tiene valores negativos. 


Finalmente mencionar el potencial osmótico, que únicamente aparece cuando existen membranas semipermeables 
en el medio poroso (por ejemplo, las raíces de las plantas que se encuentran en el suelo) y que no tendremos en cuenta. 


Cuando entre dos puntos existe una diferencia de potencial hidráulico aparece una fuerza que intenta producir 
movimiento del ﬂuido hacia posiciones de menor energía potencial, y este movimiento se producirá únicamente entre 
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puntos con distintos potenciales hidráulicos. En general, el potencial hidráulico puede variar con el tiempo, y además 
no está distribuido de forma uniforme (en ese caso no habría movimiento). Por tanto, lo que en realidad se tiene es una 
distribución espacio-temporal del potencial hidráulico, el denominado campo hidráulico 


Φ = Φ(x, y, z, t) 
(5–28) 


Se denomina superﬁcie equipotencial hidráulica al lugar geométrico de los puntos del medio permeable en los cuales 
el potencial hidráulico tiene el mismo valor. Obviamente, como el potencial hidráulico no puede tener valores distintos 
en el mismo punto en el mismo instante, las superﬁcies equipotenciales no pueden cruzarse. 


5.2.4. 
La ley de Darcy en tres dimensiones. 


Si en un punto genérico del ﬂujo permeable a través de un medio isótropo observáramos la experiencia anterior en 
tres direcciones mutuamente perpendiculares, x, y y z, se tendrían las expresiones para la ley de Darcy (5–22) en cada 
una de estas direcciones como 


vx = kix = −k dΦ 


dx 


vy = kiy = −k dΦ 


dy 


vz = kiz = −k dΦ 


dz 


que puede escribirse como 


⃗v = −k 
dΦ 


dx ⃗ux + dΦ 


dy ⃗uy + dΦ 


dz ⃗uz 


= −k⃗∇Φ. 
(5–29) 


para un medio permeable isótropo. 


Si el medio es anisótropo la situación es más compleja, y en general la conductividad hidráulica viene representada 
por un tensor con varias componentes distintas de cero, expresándose la ley de Darcy en un punto en la forma 


⃗v = 


 


 
vx 
vy 
vz 


 


 = − 


 


 
kxx kxy kxz 
kyx kyy kyz 
kzx kzy kzz 


 


 


 


 
 
 
 
 
 
 


dΦ 
dx 
dΦ 
dy 
dΦ 
dz 


 


 
 
 
 
 
 
 


= −ˆk⃗∇Φ 
(5–30) 


El tensor ˆk es un tensor simétrico de segundo orden, y por lo tanto existe un sistema de referencia (sistema de ejes 
principales) en el cual todos los elementos de fuera de la diagonal son nulos, denominándose a los elementos de 
la diagonal en ese sistema elementos principales del tensor o, en este caso concreto, conductividades hidráulicas 
principalesconductividad hidráulica!principal. En este sistema de ejes principales, r1, r2 y r3, la expresión anterior 
toma la forma 


⃗v = 


 


 
v1 
v2 
v3 


 


 = − 


 


 
k1 0 0 
0 k2 0 
0 0 k3 


 


 


 


 
 
 
 
 
 
 


dΦ 
dr1 
dΦ 
dr2 
dΦ 
dr3 


 


 
 
 
 
 
 
 


. 
(5–31) 


El vector gradiente es un vector que es siempre perpendicular a las superﬁcies equipotenciales del campo escalar en 
cada punto. Por tanto al ser, para un material isótropo ⃗v = −k⃗∇Φ la velocidad de descarga es siempre perpendi- 
cular a las superﬁcies de potencial hidráulico constante para un material isótropo. Por el contrario, para un material 
anisótropo, aunque el gradiente sigue siendo perpendicular a las equipotenciales, la velocidad no lo será al aparecer 
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la conductividad hidráulica en forma de tensor en la ecuación (5–31) con distintos valores de k. Sin embargo, sí se 
cumplirá para los materiales isótropos ﬁcticios equivalentes a uno anisótropo que se obtienen redimensionando los 
ejes y se tratan en la sección 5.2.7. 


5.2.5. 
Flujo a través de terrenos saturados homogéneos e isótropos. 


La conductividad hidráulica de un suelo saturado es constante si su estructura es estable, porque todos sus poros 
y capilares están ocupados por el agua, igual que cualquier otro material poroso en idénticas condiciones. Por el 
contrario, la conductividad hidráulica de un suelo no saturado puede cambiar continuamente en función de los cambios 
que se produzcan en el potencial de matriz. Asimismo, en algunas situaciones como durante las lluvias se tienen 
cambios temporales del contenido de agua en el suelo. Lo mismo ocurre si se está produciendo una degradación del 
medio o del ﬂuido durante el proceso de ﬂujo permeable. 


En el caso general en el que puede haber variación temporal del ﬂujo, es necesario tener en cuenta no sólo la ley de 
Darcy, sino también la ecuación de continuidad 


∂ρ 
∂t = −⃗∇(ρ⃗v). 
(3–9) 


Sin embargo cuando se tiene en cuenta que la compresibilidad de los líquidos es muy pequeña, la densidad tiene una 
variación despreciable con el tiempo y la posición, y la expresión anterior se puede escribir como 


⃗∇⃗v = 0 
(3–11) 


Estudiaremos ahora como es la ecuación diferencial que describe el ﬂujo. Para ello nos limitaremos al caso en que el 
medio permeable es estable, a régimen estacionario y a un ﬂuido incompresible como el agua. 


Si el medio es además isótropo, el tensor conductividad hidráulica es un tensor diagonal con todas las componentes 
de la diagonal iguales a un único valor de la conductividad hidráulica, y la ley de Darcy puede escribirse como 


⃗v = −k⃗∇Φ 
(5–29) 


Sustituyendo esta expresión en la ecuación anterior se tiene que la ecuación diferencial que hay que resolver para 
obtener la distribución de potencial hidráulico es 


∇2Φ = ∂2Φ 


∂x2 + ∂2Φ 


∂y2 + ∂2Φ 


∂z2 = 0. 
(5–32) 


Esta ecuación diferencial es la conocida ecuación de Laplace, ampliamente estudiada en el electromagnetismo y en 
el tratamiento de algunos problemas de conducción del calor. Como se observa en esta ecuación, la forma del campo 
hidráulico no depende de la conductividad hidráulica, sino que depende únicamente de las condiciones de contorno en 
los límites, es decir de la forma del recinto y de las características de estos límites. Obviamente, el ﬂujo sí depende de 
la conductividad hidráulica a través de la ley de Darcy. 


5.2.6. 
La ecuación de Laplace y su resolución. 


Como se decía anteriormente, la ecuación de Laplace ha sido ampliamente estudiada en Física Matemática por 
su interés en diversos campos que van desde electromagnetismo a la conducción del calor. Su resolución depende 
principalmente de la geometría del recinto espacial en el que se veriﬁca ésta (como por ejemplo un suelo, un material 
poroso como la zapata de cimentación de un ediﬁcio, un muro de carga o una presa de gravedad) y de las condiciones 
de contorno que se cumplen en los límites del recinto. Es decir, la resolución de la ecuación de Laplace depende de 
las características o los valores que ha de tener el potencial hidráulico (o su gradiente, directamente relacionado con la 
velocidad de descarga) en la superﬁcie límite del recinto en estudio. Estas condiciones de contorno vienen, en general, 
impuestas por la naturaleza física de estos límites, como por ejemplo, la presencia de una capa impermeable. 


Se dice que estas condiciones de contorno son condiciones de Dirichlet o de primera especie cuando lo que se 
especiﬁca son los valores del potencial hidráulico en la superﬁcie límite, mientras que se dice que estas condiciones 
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son condiciones de von Neumann o de segunda especie si lo que se especiﬁca son condiciones en el gradiente del 
potencial hidráulico (en la velocidad). Se dice que se tienen unas condiciones de contorno mixtas cuando sobre una 
parte de la superﬁcie límite se conocen los valores del potencial hidráulico y sobre otra parte se conoce información 
sobre la velocidad. 


Las técnicas de resolución de la ecuación de Laplace son muy conocidas en el ámbito de la Física Matemática y 
utilizan fundamentalmente el método de separación de variables en coordenadas cartesianas, esféricas o cilíndricas, 
según la simetría del recinto permeable. Estos son los casos más favorables y al mismo tiempo, los menos frecuen- 
tes. Si el recinto no tiene una simetría apropiada se utiliza, como técnica alternativa, la transformación conforme en 
variable compleja o el método de imágenes. Para la mayor parte de los problemas que son intratables analíticamente, 
es necesario recurrir a métodos numéricos de resolución de ecuaciones diferenciales, como los de diferencias ﬁni- 
tas o elementos ﬁnitos. A menudo se utilizan estos métodos también para abordar problemas que, aunque resolubles 
analíticamente, requieren un excesivo esfuerzo. 


Condiciones de contorno para la resolución de la ecuación de Laplace. 


En un medio permeable saturado aparecen distintos tipos de condiciones de contorno, todas presentes en el ejemplo 
de la presa de tierra que se muestra en la ﬁgura. Consideremos en detalle las características de cada una de ellas. 


A 


B 


C 


D 


h 


Superﬁcie ﬁltrante 


Superﬁcie impermeable 


Superﬁcie libre 


Superﬁcie de escorrentía 


Figura 5–10 Presa de tierra en la que aparecen los cuatro tipos de condiciones de contorno que estudiaremos. 


Superﬁcies impermeables. 


Una superﬁcie impermeable separa el medio permeable (saturado) de un medio impermeable con el que está en 
contacto. El ﬂuido no puede penetrar a través de una superﬁcie impermeable, ya que las características del material 
impermeable lo impiden. Por tanto, la componente de la velocidad perpendicular a la superﬁcie debe ser cero en cada 
punto. Sin embargo no se impone restricción alguna a sus componentes tangenciales. Si denominamos n y t a las 
direcciones normal y tangencial, respectivamente, en un punto de la superﬁcie, se tiene 


∂Φ 
∂n = 0. 
(5–33) 


es decir, una superﬁcie impermeable deﬁne el lugar geométrico de una serie de líneas de ﬂujo. Al mismo tiempo una 
línea de ﬂujo satisface la condición de una superﬁcie impermeable y, de hecho, las paredes de un tubo de ﬂujo, se 
consideran a todos los efectos como superﬁcies impermeables. 


Hay dos tipos de superﬁcies impermeables típicas. Una es la superﬁcie impermeable que tiene el ﬂuido por encima, 
como ocurre cuando se considera el ﬂujo a través de un medio permeable que se encuentra sobre una base de otro 
material de conductividad hidráulica mucho menor que la del medio permeable. En este caso la superﬁcie impermeable 
está deﬁnida por las líneas de corriente más bajas. También puede tenerse una estructura impermeable sobre el terreno, 
como una presa de gravedad, o la pared vertical de un muro. 
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En un material isótropo, al ser las superﬁcies equipotenciales del campo hidráulico perpendiculares a la velocidad 
de descarga, éstas deben ser perpendiculares asimismo a la superﬁcie impermeable. Ésto no es cierto para un mate- 
rial anisótropo en el que formarán ángulos distintos de 90◦. Sin embargo sí se cumple para los materiales isótropos 
equivalentes a un anisótropo que se tratarán en la sección 5.2.7. 


Superﬁcies ﬁltrantes. 


Las superﬁcies ﬁltrantes corresponden a los límites del medio permeable en contacto con los límites del ﬂuido libre. 
En estas superﬁcies ﬁltrantes el valor del potencial hidráulico será el del ﬂuido libre en contacto con la misma, es 
decir, 
Φ = p 


ρg + z = cte. 
(5–34) 


Las superﬁcies ﬁltrantes en contacto con una masa de ﬂuido libre son por tanto superﬁcies equipotenciales, cuyo 
potencial hidráulico es igual a la altura de carga total del ﬂuido en contacto con el medio permeable. Obviamente, 
distintas superﬁcies ﬁltrantes del medio permeable en contacto con distintas masas de ﬂuido libre con distintas alturas 
de carga totales tendrán distintos potenciales hidráulicos. 


Para materiales isótropos la velocidad de descarga será perpendicular a las superﬁcies equipotenciales del campo 
hidráulico y por tanto a las superﬁcies ﬁltrantes. Para materiales anisótropos esto no será cierto, a no ser que se trate 
de los materiales isótropos equivalentes a uno anisótropo que se tratan en la sección 5.2.7. 


Superﬁcie de escorrentía. 


La superﬁcie de escorrentía representa un límite en el cual el ﬂujo emergente del medio permeable sale a una zona 
libre de líquido y de terreno, es decir, sale a la atmósfera. Como la presión a la salida es la atmosférica, se cumple en 
esta superﬁcie que el potencial hidráulico es 
Φ = z 
(5–35) 


En principio no hay restricciones en ⃗v, salvo que obviamente debe tener una componente perpendicular a la superﬁcie. 


Superﬁcie libre o línea de ﬁltración. 


La superﬁcie libre está compuesta por las líneas de ﬂujo superiores en el terreno permeable y coincide con el nivel 
freático. Esta superﬁcie se encuentra dentro del medio permeable y separa la región de ﬂujo saturado de la región del 
medio permeable a través de la cual no hay ﬂujo. Se diferencia de la superﬁcie impermeable en que en este caso son 
las características del ﬂujo las que hacen que no se atraviese esta superﬁcie, en vez de las características del material 
de la superﬁcie, que para la superﬁcie libre es el mismo que en la zona en la que sí hay ﬂujo. Su determinación es uno 
de los principales problemas que surgen a la hora de estudiar de forma analítica el problema del ﬂujo permeable. Por 
un lado esta superﬁcie debe comportarse como una superﬁcie impermeable, al estar compuesta por líneas de ﬂujo, y 
por otro lado es evidente que la presión en cada punto de esa superﬁcie es constante e igual a la presión atmosférica, 
por lo que el potencial hidráulico será 
Φ = z 
(5–36) 


lo que muestra que el potencial hidráulico a lo largo de una superﬁcie libre varía linealmente con la altura. 


5.2.7. 
Flujo permeable a través de medios anisótropos 


Hemos considerado el caso de un ﬂujo isótropo en el que el coeﬁciente de permeabilidad del medio en un punto es 
independiente de la dirección del ﬂujo en ese punto. Si además el medio tiene el mismo coeﬁciente de permeabilidad 
en todos los puntos del recinto se dice que es homogéneo e isótropo. Cuando el coeﬁciente de permeabilidad depende 
de la dirección y la dependencia es la misma para todos los puntos se dice que el medio es homogéneo y anisótropo. 


La mayor parte de los terrenos son anisótropos en alguna medida. Por ejemplo, los terrenos sedimentarios están 
habitualmente formados por ﬁnas capas alternas y además, la tensión de compresión en la dirección horizontal es 
distinta de la que se tiene en la dirección vertical. Generalmente en terrenos sedimentarios homogéneos naturales, el 
coeﬁciente de permeabilidad es mayor en la dirección horizontal que en la vertical. En esta sección se abordará lo 
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x 


y 


b 


L 
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k2 


kn−2 


kn−1 


kn 


d1 


d2 


dn−1 


dn−2 


dn 


... 


Paralelo 


Perpendicular 


Figura 5–11 Flujo paralelo y perpendicular a través de un material permeable formado por capas de distintos materiales. 


siguiente: 


Flujo paralelo y perpendicular a través de un medio anisótropo compuesto de capas isótropas. 


Consideremos una sección vertical de un terreno estratiﬁcado formado por n capas isótropas, cada una de ellas de 
espesor d1, d2, ..., dn, con unas conductividades hidráulicas k1, k2, ..., kn, respectivamente y analicemos dos casos 
particulares interesantes. 


Flujo paralelo: 
Consideremos en primer lugar el ﬂujo paralelo a las capas. El caudal ﬁltrado será la suma de los 
caudales ﬁltrados a través de cada una de las capas, 


Qi = dibki 
hi − hf 


L 
⇒ 
Q = 


n 


i=1 
Qi = bhi − hf 


L 


n 


i=1 
kidi 


donde h1 − h2 es la pérdida de carga en las capas en una distancia L. La velocidad promedio en la dirección x será 


vx = Q 


S = 


n 


i=1 


kidi 


d 
hi − hf 


L 


por lo que se puede escribir una conductividad hidráulica efectiva cuya componente x es 


kx = 


n 


i=1 


kidi 


d 
(5–37) 


Flujo perpendicular: 
Consideremos ahora el ﬂujo perpendicular a las capas. Para que se satisfaga la ecuación de 
continuidad la componente y de la velocidad debe ser la misma en todas las capas. Como la pérdida total de carga 
debe ser igual a la suma de las pérdidas de carga en cada capa, se tiene 


∆hi = di 


ki 
v 
⇒ 
∆h = 


n 


i=1 
∆hi = v 


n 


i=1 


di 
ki 


y 


v = hi − hf 


n 


i=1 


di 
ki 


= 
d 


n 


i=1 


di 
ki 


hi − hf 


d 
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que corresponde a una conductividad hidráulica efectiva de componente y 


ky = 
d 


n 


i=1 
di/ki 


(5–38) 


Flujo a través de un medio anisótropo general. 


Para estudiar lo que ocurre en un medio anisótropo consideraremos el problema descrito en los ejes principales del 
tensor conductividad hidráulica. En esos ejes, en función de las conductividades hidráulicas principales k1, k2 y k3, 
las componentes de la velocidad de descarga que se obtienen de la Ley de Darcy quedan: 


v1 = −k1 
∂Φ 
∂r1 
(5–39a) 


v2 = −k2 
∂Φ 
∂r2 
(5–39b) 


v3 = −k3 
∂Φ 
∂r3 
(5–39c) 


En el sistema de ejes principales la ecuación de continuidad se escribe como 


⃗∇⃗v = ∂v1 


∂r1 
+ ∂v2 


∂r2 
+ ∂v3 


∂r3 
= 0 
(5–40) 


En un medio homogéneo como el que estamos considerando, [kij] no depende del punto, y por tanto es constante 
frente a las derivadas espaciales, quedando, al sustituir la ley de Darcy en la ecuación de continuidad 


k1 
∂2Φ 
∂r2 
1 
+ k2 
∂2Φ 
∂r2 
2 
+ k3 
∂2Φ 
∂r2 
3 
= 0. 
(5–41) 


Si consideramos el siguiente cambio de variable, en función de una constante arbitraria k con dimensiones de conduc- 
tividad hidráulica, 


q1 = 
k 


k1 
r1 


q2 = 
k 


k2 
r2 


q3 = 
k 


k3 
r3 


⇐⇒ 


r1 = 



k1 
k q1 


r2 = 



k2 
k q2 


r3 = 



k3 
k q3 


(5–42) 


la ecuación (5–41) se puede reescribir como 


k ∂2Φ 


∂q2 
1 
+ k ∂2Φ 


∂q2 
2 
+ k ∂2Φ 


∂q2 
3 
= 0. 
(5–43) 


que, eliminando k se puede escribir como 


∂2Φ 
∂q2 
1 
+ ∂2Φ 


∂q2 
2 
+ ∂2Φ 


∂q2 
3 
= ∇2 
qΦ = 0. 


que es de nuevo la ecuación de Laplace, pero ahora en el sistema de coordenadas reescalado. Así, mediante una simple 
expansión o contracción de las coordenadas espaciales podemos transformar una región del ﬂujo anisotrópa en una 
región isótropa ﬁcticia en la que la ecuación de Laplace es válida y en la que en principio podríamos obtener una 
solución para ésta. Una vez se haya resuelto el problema para esta región ﬁcticia, para obtener la solución para el ﬂujo 
real que se tiene basta con invertir la transformación. 


El hecho de que se cumpla la ley de Darcy (5–31) no puede verse afectado por la elección del sistema de ejes y 
en concreto por el redimensionado de los mismos que acabamos de tratar. No solo eso, si en la ley de Darcy (5–31) 
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realizamos el cambio de variables (5–42) tanto en las componentes dri/dt de la velocidad como en las ∂Φ/∂ri del 
gradiente la ley de Darcy queda, para el material anisótropo en el sistema de coordenadas reescalado en la forma 
⃗vq = −k⃗∇qΦ similar a la ley de Darcy para un material isótropo (5–29) en función del parámetro k. Esto implica 
que en el material anisótropo analizado en estos ejes como un material isótropo ﬁcticio la velocidad de descarga es 
también perpendicular a las superﬁcies equipotenciales del campo hidráulico expresado también en estos ejes. 


5.2.8. 
La red de ﬂujo. 


Se denomina red de ﬂujo a una representación gráﬁca de las líneas de corriente y de sus correspondientes superﬁ- 
cies equipotenciales en una región del ﬂujo. La construcción de la red de ﬂujo es uno de los métodos más valiosos 
empleados en el análisis de los problemas de ﬂujo en dos dimensiones, aunque a menudo requiere tediosos ajustes de 
prueba y error. Al construir la red de ﬂujo se llena la zona en permeación con una serie de líneas de ﬂujo y otra serie 
de superﬁcies equipotenciales que forman un ángulo recto entre unas y otras en los puntos en que se cruzan. 


Nótese que la construcción de la red de ﬂujo sólo proporciona una solución aproximada al problema, aunque puede 
ser buena incluso en el caso de geometrías complejas a partir de una serie de reglas sencillas (inicialmente debidas 
a Phillip Forchheimer en torno a 1900 y después mejoradas por Arthur Casagrande en 1937). El método habitual de 
construcción de la red de ﬂujo es el siguiente: 


1. Dibujar a escala los límites del recinto de forma que puedan hacerse terminar en dichos límites todas las líneas de 
corriente y equipotenciales. En caso de tratarse de un material anisótropo homogéneo debe primero reescalarse 
el problema en la forma descrita en la sección anterior. 


2. Trazar de forma aproximada unas pocas líneas de ﬂujo, teniendo en cuenta las condiciones en los límites del 
recinto. Es decir, entre otras cosas, que las líneas de corriente han de ser tangentes a las superﬁcies impermeables 
y perpendiculares a las superﬁcies ﬁltrantes. No debe perderse la noción de que éstas son sólo unas pocas de 
entre el continuo de líneas de corriente. Como una ayuda a la hora de espaciar las líneas de corriente, obsérvese 
que la distancia entre líneas de corriente consecutivas aumenta en la dirección de mayores radios de curvatura. 
Es conveniente empezar en las zonas más sencillas, por ejemplo en aquéllas en las que las líneas de ﬂujo 
son paralelas entre sí, lo mismo que las equipotenciales, de forma que los “cuadrados” de la red de ﬂujo son 
realmente cuadrados, y después seguir hacia las zonas de mayor complejidad. 


3. Trazar las equipotenciales correspondientes, considerando que deben intersectar a ángulos rectos a todas las 
líneas de corriente para materiales isótropos o anisótropos reescalados. Al inscribir un círculo en uno de estos 
“cuadrados curvilíneos” éste debe tocar cada lado sólo una vez. Estos “cuadrados” deben cambiar de tamaño 
gradualmente. 


4. Aﬁnar las líneas de corriente y equipotenciales trazadas para que se satisfagan bien las condiciones anteriores. 
Asimismo debe tenerse en cuenta la simetría del problema por si es posible simpliﬁcar la representación. 


5. Se recomienda como una prueba adicional comprobar que las diagonales de los “cuadrados” formados por líneas 
de corriente y equipotenciales forman curvas sin variaciones excesivamente bruscas que se cruzan formando 
ángulos rectos. 


5.2.9. 
Ejemplos de redes de ﬂujo en distintos sistemas 


Presas de tierra 


Cuando las presas de tierra están en las condiciones que se muestran en la ﬁgura 5–10, con una superﬁcie de 
escorrentía de esa forma, puede producirse perforación de la presa como consecuencia de una erosión interna que 
poco a poco va arrastrando la tierra y produciendo un canal por el interior de la presa. Esta erosión, denominada 
tubiﬁcación, que comienza en la parte aguas abajo y va propagándose por el interior de la presa a través de líneas de 
conductividad hidráulica máxima puede evidenciarse muchos años después de que ésta se ha construido. 


Para minimizar el riesgo de perforación se utilizan distintos sistemas, como paredes impermeables enterradas o la 
colocación de un núcleo impermeable en la presa. Uno de los sistemas consiste en controlar el efecto de la ﬁltración 
mediante un dren interno como se muestra en la ﬁgura 5–12 diseñado de forma que la superﬁcie libre está siempre por 


Agustin Martin Domingo 


5.2. Fluidos en medios permeables. 
91 


h 


base impermeable 


superﬁcie ﬁltrante 


superﬁcie libre 


superﬁcie de escorrentía 


dren 


Figura 5–12 Red de ﬂujo en una presa de tierra típica con un dren interno. 
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Φ0 ≃ 0 m 
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Φ = 0,9 h 


Φ10 = h 


Figura 5–13 Red de ﬂujo en un presa de gravedad con su base al nivel del suelo.Las superﬁcies equipotenciales se muestran con líneas de trazos y 
las líneas de ﬂujo con líneas continuas. 


h 


Φ0 ≃ 0 m 


Φ1 
Φ9 


Φ10 = h 


Figura 5–14 Red de ﬂujo en una presa de gravedad ligeramente enterrada y con una pantalla impermeable y una base impermeable a una cierta 
profundidad. 
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h0 


h1 


p 
ρg 


Figura 5–15 Distribución de presiones bajo una presa 


debajo de la superﬁcie externa de la presa y por tanto no haya superﬁcie de escorrentía más que en el dren, donde se 
controla que la tierra no pase al mismo mediante rejillas y ﬁltros de tamaño de poro adecuado. 


Presas de gravedad 


En la ﬁgura 5–13 se muestra un ejemplo de red de ﬂujo bajo una presa de gravedad sencilla cuya base está al 
mismo nivel del terreno, que se supone homogéneo e isótropo. Como se ve en la ﬁgura, en estas condiciones el 
ﬂujo es simétrico y los potenciales hidráulicos van desde un valor máximo igual a h en la superﬁcie ﬁltrante aguas 
arriba a un valor mínimo aproximadamente nulo (suponiendo un caudal muy pequeño) aguas abajo. Las superﬁcies 
equipotenciales se muestran con líneas de trazos y las líneas de ﬂujo con líneas continuas. 


Puede mejorarse notablemente el comportamiento de esa presa frente a las ﬁltraciones mediante dos sencillas modi- 
caciones en el diseño. Por una parte se entierra parcialmente la presa, se forma que su parte inferior queda por debajo 
del nivel del terreno y por otra se coloca una pantalla impermeable en su parte inferior como se muestra en la ﬁgura 
5–12. Estas medidas hacen que el peligro de perforación (tubiﬁcación) se reduzca notablemente así como que las 
velocidades de salida aguas abajo se reduzcan como consecuencia del alargamiento de los recorridos de las líneas de 
ﬂujo, dando como resultado una disminución de las ﬁltraciones bajo la presa. 


5.2.10. 
Subpresiones 


La presencia de ﬁltraciones bajo una estructura como la de una presa introduce un elemento adicional a tener en 
cuenta en cuanto a su estabilidad, debido a la aparición de una serie de presiones en la base del elemento, denominadas 
subpresiones. Estas subpresiones producen por una parte un empuje y por otra parte un momento de vuelco a causa de 
la no homogeneidad de su distribución, debiendo ambos efectos ser tenidos en cuenta. 


Por ejemplo, en una presa en la que no se han tomado medidas para aminorar el efecto de estas subpresiones la 
distribución de presiones en su base puede ser, en la práctica, de la forma que se muestra en la ﬁgura 5–15, apro- 
ximadamente lineal. Esto hace que, a la hora de estimar el momento de vuelco sea necesario considerar, además de 
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los momentos debidos a las fuerzas actuando sobre las paredes verticales el efecto de las presiones sobre la base. 
Asimismo, si el suelo está por debajo del nivel freático puede ser necesario tomar medidas para evitar el efecto del 
empuje sobre éste. 


Para amortiguar el efecto de las subpresiones se instalan, cerca de la zona aguas arriba, drenajes verticales dentro de 
la presa desde la base de la misma a una galería de desagüe a un nivel superior al nivel aguas abajo. Este mismo sistema 
se utiliza en losas grandes, donde es preferible taladrar verticalmente el piso para aliviar el efecto de las subpresiones 
y evacuar el agua ﬁltrada. En otros casos se utilizan tirantes anclados bajo el piso para contrarrestar el efecto de las 
subpresiones. 
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